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PREFACE 


Reviews  of  experimental  techniques  are  necessary  from  time  to  time  due  to  the  rapid  development  of  the  hardware  and 
software  of  the  respective  methods.  Flow  visualization  is  an  important  diagnostic  tool  for  wind  tunnels,  water  tunnels,  and 
other  facilities,  since  it  provides  an  over- all  view  of  the  flow  in  a  whole  field,  sometimes  even  allowing  for  taking  quantitative 
data.  Flow  visualization  has  had  a  considerable  benefit  from  the  development  of  laser  optics  and  computational  methods. 
This  development  is  reviewed  together  with  a  presentation  of  the  physical  background  of  the  existing,  sometimes  “classical” 
methods  and  principles.  Very  essential  in  such  a  review  is  the  presentation  of  examples.  I  am  grateful  to  all  authors  and 
institutions  who  have  provided  sample  photographs  to  this  report  and  thereby  helped  to  facilitate  the  understanding  of  the 
particular  techniques.  I  am  also  thankful  to  Mrs  H.Miiller  who  has  typed  the  manuscript. 

Essen,  May  1987 
W.Merzkirch 


*  * 


* 


Les  materiels  et  logiciels  auxquels  font  appel  les  techniques  experimentales  evoluent  a  un  rythme  tel  qu’il  est 
indispensable  de  faire  periodiquement  le  point  de  la  situation.  Que  ce  soit  en  soufflerie  aerodynamique  ou  hydraulique,  ou 
toute  autre  installation  similaire,  la  visualisation  des  ecoulements  constitue  un  outil  essentiel  de  diagnostic  car  elle  procure 
une  vue  globale  d’un  ecoulement  complet  et  permet  meme  parfois  de  recueillir  des  donnees  quantitatives. 

Le  developpement  de  l’optique  laser  et  des  methodes  de  calcul  a  grandement  contribue  a  Famelioration  de  la 
visualisation  des  ecoulements.  L’etude  presentee  porte  sur  ce  developpement  et  recense  les  fondements  physiques  des 
methodes  et  principes  existants,  fussent-ils  “classiques”.  Elle  inclut  nombre  d’exemples  sans  lesquels  elle  ne  saurait  etre 
complete,  c’est  pourquoi  je  remercie  tous  les  auteurs  et  organismes  dont  les  photographies  iliustrent  cette  etude  et  facilitent 
la  comprehension  des  techniques  particulieres  auxquelles  elles  se  rapportent. 


Essen,  mai  1987 
W.Merzkirch 
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1  .  INTRODUCTION 


The  role  of  flow  visualization  in  experimental  fluid-mechanical  research  has  been 
appraised  many  times,  and  a  number  of  reviews  or  comprehensive  descriptions,  either  of 
the  whole  field  or  particular  applications,  are  available;  see,  e.g.,  [1.1-1. 8].  The  de¬ 
velopment  of  new  techniques  and  instruments  make  it  necessary  to  revise  the  visualiza¬ 
tion  methods  from  time  to  time.  The  majority  of  the  recent  technical  developments  in 
this  area  result  from  the  progress  made  in  laser  optics  and  the  wide-spread  use  of  com¬ 
puters  for  analyzing  experimental  data.  Laser-induced  fluorescence,  speckle  photography, 
image  processing,  and  optical  tomography  are  examples  of  such  modern  developments. 

The  application  of  novel  techniques  can  provide  new  insight  into  existing  problems  in 
fluid  mechanics.  As  an  example,  two  references  are  quoted  which  discuss  the  state  of  the 
art  of  research  in  turbulence:  "Flow  visualization  ...  appears  to  have  found  a  kind  of 
rebirth  in  recent  years.  ...  Flow  visualization  is  essential  not  only  in  obtaining  a 
clear  perception  of  the  flow  ...  but  also  in  planning  sophisticated  measurements.  ... 
Still  another  facet  of  the  second  revolution  (i.e.:  in  turbulence  research)  in  which 
flow  visualization  plays  the  central  role  is  image  processing"  [1.9] ,  and:  "It  is  iron¬ 
ical  that  coherent  structures  ...,  the  Brown-Roshko  vortices,  and  longitudinal  vortices 
in  the  sublayer  of  a  turbulent  boundary  layer  were  all  found  by  the  most  primitive  of 
experimental  methods,  flow  visualization"  [1.10]. 

The  latter  reference  brings  up  the  question  on  what  kind  of  information  is  availa¬ 
ble  from  a  flow  photograph,  and  how  does  this  compare  with  data  from  different  experi¬ 
mental  sources.  Again,  this  problem  has  been  the  subject  of  extensive  discussions  and 
reviews,  and  the  situation  may  well  be  illustrated  by  quoting  Hussain  [1.9]:  "Flow  vi¬ 
sualization  presents  excessive  information  but  very  little  hard  data,  and  hot-wires  or 
LDAs  give  some  hard  data  but  very  little  flow  physics". 

This  report  reviews  the  methods  of  flow  visualization  with  particular  emphasis  on 
their  application  in  wind  tunnels,  water  tunnels,  and  experiments  related  to  propulsion 
research.  The  definition  of  a  method  to  be  included  is  that  it  provides  a  visual  pat¬ 
tern  of  the  flow  in  form  of  a  whole-field  record,  e.g.,  a  photograph.  It  will  be  shown 
that,  in  contrast  to  the  references  quoted  previously,  some  of  the  newly  developed  tech¬ 
niques  are  not  really  primitive  and  involve  complicated  optical  instruments,  and  they 
may  even  provide  "hard  data"  on  the  flow  field. 

The  subject  of  many  flow  studies  is  the  interaction  of  a  solid  body  with  a  fluid 
flow.  Therefore,  the  discussion  of  methods  starts  with  those  that  visualize  the  flow 
pattern  on  a  solid  surface.  The  following  chapters  deal  with  methods  visualizing  the 
whole  field,  either  by  introducing  a  foreign  substance  (smoke,  dye)  or  by  making  use  of 
refractive  index  changes  in  the  flowing  fluid  (optical  methods) .  Emphasis  is  placed  on 
the  discussion  of  the  physical  background  of  a  method  rather  than  on  technical  details. 
The  latter  may  be  found,  e.g.,  in  laboratory  reports  some  of  which  are  referred  to  in 
the  cited  bibliography  of  each  chapter. 
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2.  SURFACE  FLOW  PATTERNS 


The  rates  of  shear  forces,  pressure  forces,  and  heating  loads  excerted  by  a  fluid 
flow  to  the  wall  of  a  solid  body  are  subject  of  many  studies.  Visualization  of  the  flow 
pattern  very  close  to  or  at  the  body  surface  can  be  helpful  for  such  studies.  Many  meth¬ 
ods  of  surface  flow  visualization  are  making  use  of  a  coating  of  the  surface  with  a  ma¬ 
terial  that  interacts  with  the  flow  and  thereby  produces  a  visible  information.  Depend¬ 
ing  on  whether  the  interaction  is  mechanical,  chemical,  or  thermal,  information  on  the 
following  quantities  might  be  derived: 

wall  pressure,  skin  friction,  or  simply  flow  direction; 

-  mass  transfer  to  the  wall; 

-  wall  heat  transfer  or  wall  temperature. 

In  the  following  paragraphs,  the  methods  for  visualizing  surface  flow  patterns  will 
be  discussed  in  the  indicated  order  of  sequence.  The  physical  processes  on  which  these 
methods  depend:  transfer  of  momentum,  mass,  and  heat,  are  strongly  influenced  by  the 
circumstances  whether  the  flow  is  laminar  or  turbulent,  attached  or  separated .  Therefore, 
in  the  simplest  way  of  their  application,  the  methods  just  serve  to  detect  the  loca¬ 
tions  of  transition,  separation,  or  attachment. 


2.1  OIL  FILM  VISUALIZATION 


Since  several  decades  of  wind  tunnel  practice  oil  film  visualization  has  become  a  stand¬ 
ard  wind  tunnel  technique.  A  number  of  reviews  with  many  practical  advices  are  available, 
e.g.,  Maltby  [2.1]  or  Settles  and  Teng  [2.2].  The  surface  of  interest  is  coated  with  a 
thin  layer  of  oil  in  which  a  finely  powdered  pigment  is  solved.  When  the  wind  tunnel  is 
turned  on,  the  oil  is  carried  away  with  the  air  stream,  and  the  dry  pigment  remains  on 
the  surface  where  it  forms  a  streaky  pattern  indicating  the  direction  of  the  flow  close 
to  the  surface.  Particular  phenomena  in  the  remaining  pattern  are  attributed  to  flow 
separation  from  the  surface  or  flow  reattachment.  A  test  model  to  which  the  mixture  of 
oil  and  pigment  was  applied  can  be  taken  out  of  the  test  section,  after  the  tunnel  has 
been  stopped,  and  the  pattern  on  the  model  surface  can  be  inspected  outside  the  wind  tun¬ 
nel  (Fig.  2.1). 

Many  instructions  have  been  given  in  the  literature  on  how  to  prepare  an  oil/pig¬ 
ment  mixture  appropriate  for  a  specific  test  condition;  e.g.,  [2.2]  to  [2.4].  A  large 

amount  of  practical  experience  is  collected  and  described  in  unpublished  laboratory  re¬ 
ports  or  manuals.  The  reader  might  find  as  many  different  recipees  as  he  consults  re¬ 
ports.  The  object  is  to  prepare  a  mixture  of  such  a  consistency  that  it  will  run  easily 
under  the  given  test  conditions  and  leave  behind  the  desired  streaks  of  the  pigment. 
Ideally,  the  mixture  should  not  begin  to  run  until  the  selected  wind  speed  is  reached, 
and  after  an  appropriate  time  of  running,  the  pattern  should  be  sufficiently  dry  to  be 
unaffected  by  the  unsteady  air  flow  when  the  tunnel  is  stopped. 
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Fig.  2.1:  Oil  film  pattern  on  the  upper  surface  of  the  right  wing  of  an  orbiter  model 
at  M  =  7.4  wind  tunnel  flow  and  18°  angle  of  attack.  Attached  flow  in  the 
front  part  of  the  surface.  Contour  lines  indicate  the  formation  of  a  conical 
trailing  vortex  and  its  interaction  with  the  separated  flow  on  the  rear  side 
of  the  wing  (NASA  Ames  Research  Center,  Moffett  Field,  CA) . 

From  the  latter  it  follows,  that  the  lower  the  air  velocity,  the  lower  should  be 
the  viscosity  of  the  oil.  Liquids,  that  have  been  used  als  solvents  or  carriers  for  the 
pigment,  are,  in  ascending  order  of  viscosity:  alcohol,  kerosene,  light  Diesel  oil, 
light  transformer  oil.  Viscosity  and  surface  tension  of  these  liquids  can  be  adjusted  by 
additives.  Two  special  situations  for  modern  wind  tunnel  testing  should  be  mentioned: 

The  low  freestream  pressures  associated  with  high-speed  wind  tunnels  are  unfavorable  for 
using  standard  oils,  because  these  oils  are  characterized  by  high  vapor  pressures;  then, 
vacuum-pump  oil  is  the  right  choice.  Propane  (C3Hg) ,  which  has  a  low  melting  point,  has 
been  employed  in  a  cryogenic  wind  tunnel  [2.5]  where  the  mentioned  standard  oils  would 
freeze . 

The  pigment  in  the  mixture  should  provide  a  clear  pattern  against  the  model  sur¬ 
face  as  the  background.  A  white  powder  like  titanium  dioxide  (Ti02)  or  china  clay  can 
be  used  on  a  dark  model  (Fig.  2.2) .  Lampblack,  a  fine  powder  that  mixes  well  with  oil 
is  most  suitable  on  a  light  model  surface;  see,  e.g.,  Meznarsic  and  Gross  [2.6],  or 
Sparrow  and  Comb  [2.7].  The  dispersion  of  the  pigment  in  oil  may  be  improved  by  adding 
a  few  drops  of  oleic  acid,  see  Keener  [2.8] .  Fluorescent  pigments,  when  illuminated 
with  ultraviolet  light,  can  provide  a  very  brillant  surface  pattern;  only  the  surface 
pattern  to  be  visualized  is  seen,  whereas  model  and  background  remain  invisible.  This  is 
of  particular  interest  if  the  model  geometry  is  complex;  see,  e.g.,  Gessner  and  Chan 
[2.9]  . 


The  pigment  pattern  remaining  on  the  model  surface  can  be  photographed  for  the  pur¬ 
pose  of  recording  it.  If  the  model  surface  is  curved,  a  problem  comes  up  in  relating  the 
plane  record  of  the  curved  surface  pattern  to  a  surface  coordinate  system.  As  a  solution 
to  this  problem,  Sparrow  and  Comb  [2.7]  cover  the  model  surface  with  white,  plasticized, 
self-adhering  contact  paper,  and  they  apply  the  oil-pigment  (lampblack)  mixture  to  the 
paper.  After  the  visible  pattern  is  formed,  the  paper  is  separated  from  the  wall  and 
laid  flat  for  being  photographed.  For  the  same  purpose,  i.e.  for  eliminating  camera 
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Fig.  2.2:  Oil  film  pattern  on  the  plane  side  walls  of  a  cascade  wind  tunnel  at  Mach 

number  M  =  0.4  and  with  weak  blowing  at  the  front  edge  of  the  blades  (DFVLR- 
Institut  fur  Antriebstechnik ,  Koln,  F.R. Germany) . 


parallax  when  photographing  the  surface  pattern,  Settles  and  Teng  [2.6]  apply  transpar¬ 
ent  adhesive  tapes  to  the  model  after  the  wind  tunnel  experiment.  The  very  thin  coating 
of  pigment  remaining  on  the  surface  is  lifted  off  and  preserved  for  taking  a  record. 

The  interpretation  of  the  observed  pattern  is,  of  course,  strongly  related  to  the 
question  on  how  the  pattern  is  formed.  Several  explanations  have  been  given;  but  the 
problem  has  never  been  investigated  systematically.  It  has  been  said  that  the  fine  par¬ 
ticles  of  the  pigment  coagulate,  and  that  the  striations  are  caused  by  the  wake  behind 
such  coagulations  in  the  flowing  oil  film.  Hornung  [2.10]  observed  that  the  flowing  oil 
film  becomes  unstable,  breaks  up  to  form  droplets  which,  due  to  their  large  resistance, 
shoot  at  high  speed  across  the  oil  film,  thus  producing  the  striations.  Murai  et  al. 
[2.11]  verified  longitudinal  instabilities  in  the  thin  oil  film  which  could  be  respon¬ 
sible  for  the  observed  pattern. 

Information  on  the  accuracy  of  the  oil  film  technique  can  be  obtained  if  one  ana¬ 
lyzes  the  flow  pattern  in  the  oil  film.  The  only  available  analysis  has  been  developed 
by  Squire  [2.12].  He  investigates  the  problem  of  how  the  air  flow  close  to  the  wall  is 
disturbed  due  to  the  presence  of  the  oil  film.  The  boundary  conditions  for  the  air  flow 
now  have  to  be  fulfilled  at  the  oil/air  interface  where  the  velocity  of  the  two  fluids 
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is  the  same  and  not  zero.  The  velocity  u  at  this  interface  is 


u (y=h) 


**  '  y=h 


(2.1) 


where  x  is  the  streamwise  coordinate  in  the  plane  of  the  wall,  y  is  the  coordinate  nor¬ 
mal  to  the  wall,  h  is  the  oil  film  thickness,  [i  viscosity,  and  p  pressure.  Since  the  ra¬ 
tio  ^a:Lr/^oil  is  of  the  order  1CT2  to  10-4,  u(y=h)  is  normally  a  very  small  quantity, 
i.e.,  the  change  in  the  conditions  for  the  air  flow  can  be  neglected,  with  the  exception 
of  one  situation,  namely,  when  the  pressure  change  in  streamwise  direction,  dp/dx,  be¬ 
comes  very  large.  This,  however,  is  the  case  at  the  location  of  flow  separation  or  at¬ 
tachment,  particularly  if  the  separation  is  enforced  by  a  shock  wave.  Squire’s  analysis 
gives  no  information  on  how  big  the  error  in  the  measurement  will  be,  e.g.  in  terms  of 
the  difference  between  the  indicated  and  the  real  location  of  separation.  But,  since 
many  oil  film  experiments  are  used  for  visualizing  the  positions  of  flow  separation 
and  attachment.  Squire's  results  are  a  strong  warning  that  such  measurements  must  be 
carefully  interpreted. 

The  applications  of  the  oil  film  technique  for  aerodynamic  testing  are  so  numer¬ 
ous,  that  only  a  few  references  can  be  listed,  which  are  representative  for  many  others. 
Fow  low-speed  flows,  see  George  [2.13];  high-speed  configurations:  Reding  and  Ericsson 
[2.14];  separation  by  shock  impingement:  Settles  et  al.  [2.15]. 


2.2  OIL  FILM  SKIN  FRICTION  MEASUREMENTS 

In  various  publications,  e.g.,  [2.16-2.18],  Tanner  demonstrates  that  it  is  possible 

to  relate  the  height  of  the  oil  film,  h,  to  the  skin  friction  or  the  wall  shear  stress 
of  the  air  flow,  x  .  Suppose  that  the  oil  has  been  deposited  on  the  body  surface  in  form 
of  a  small  drop  or  along  a  straight  line  normal  to  the  main  flow  direction.  The  air  flow 
is  started  at  time  t  =  0,  the  oil  drop  develops  into  a  thin  film  whose  thickness  h  var¬ 
ies  with  position  x  and  time  t  according  to  (see  Fig.  2.3) 


incident  I  ight 


/ 


/reflected  light 


Fig.  2.3: 

Measurement  of  the  oil 
film  thickness  h  at  time 
t  >  0  by  means  of  optical 
interference  of  two  re¬ 
flected  light  rays  (ac¬ 
cording  to  [2.16]). 
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x 


h(x,t) 


T, 


w 


dx 


1/2 


(2.2) 


With  known  value  of  the  viscosity  of  the  oil,  it  is  possible  to  determine  the  skin 

friction  T  ,  once  the  oil  film  thickness  can  be  measured.  Since  the  film  is  very  thin, 
w 

Tanner  measures  h  by  interferometric  means  (Fig.  2.3).  Interference  is  generated  between 
light  rays  reflecting  from  the  surface  of  the  oil  film  and  from  the  body  surface,  res¬ 
pectively.  With  the  surface  being  plane  and  flat,  it  is  possible  to  take  an  interfero- 
gram  of  the  whole  surface  for  one  particular  instant  of  time,  t^  >  0  (Fig.  2.4) .  The  vis¬ 
ible  interference  fringes  are  curves  of  equal  height  h,  and  the  photograph  includes  in¬ 
formation  on  the  distribution  of  the  wall  shear  stress,  t  ,  in  the  observed  surface. 

w 


Fig.  2.4:  Curves  of  equal  height  visualized  by  optical  interferometry  in  the  oil  film 


traces  that  have  developed  from  a  number  of  oil  drops  on  a  plane  surface 
(planview) .  Main  flow  direction  of  the  air  is  from  left  to  right.  The  air  is 
forced  to  flow  around  a  vertical  wedge.  The  oil  is  collected  in  the  center 
of  two  vortices  behind  the  wedge,  as  indicated  by  the  dark  spots  (L.H.  Tanner, 
Brighton  Polytechnic,  England) . 


For  the  application  to  curved  body  surfaces.  Tanner  describes  a  modification,  which 
no  longer  delivers  information  on  a  whole  field,  but  measures  the  oil  film  thickness  at 
a  point  and  as  a  function  of  time.  Two  focused  laser  beams  that  are  reflected  from  the 
solid  surface  and  from  the  oil  surface,  respectively,  interfere  with  oneanother,  as  in¬ 
dicated  in  Fig.  2.3.  A  photodiode  records  the  interference,  i.e. ,  alternate  bright  and 
dark  signals,  according  to  the  variation  with  time  of  the  film  thickness  h  at  a  given 
position  x.  Wall  shear  stress  measurements  have  been  performed  with  this  "skin  friction 
probe"  in  supersonic  flow,  [2.19,2.20]. 

2.3  WALL  TUFTS 

A  relatively  simple  means  for  obtaining  an  idea  on  the  direction  of  flow  close  to  a  sol¬ 
id  wall  is  to  attach  one  end  of  short  tufts  to  the  body  surface.  In  laminar,  attached 
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Fig.  2.5:  Fluorescent  minitufts  on  the  model  of  a  transport  airplane  in  a  wind  tunnel 
(J.P.  Crowder,  Boeing  Aero  Laboratories) . 


flow,  these  tufts  may  well  indicate  the  local  flow  direction.  When  the  flow  becomes  un¬ 
steady  or  turbulent,  the  tufts  perform  a  certain  unsteady  motion,  and  this  may  be  taken 
as  an  indication  that  the  wall  boundary  layer  has  become  turbulent.  A  more  violent  mo¬ 
tion  of  the  tufts,  or  a  tendency  to  lift  from  the  surface,  may  indicate  a  separated  flow 
regime. 

The  choice  of  tuft  size  and  material  depends  on  the  flow  conditions  and  the  size  of 
the  model  to  be  tested.  Tufts  of  ordinary  yarn  and  several  centimeters  long  have  been 
used  on  full-scale  car  models  in  a  wind  tunnel,  and  on  full-scale  airplanes  in  flight. 

It  is  obvious  that  tufts  of  this  size  and  the  corrugation  caused  by  the  glue  or  the 
device  fixing  the  tufts  affect  the  surface  conditions  of  the  flow.  Furthermore,  it  has 
been  observed  that,  due  to  instabilities  induced  at  particular  wind  speeds,  the  tufts 
may  perform  a  self-excited  flagging  motion,  which  obscures  the  indication  of  flow  sepa¬ 
ration.  For  in-flight  tests,  Crowder  and  Robertson  [2.21]  use  light,  rigid,  narrow,  co¬ 
nical  elements,  which  are  attached  to  the  wind  surface  with  a  short  string  at  their 
apex.  In  contrast  to  the  tufts,  these  elements,  several  centimeters  long,  behave  stable 
in  attached  flow. 

In  order  to  minimize  interference  between  tufts  and  flow,  Crowder  [2.22]  has  de¬ 
veloped  "minitufts"  made  from  thin  nylon  monofilament,  with  a  diameter  of  about  20  urn. 

In  order  to  enhance  the  visibility,  the  nylon  is  treated  with  fluorescent  dye,  and  the 
tufts  are  observed  or  photographed  with  UV  illumination  (Fig.  2.5).  The  surface  flow 
pattern  on  rotating  propellers  [2.23]  as  well  on  models  in  a  water  tunnel  [2.24]  has 
been  successfully  visualized  with  these  minitufts. 


2.4  SURFACE  FLOW  VISUALIZATION  DEPENDING  ON  MASS  TRANSFER 

Sublimation  evaporation  of  a  coating  material  into  the  adjacent  air  flow  is  a 
means  for  visualizing  the  mass  transfer  between  the  solid  surface  and  the  flowing  gas 
Use  of  this  effect  is  made  by  the  surface  techniques  employing  coatings  of  either 
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Fig.  2.6:  Visualization  of  the  location  of  laminar-to-turbulent  transition  on  the  low 
pressure  surface  of  a  turbine  blade  model.  Visualization  is  by  means  of  azo¬ 
benzene  whose  sublimation  is  more  rapid  in  the  turbulent  boundary  layer,  this 
region  now  appearing  dark  (DFVLR-Institut  fur  Antriebstechnik ,  Koln,  F.R. Ger¬ 
many)  . 


naphtalene  {see  several  publications  by  Sparrow  et  al.,  e.g.,  [2.25]),  or  azobenzene 

[2.26].  The  rates  of  sublimation  or  mass  transfer  are  different  depending  on  whether  the 
flow  is  laminar  or  turbulent,  so  that  the  location  of  transition  can  be  made  visible 
(Fig.  2.6).  The  naphtalene  technique  can  even  provide  some  degree  of  quantitative  infor¬ 
mation:  The  surface  is  coated  with  a  layer  of  dense  naphtalene,  the  outer  surface  of  the 
layer  being  smooth.  Naphtalene  sublimates  into  the  outer  air  flow,  and  the  resulting 
profile  of  the  originally  smooth  surfaces  reflects  low  and  high  mass  transfer  rates.  A 
measurement  of  the  naphtalene  profile  allows  for  obtaining  values  of  the  local  mass 
transfer  coefficient. 

A  very  distinct  surface  visualization  is  achieved  with  methods  in  which  the  mass 
transfer  results  in  a  chemical  reaction  causing  a  color  change  on  the  body  surface.  For 
the  technique  described  by  Kottke  [2.27,2.28],  the  wall  is  coated  with  a  thin  wet  layer 
(filter  paper  or  gel)  containing  an  aqueous  solution  of  MnCl2  and  H2°2 *  A  reacting  gas, 
NH^ ,  is  added  in  form  of  a  short  pulse  and  at  low  concentration  to  the  main  airstream. 
NH^  is  adsorbed,  according  to  the  local  partial  concentration  differences,  by  the  wet 
layer,  and  a  reaction  takes  place  in  which  MnO^  is  formed  as  an  end  product.  While  the 
original  coating  is  bright,  MnC>2  is  dark,  so  that  the  observed  color  intensity  of  MnC>2 
is  a  measure  of  the  local  mass  transfer  rate  (Fig.  2.7).  Kottke  performed  quantitative 
measurements  by  photometric  means.  He  observerd  that  the  measured  position  of  maximum 
mass  transfer  in  a  reattaching  flow  is  not  consistent  with  the  expected  line  of  flow 
attachment.  This  is  just  another  indication  of  the  difficulties  that  arise  in  inter¬ 
preting  surface  patterns  produced  by  separating  flows. 


2.5  VISUALIZATION  OF  WALL  TEMPERATURE 


The  derivation  of  aerodynamic  heat  transfer  coefficients  requires  the  recording  of 
the  time  history  of  a  number  of  isotherms  on  the  surface  under  study.  Such  isotherms 
can  be  visualized  by  means  of  surface  coatings  whose  pattern  or  color  reacts  to  changes 
of  the  surface  temperature.  Two  coatings  of  this  kind  will  be  discussed  here:  tempera- 
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Fig.  2.7:  Visualization  of  mass  transfer  to  the  surface  of  a  car  model  in  a  wind  tun¬ 
nel.  White  areas  indicate  separated  flow  regimes  with  low  mass  transfer 
rate  (V.  Kottke ,  Universitat  Stuttgart,  F. R. Germany ) . 


ture  sensitive  paints  and  liquid  crystals,  in  addition  to  a  method  that  does  not  rely  on 
a  coating:  infrared  photography.  The  three  methods  have  different  ranges  of  application 
and  sensitivity.  The  paints  are  used  at  temperatures  of  several  hundred  Kelvin  in  high- 
enthalpy  flow  facilities.  Liquid  crystals  are  sensitive  in  a  range  of  only  a  few  Kelvin, 
and  their  application  is  at  temperatures  not  too  far  from  room  temperature.  Applicabili¬ 
ty  and  sensitivity  of  the  infrared  camera  is  somewhere  between  the  two  former  methods. 

Temperature-sensitive  paints  have  been  used  since  many  years  of  high-speed  wind 
tunnel  testing,  [2.29,2.30].  These  paints  consist  of  several  components  each  of  which 
undergoes  a  visible  change  of  its  internal  structure  or  phase  at  a  specific  temperature. 
The  edge  between  two  colors  is  a  curve  of  constant  surface  temperature.  If  the  model  is 
made  of  a  material  with  low  thermal  conductivity,  the  indicated  temperature  may  be  taken 
as  the  adiabatic  wall  temperature.  The  edges  between  colors  move  on  the  body  surface  in 
downstream  direction  with  increasing  time,  until  the  whole  body  surface  might  have 
reached  a  uniform  temperature.  A  number  of  different  paints  are  commercially  available. 

The  edge  between  two  colors  of  the  paint  is  to  a  certain  degree  diffusive,  so  re¬ 
ducing  the  local  resolution  of  the  mapping.  More  accurate  in  this  respect  is  a  one-com¬ 
ponent  paint  {"phase  change  paint")  which  changes  from  solid  to  liquid  at  a  precise  sur¬ 
face  temperature  [2.31-2.34].  Wax  can  be  the  major  component  of  such  coating.  The  ob¬ 
served  melt  line  is  an  isotherm.  Keyes  [2.35]  reports  that  the  coating  can  be  blown  off 
due  to  high  shear  rates  in  the  flow,  and  this  phenomena  could  be  misinterpreted  to  be  a 
phase  change.  In  the  cases  of  both  temperature  sensitive  paint  and  phase  change  paint, 
the  process  causing  the  visible  color  change  is  not  reversible,  so  that  this  coating  of 
the  surface  can  be  used  only  once. 

A  reversible  color  change  as  a  function  of  the  surface  temperature  can  be  obtained 
with  a  coating  of  liquid  crystals.  Klein  [2.36]  gave  a  first  description  of  the  use  of 
liquid  crystals  for  the  study  of  aerodynamic  surface  heating.  Temperature  changes  in 
such  crystals  disturb  the  intermolecular  forces  and  generate  a  shift  in  the  molecular 
structure,  thus  causing  a  shift  in  wavelength  of  the  light  scattered  from  the  crystals. 
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Fig.  2.8:  Surface  temperature  mapping  by  liquid  crystals  in  a  hypersonic  wind  tunnel. 

Fins  on  a  cylindrical  body  whose  axis  is  in  the  flow  direction  produce  ob¬ 
lique  shock  waves*.  The  intersection  of  these  shock  waves  with  the  cylindri¬ 
cal  surface  is  visible  due  to  the  temperature  rise  through  the  shock  (H, 
Scholer,  DFVLR-Institut  fur  Experimentelle  Stromungsmechanik ,  Gottingen, 

F. R. Germany) . 


An  extensive  literature  is  available  on  liquid  crystals  and  their  application  in  differ¬ 
ent  fields.  Among  the  various  types,  cholesteric  liquid  crystals  have  been  proven  to  be 
most  suitable  for  studies  in  fluid  mechanics  and  convective  heat  transfer.  Unlike  tem¬ 
perature  sensitive  paints,  cholesteric  liquid  crystals  react  to  changes  of  temperatures 
only  within  a  few  degrees  Kelvin.  The  crystals  are  commercially  available,  for  a  great 
range  of  both  bandwidth  and  mean  temperature  at  which  they  are  sensitive.  Within  the 
bandwidth  of  their  sensitivity,  up  to  ten  color  hues  may  be  discriminated  by  the  eye. 

The  visible  edge  between  two  colors  is  a  curve  of  constant  temperature  (Fig.  2.8).  The 
value  of  the  temperature  must  be  found  by  calibration  for  constant  viewing  direction. 

A  thin  sheet  of  liquid  crystals  on  a  solid  surface  might  be  affected  by  shear 
forces  or,  if  the  coating  is  used  in  water,  by  contamination.  In  order  to  protect  the 
crystals,  they  can  be  made  available  in  the  form  of  sheets  consisting  of  the  liquid 
crystals  laid  on  plastic  material  and  covered  by  a  transparent  protecting  layer  [2.37, 
2.38],  or  they  are  "microencapsulated",  i.e.,  enclosed  in  plastic  spheres  of  10  to  50  urn 
in  diameter  [2.39].  Dispersing  the  encapsulated  liquid  crystals  in  water  gives  an  emul¬ 
sion  which  can  be  brushed  or  sprayed  on  the  model  surface.  Cholesteric  liquid  crystals 
have  been  applied  to  the  flow  in  wind  tunnels  [2.39,2.40]  and  in  water  tunnels  [2.41, 
2.42]  . 


The  surface  of  a  solid  body  emits  infrared  (IR)  radiation  whose  intensity  is  a 
function  of  the  surface  temperature.  This  radiation  can  be  received  with  an  IR  camera, 
which,  in  connection  with  appropriate  computer  hardware  and  software,  converts  the 
signal  into  a  visual  pattern.  Several  successful  attempts  have  been  made  for  mapping 
the  surface  temperature  distribution  of  models  in  wind  tunnels  or  in  other  flow  en¬ 
vironments  by  means  of  an  IR  camera;  see,  e.g.  [2.43,2.44].  Most  of  the  reported  experi¬ 
ments  have  been  performed  with  an  IR  camera  which  is  commercially  available?  it  employs 
a  detector  sensitive  to  the  wavelength  range  between  3  and  5.8  fim.  The  detector,  which 
must  be  cooled  by  liquid  nitrogen,  produces  an  electric  signal  proportional  to  the  re¬ 
ceived  total  energy.  The  total  radiant  energy  emitted  by  the  surface  under  investigation 
is,  according  to  the  Stef an-Boltzmann  law  for  a  gray  body,  proportional  to  the  4th  power 
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Fig.  2.9: 

Infrared  (IR)  camera  display  of  the 
temperature  distribution  on  a  surface 
on  which  a  gas  jet  is  impinging.  The 
visible  pattern  are  surface  tempera¬ 
ture  contours.  Mean  jet  flow  veloci¬ 
ty  is  45  m/s  (R.H.  Page,  C.  Ostowari, 
Texas  A&M  University) . 


of  the  absolute  temperature.  Gauffre  and  Fontanella  [2.45]  give  an  overview  on  IR  cam¬ 
eras  used  for  aerodynamic  experiments,  while  Boylan  et  al.  [2.43]  describe  ways  of  ca¬ 
libration  of  the  system.  The  surface  temperature  pattern  can  be  displayed  in  real  time 
on  a  monitor  (Fig.  2.9). 

A  way  for  determining  the  aerodynamic  heat  transfer  coefficient  h  has  been  devised 
by  Jones  and  Hunt  [2.31],  For  this  purpose,  it  is  necessary  to  record  the  isotherms  as 
a  function  of  time,  e.g.,  by  taking  a  movie  film.  The  method  is  based  on  the  transient, 
one-dimensional  analysis  of  heat  flow  from  the  model  surface  into  the  body.  The  differ¬ 
ential  equation 


3  T  3  2  T 

__  =  a*  ~ — r 

3  t  3  y2 


(2.3) 


(with  T  =  absolute  temperature,  t  =  time,  y  =  coordinate  normal  to  the  model  surface, 
a  =  k/pCp  =  thermal  diffusivity  of  the  model  material,  k  =  thermal  conductivity,  p  = 
density  of  model  material,  and  c^  its  specific  heat)  is  solved  for  the  following  ini¬ 
tial  and  boundary  conditions : 

T(y,0)  =  T^  is  the  initial  temperature  of  the  model.  In  practice,  t =  0  is  the  instant 
at  which  the  model  is  injected  into  the  wind  tunnel. 

T(°°,t)  =  T^:  thermal  diffusion  within  the  model  material  is  slow. 

_  3  -l^t).  _  h.  (Taw  -  KO/t)) ,  with  Taw  =  adiabatic  wall  temperature,  h  =  aerodynamic  heat 

transfer  coefficient.  This  condition  describes  the  development  of  the  surface  temper¬ 
ature  (at  y  =  0)  due  to  aerodynamic  heating. 

With  a  number  of  additional  simplifying  assumptions ,  Jones  and  Hunt  develop  the  follow¬ 
ing  solution  for  the  heat  flow  equation  (2.3): 


T*  =  1  -  exp ( 32 ) • erfc  3 


(2.4) 


with  the  abbreviations 
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Tind  temPerature  °f  an  isotherm  indicated  by  the  specific  coating.  The  determina¬ 

tion  of  h  is  as  follows:  With  given  Tind  for  a  specific  indicator  one  determines  the 
value  of  T*  and,  from  a  plot  T*  versus  3  [2.31],  the  corresponding  value  of  3.  With 
known  values  k,  a  and  measured  time  t  (for  the  isotherm  to  appear  at  a  specific  position 
on  the  surface,  as  measured  from  the  movie  film) ,  one  determines  h  =  k* 3/ /at .  The  time 
required  for  the  isotherms  to  move  over  the  model  surface  should  be  large  compared  with 
the  time  necessary  for  establishing  a  steady  flow  around  the  model,  and  it  must  be  short 
compared  with  the  thermal  diffusion  time  of  the  model  material.  Segletes  [2.32]  points 
out  that  the  one-dimensional  heat  transfer  model  used  in  the  above  analysis  may,  in  some 
cases,  lead  to  large  errors.  He  presents  correction  factors  that  improve  the  results  of 
the  1  — D  analysis  for  particular  corner  flows. 
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3.  DIRECT  INJECTION  METHODS 


The  injection  of  a  foreign,  visible  material  is  a  traditional  means  of  flow  visualization. 
Most  commonly,  dyes  are  used  in  water,  and  smoke  in  air  flows.  The  material  is  released 
from  a  particular  position,  thus  marking  a  specific  line  in  the  flow  downstream  of  that 
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position?  or,  the  material  is  filling  a  complete  regime  in  the  flow,  and  it  marks  the 
boundary  and  therefore  the  geometric  extension  of  such  a  regime,  provided  that  there  is 
little  mass  exchange  across  the  boundary.  The  latter  is  of  interest  for  visualizing  sep¬ 
arated  flow  regimes. 

The  particles  of  the  injected  material  are  required  to  be  small  so  that  they  can 
well  follow  the  flow.  However,  the  motion  of  individual  particles  is  not  studied  here. 

The  number  density  of  the  particles  in  the  fluid  is  so  high,  that  the  individual  par¬ 
ticles  lose  their  identity  in  the  observations,  and  instead,  one  observes  a  continuous 
ensemble  of  the  "marked"  fluid,  e.g.  in  the  form  of  the  aforementioned  lines.  In  steady 
flow  the  line  is  a  streamline,  and  the  interpretation  of  a  pattern  of  such  lines  is 
straightforward.  For  an  unsteady  flow  field,  the  marked  contours  are  named  "streaklines", 
defined  as  the  instantaneous  locus  of  all  fluid  elements  having  passed  a  specific  posi¬ 
tion,  i.e.,  the  position  where  the  dye  or  smoke  is  injected.  Difficulties  can  arise  in 
the  interpretation  of  the  streakline  patterns  of  an  unsteady  flow;  see,  e.g.,  [3.1]. 


3 . 1  SMOKE  FLOW  VISUALIZATION 


The  visualization  of  the  flow  in  wind  tunnels  by  means  of  smoke  is  now  a  standard 
experimental  tool  for  the  work  in  these  facilities.  The  advancement  of  this  technique 
is  closely  related  to  the  development  and  the  history  of  wind  tunnels  [3.2].  An  essen¬ 
tial  portion  of  the  progress  and  refinement  of  the  smoke  technique  is  due  to  the  work 
of  Brown  [3.3]  at  the  University  of  Notre  Dame.  He  systematically  developed  the  genera¬ 
tion  of  appropriate  smokes  as  well  as  the  performance  of  suitable  wind  tunnels,  later 
referred  to  as  "smoke  tunnels".  Reviews  of  this  visualization  method  have  been  given 
by  Maltby  and  Keating  [3.4]  and  by  Mueller  [3.5]. 

The  term  "smoke"  is  used  here  in  a  wide  sense,  not  only  restricted  to  combustion 
products,  and  the  discussion  will  include  steam,  vapor,  mist,  and  aerosols.  From  the 
usual  requirement  that  one  wants  a  tracer  material  to  fulfill,  namely 

-  being  neutrally  buoyant, 

-  being  non-toxic,  and 

-  having  low  mixing  rates  with  the  main  fluid, 

smoke  fulfills  none  of  these.  The  use  of  smoke  is  therefore  a  compromise,  and  no  better 
technical  solution  has  been  found  yet. 

The  basic  types  of  producing  smoke  are:  burning  or  smoldering  tobacco,  wood  or 
straw;  vaporizing  mineral  oils?  producing  mist  as  the  result  of  the  reaction  of  various 
chemical  substances?  and  condensing  steam  to  form  a  visible  fog.  Except  of  the  latter, 
all  these  substances  are  toxic  to  some  degree.  The  density  of  these  tracer  materials  is 
much  larger  than  the  density  of  air,  but  since  the  particle  size  of  the  tracers  is  very 
small,  mostly  below  1  mu,  sedimentation  effects  are  minimized.  Diffusion  of  the  various 
types  of  smoke  into  the  ambient  air  is  much  higher  than  diffusion  of  dye  into  water,  so 
that  streaklines  can  be  identified  over  a  reasonable  distance  only  in  laminar  flow  (Fig 
3.1). 

A  number  of  possible  smokes  and  ways  of  their  generation  (often  called  smoke-gene¬ 
rators)  ,  together  with  relevant  references  are  listed  in  Table  3.1.  Cigarette  smoke  has 
the  finest  tracer  particles  among  all  smokes.  Kerosene  mist  probably  is  the  most  fre¬ 
quently  used  "smoke"  for  wind  tunnel  studies.  The  generation  of  TiC>2  mist  from  TiCl4 
droplets  deposited  on  the  surface  of  a  test  model  is  not  a  direct  injection  method, 
but  rather  a  method  of  chemical  controlled  production  of  the  trace  material  (see  also 
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a  =  6°  a  =  14° 

Fig.  3.1:  Smoke-line  visualization  of  the  flow  over  an  airfoil  at  Re=4o.ooo  and 
various  angles  of  attack,  a  (T.J.  Mueller,  University  of  Notre  Dame, 
Indiana,  USA) . 


Table  3 . 1 


Smokes  and  methods  of  their  generation 


cigarette,  smoke 
burning  pine  wood 
kerosene  mist 


Ti02  mist 
water  fog 


properties  see  [3.6] 
smoke  generator  see  [3.8] 

Preston-Sweeting  smoke  generator  [3.], 
see  also  [3.4]  [3.9]  ? 

University  of  Notre  Dame  smoke  generator  [3.5] 
smoke  wire,  see  chapter  5 

generator  based  on  the  reaction  TiCl4  +  2  H20f  see  [3.10]  ; 

injection  of  steam  and  liquid  nitrogen,  see  [3.14  -  3.16]; 
by  dropping  C02  pellets  in  hot  water;  generator  see  [3.17] 
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chapter  4)  .  The  smoke  wire,  an  electrically  controlled  device  (see  chapter  5)  ,  allows 
for  the  generation  of  very  fine  smoke  lines.  Direct  injection  of  the  smoke  is  performed 
by  means  of  a  smoke  pipe  (Fig.  3.2)  or  a  system  of  such  pipes  (’’rake")  . 


Fig.  3.2:  Smoke  is  injected  with  a  smoke  pipe  into  the  separated  flow  behind  an  exter¬ 
nal  mirror  at  the  front  door  of  a  car.  The  smoke  fills  the  entire  separated 
flow  regime.  The  experiment  has  been  performed  in  a  full-scale  windtunnel 
(Volkswagenwerk  AG,  Wolfsburg,  F.R.  Germany) . 

The  application  of  smoke  to  a  recirculating  wind  tunnel  creates  the  problem  that, 
after  a  certain  time  of  operation,  the  tunnel  is  completely  filled  with  smoke.  For  this 
reason,  and  also  because  any  smoke  or  aerosol  is  toxic  to  some  extent,  there  is  an  inte¬ 
rest  in  using  steam  or  water  fog  as  the  trace  material.  Steam  in  combination  with  a  coo¬ 
ling  agent  can  be  introduced  in  the  air  flow  where  a  visible  fog  is  produced  [3.14] . 

After  heat  exchange  with  the  surrounding  air,  the  fog  disappears,  leaving  the  main  air 
stream  clean  in  the  recirculating  system.  An  air  drying  device  might  be  needed.  Liquid 
nitrogen  may  serve  as  the  cooling  agent,  so  that  steam  and  liquid  nitrogen  are  expelled 
together  into  the  air  stream  through  a  mixing  nozzle  [3. 14 ,  3.16].  From  preliminary  ex¬ 
periments  performed  in  relatively  small  wind  tunnels,  it  has  been  reported  that  the  visi¬ 
bility  of  such  fog  compares  with  that  of  Ti02  mist  or  cigarette  smoke;  but  the  fog  "par¬ 
ticles"  are  larger  than  the  particles  of  the  two  former  materials. 

For  the  reasons  discussed  in  the  latter  paragraph,  many  smoke  flow  experiments  are 
performed  in  wind  tunnels  exhausting  directly  into  the  atmosphere.  Injection  of  the  smoke 
ahead  of  the  contraction  and  a  large  contraction  ratio  serve  to  stabilize  the  smoke  lines 
In  principle,  there  is  no  upper  limit  in  air  speed  for  the  application  of  smoke,  and 
smoke  lines  have  been  generated  even  in  supersonic  flow  [3.18],  Illumination  in  the  smoke 
experiments  is  mostly  conventional,  i.e.,  front  light  illumination  under  a  certain  angle 
with  respect  to  the  viewing  direction  with  conventional  sources,  e.g.,  mercury  lamps, 
halogen  lamps,  spot  lights.  More  recently,  illumination  in  form  of  a  plane,  thin  light 
sheet  is  applied  in  order  to  visualize  flow  structures  in  a  certain  cross  section  of  the 
flow  field  (see  chapter  5  )  . 

Beyond  the  use  in  wind  tunnels,  smoke  flow  visualization  has  been  applied  in  any 
kind  of  air  flows,  even  for  large  scale , outdoor  experiments  or  studies.  Its  special  role 
for  detecting  separated  flow  regimes  (Fig.  3.3),  vortices,  and  coherent  structures  in 
shear  layers  (Fig.  3.4)  should  be  emphasized. 
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Fig.  3.3:  Smoke  is  introduced  into  the  attached  boundary  layer  of  the  inclined  back  of 
a  car.  It  separates  with  the  boundary  layer  and  fills  the  entire  wake  (Volks- 
wagenwerk  AG,  Wolfsburg,  F.R.  Germany). 


Fig.  3.4:  Smoke  in  the  separated  flow  regime  behind  a  vertical  fence  visualizes  coherent 
structures  in  the  separated  shear  layer.  Reynolds  number  based  on  the  fence 
height  is  1.5  x  1o3  (H.H.  Fernholz,  Technische  Universitat  Berlin,  R.F .Germany) . 

3.2  VISUALIZATION  BY  DYE  IN  WATER 


The  visualization  of  water  flows  by  means  of  dye  is  as  popular  as  the  use  of  smoke 
in  air  flows.  The  various  dye  techniques  have  been  reviewed  by  Clayton  and  Massey  [3.19]  , 
and  by  Werle  [3.20],  The  particular  situation  of  dye  visualization  in  towing  tanks  has 
been  described  by  Gad-el-Hak  [3.21],  Unlike  smoke,  appropriate  dyes  can  fulfill  almost 
any  requirements  of  an  ideal  tracer. 


-  neutral  buoyancy, 

-  high  stability  against  mixing 

-  good  visibility. 
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Most  materials  are  non-toxic,  and  even  if  they  are  not  so,  the  toxicity  is  kept  with¬ 
in  the  liquid  flow  system  and  will  have  little  exchange  with  the  environment.  As  the  name 
says,  dye  visualization  includes  colour  as  an  additional  dimension  for  discrimination. 
Neutral  buoyancy  can  be  attained  by  mixing  the  dye  material  with  a  fluid  whose  specific 
weight  is  lower  than  that  of  water,  e.g.,  alcohol.  The  mixture  normally  is  not  a  true  so¬ 
lution,  and,  under  the  action  of  inertial  or  centrifugal  forces,  which  act  in  different 
ways  on  the  components  of  the  mixture,  the  dye  tracers  might  not  indicate  the  true  flow 
direction.  Such  a  discrepancy  is  avoided  by  complete  dilution  of  the  dye  in  the  water 
(food  coloring, ink) .  A  number  of  possible  dyes  are  listed  in  Table  3.2,  together  with  re¬ 
levant  references.  Milk,  which  can  be  mixed  with  certain  dyes,  particularly  food  coloring, 
preserves  the  dye  filaments  from  a  rapid  diffusion  into  the  water.  At  the  same  time,  milk 
improves  the  reflectivity  of  the  dye. 


Table  3 . 2 


Dyes  for  water  flow  visualization  and  respective  references 


food  coloring 

food  coloring  and  milk 

ink 

potassium  permanganate 
fluorescent  ink 
Rhodamine 
Fluoresceine 


[3.22  -  3.26] 
[3.21;  3.27?  3.28] 

[3.29] 

[3.30] 

[3.31] 

[3.32  -  3.34] 

[3.33  -  3.35] 


As  in  the  case  of  smoke  tracers,  dye  is  released  from  small  ejector  tubes  or  from  small 
orifices  in  the  wall  of  a  test  model.  Hypodermic  tubes  or  syringes,  or  small  Pitot  tu¬ 
bes  can  serve  as  dye  ejectors.  Conventional  front  light  illumination  is  sufficient  in 
most  cases. 

Special  lighting  devices  are  necessary  for  fluorescing  dyes.  These  dyes  must  be  excited 
by  a  light  source  which  radiates  in  the  spectral  range  of  greatest  fluorescence  of  the 
dye,  e.g.,  a  mercury  lamp  or  an  Argon  ion  laser.  The  exciting  illumination  should  be 
normal  or  oblique  to  the  direction  of  observation,  because  one  records  the  scattered 
fluorescent  radiation.  Table  3.2  includes  two  trademarks  of  fluorescing  dyes.  Rhodamine- 
6G ,  the  most  frequently  employed  material  of  this  type,  emits  a  yellow  radiation,  when 
excited  by  an  Ar+  laser,  Rhodamine-B  radiates  dark-red,  Fluoresceine  green. Rhodamine  is 
soluble  in  methanol,  which  might  be  of  interest  for  adjusting  the  specific  weight  of  the 
tracer  material.  Recent  developments  in  laser- induced  fluorescence  allow* for  determining 
the  concentration  of  the  dyed  fluid  from  the  measurement  of  radiation  intensity  [3.35]. 

The  flow  phenomena  visualized  by  means  of  dye  in  water  tunnels  are  the  same  as  those 
studied  by  smoke  visualization  in  air  (Figs. 3.5,  3.6).  The  simulation  of  aerodynamic  flow 
problems  in  water  tunnels  has  a  specific  experimental  advantage,  see,  e.g., [3.26;  3.36]; 
For  the  same  unit  Reynolds  number  and  model  scale,  the  velocity  in  water  is  1/15  of  that 
in  air,  so  that  these  phenomena  can  be  observed  at  a  relatively  low  speed. 

A  special  visualization  technique  has  been  developed  for  towing  tanks  [3.21;  3.32]. 
Thin,  horizontal  dye  layers  are  produced  in  the  water  before  the  model  is  set  into  motion. 
The  flow  around  the  model  displaces  the  colored  layers  whose  movement  can  be  observed 
and  recorded.  The  dye  layers  can  be  stabilized  by  a  weak  saline  stratification  in  the 
tank.  The  time  scale  for  diffusion  of  the  dye  layers  into  the  stratified  salt-water  is 
large  compared  to  the  time  needed  for  one  test  run,  so  that  the  system  of  horizontal  dye 
layers  can  be  used  for  several  experiments. 
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Fig.  3.5:  Dye  line  visualization  of  trailing  vortices  separating  from  a  delta  shaped 
wing  in  a  water  tunnel.  (H.Werle,  ONERA, 

France) 


Fig.  3.6:  Dye  visualization  of  boundary  layer  transition  and  separation  on  a  sphere 
in  a  water  tunnel  at  Re  =  3  .  1o5.  (H.  Werle,  ONERA,  France) 
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4.  ELECTRICAL,  CHEMICAL  AND  OPTICAL  CONTROL  OF  TRACER  PRODUCTION 


The  release  of  smoke  and  dye  from  injection  tubes  suffers  from  the  disturbance  of 
the  flow  by  these  mechanical  devices  and  a  lack  of  control  in  time  of  tracer  production. 
Another  shortcoming  is  that  these  devices  are  only  point  sources  of  the  trace  material. 
It  would  be  advantageous  if  the  material  could  be  produced  along  a  prescribed  curve  in 
the  flow  at  the  same  instant  of  time,  i.e.  released  from  a  line  source.  Orientation  of 
the  line  source  normal  to  the  main  flow  direction  would  result  in  the  visualization  of 
velocity  profiles.  It  will  be  shown  that  such  arrangements  can  be  realized  for  both  air 
and  water  flows.  The  line  source  is  either  defined  by  a  thin  wire,  and  the  tracer  pro¬ 
duction  is  controlled  by  an  electric  signal  along  the  wire;  or,  the  reaction  generating 
the  tracer  material  is  controlled  without  the  need  of  a  wire  by  a  thin  optical  or  elec¬ 
tromagnetic  beam  which  propagates  along  a  straight  line  through  the  fluid.  Table  4.1 
gives  a  survey  of  these  methods  which  will  be  discussed  in  the  following  subsections. 

Table  4.1 

Survey  of  methods  for  the  visualization  of  velocity  profiles 


fluid  method  tracer  tracer  generated  by  general  properties 


water 

tellurium 

method 

tellurium 

cloud 

wire  electrode 

non-reversible  tracer  production 
(contamination  of  fluid) ;  not 
neutrally  buoyant 

water 

hydrogen 

bubble 

method 

hydrogen 

bubbles 

wire  electrode 

no  contamination  of  fluid;  not 
neutrally  buoyant 

water 

thymol  blue 
method 

pH  indica¬ 
tor 

wire  electrode 

true  solution  (neutrally  buoyant) ; 
reversible  reaction  (no  contami¬ 
nation) 

organic 

fluid 

photochro- 

mism 

photochro- 
mic  dye 

UV  laser  beam 

neutrally  buoyant;  reversible  re¬ 
action  (no  contamination) ;  optical 
control  (no  wire  needed) 

air 

smoke  wire 

smoke 

pulse-heated 

wire 

contamination  by  smoke 

air 

spark  tra¬ 
cer  method 

line  of 
ionized  (lu¬ 
minous)  air 

electric  discharge 

no  contamination;  tracer  lighter 
than  working  fluid:  large  errors 
possible  in  accelarated  flow 

4 . 1  HYDROGEN  BUBBLE  METHOD 

Electrolysis  of  water  is  the  basic  principle  of  this  method.  It  is  known  that  hydro¬ 
gen  bubbles  are  formed  at  the  cathode  and  oxygen  bubbles  at  the  anode,  when  a  DC  voltage 
is  applied  between  two  electrodes  in  water.  The  hydrogen  bubbles  develop  with  much  smal¬ 
ler  size  than  the  oxygen  bubbles,  and  only  the  hydrogen  bubbles  are  used  as  flow  tracers. 

A  thin  wire  normal  to  the  mean  flow  direction  usually  serves  as  the  cathode.  A  short 
electric  pulse  in  the  electrolytic  circuit  generates  a  line  of  hydrogen  bubbles  along  the 
wire.  This  line  or  row  is  carried  away  with  the  flow,  and  it  deforms  according  to  the 
local  velocity  profile,  which  becomes  visible  when  the  bubble  line  is  seen  from  the  side 
(Fig.  4.1) .  By  pulsing  the  voltage  at  a  constant  frequency,  one  produces  several  succes¬ 
sive  bubble  lines.  The  correspondingly  marked  flow  curves  are  separated  by  a  constant 
time  interval,  and  they  are  called  "time  lines". 

The  hydrogen  bubble  flow  visualization  apparatus  is  commercially  available.  Plati¬ 
num  or  stainless  steel  wires  with  diameters  of  the  order  of  0.01  to  0.02  mm  are  used  as 
the  cathode.  The  anode  of  arbitrary  shape  is  placed  at  some  other  location  in  the  flow 
channel.  The  diameter  of  the  wire  determines  the  bubble  size.  A  rule  of  thumb  is  that  the 
bubble  diameter  is  of  the  order  of  the  cathode  wire's  diameter.  The  bubble  size  depends 
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a:  Generation  of  hydro¬ 
gen  bubble  time  lines 
in  the  boundary  layer 
on  a  flat  plate  (adop¬ 
ted  from  Smith  and 
Paxson  [4.1])  . 


platinum  wire 
(cathode) 


Ml iw  I'.nes 


pkste-  lanotse) 


Fig.  4.1  b:  Operation  of  the 

system  with  a  verti¬ 
cal  wire. 


also,  but  to  a  minor  extent,  on  the  conductivity  of  the  liquid  and  on  the  applied  voltage. 
Normal  tap  water  is  sufficient  for  the  electrolytic  process;  values  of  100  to  2oo  V  are 
typical  for  the  applied  voltage.  It  is  desirable  to  have  the  bubbles  as  small  as  possible. 
Then,  the  Reynolds  number  associated  with  the  rising  motion  (buoyancy)  of  the  bubbles  is 
very  low,  and  the  rise  velocity  determined  from  Stokes'  law  is  very  small,  thus  minimizing 
the  error  due  to  buoyancy. 

In  the  "time  line"  operation,  one  marks  a  line  of  fluid  elements  whose  position  co¬ 
incided  at  a  given  instant  with  the  position  of  the  wire.  A  system  of  time  lines  does  not 
provide  information  on  the  particle  trajectories.  These  trajectories  can  be  made  visible 
by  insulating  short  sections  of  the  wire.  The  lime  lines  are  then  interrupted  at  the  re¬ 
spective  positions,  and  connecting  these  interruptions  yields  the  direction  of  the  flow. 

If  one  uses  a  kinked  wire,  the  bubbles  will  only  develop  from  the  kinks;  and  with  a  conti¬ 
nuous  operation  of  the  electric  system  one  may  thus  visualize  streaklines  (or  streamlines) 
instead  of  velocity  profiles. 


From  several  sources  it  follows  that  an  angle  of  65°  between  the  illumination  system 
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and  the  direction  of"  observation  is  most  suitable.  This  might  have  to  do  with  the  scat¬ 
tering  characteristics  of  the  hydrogen  bubbles.  Recording  can  be  performed  with  conven¬ 
tional  photography  or  cinematography.  Video  recording  is  of  advantage  for  the  further  pro¬ 
cessing  of  the  flow  pictures,  particularly  if  the  flow  pictures  consist  of  such  well  de¬ 
fined  structures  like  the  time  lines,  which  can  easily  by  recognized  and  analyzed  in  re¬ 
spective  image  processing  systems  [4.1,  4.2].  The  period  of  time  during  which  the  bubb¬ 
les  persist  and  can  be  observed  in  the  flow  is  limited  by  the  dissolution  of  the  hydro¬ 
gen  bubbles  in  the  fluid.  Diffusion  of  the  bubbles  increases  with  Reynolds  number  and  is 
very  rapid  in  turbulent  flows.  The  application  of  the  method  is  therefore  restricted  to 
low-speed  flows,  maximum  velocities  being  of  the  order  of  20  to  30  cm/s. 

The  hydrogen  bubble  method  has  been  developed  in  the  early  sixties  [4.3  -  4.5]  from 
an  earlier  technique , Wortmann ' s  tellurium  method  [4.6],  which  uses  almost  the  same  appa¬ 
ratus,  but  with  a  cloud  of  tellurium  as  the  tracer  separating  in  the  electrolytic  reac¬ 
tion  from  a  tellurium  wire  electrode.  The  disadvantage  of  the  tellurium  is  that  it  does 
not  dissolve  in  the  water  and,  after  a  certain  time  of  operation,  will  contaminate  the 
water  in  a  recirculating  system.  Recent  progress  on  the  hydrogen  bubble  method  has  been 
summarized  by  Matsui  et  al.  [4.7].  The  method  has  become  a  standard  experimental  tool  for 
water  tunnels  and  towing  tanks,  and  its  applications  are  widespread.  It  has  enabled  ex¬ 
perimentalists  to  make  considerable  progress  in  studying  boundary  layer  instabilities 
and  in  recognizing  the  structure  of  turbulent  flows,  see,  e.g.  [4.8  -  4.13]  and  Fig.  4.2. 


(a)  y/R  =  1 


(b)  y/R  =0.75 


(c)  y/R  =  0.6 


(d)  y/R  =  0.5 


(e)  y/R  =  0.375 


(f)  y/R  =  0.25 


Fig.  4.2:  Time-line  visualization  of  the  vortex  system  developing  from  a  hemisphere 

(of  radius  R)  set  onto  a  flat  plate.  Photographs  are  top  views  with  the  hori¬ 
zontal  cathode  wire  at  different  heights  y  above  the  plate.  Flow  is  from  left 
to  right  (C.R.  Smith,  Lehigh  University,  Bethlehem,  PA,  USA.  ) 

Schraub  et  al.  [4.4]  presented  an  extensiv  discussion  of  the  error  sources  when  velocity 
fields  are  mapped  with  the  hydrogen  bubble  method.  The  most  significant  error  sources 
in  brief: 

-  uncertainties  in  measuring  the  distance  Ax  between  the  time  lines; 

-  averaging  the  velocity  over  the  distance  Ax? 


are 
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-  velocity  defect  of  the  bubbles  due  to  buoyancy,  particularly  in  accelerated  flows; 
velocity  defect  in  the  fluid  as  caused  by  the  wake  of  the  cathode  wire. 

4.2  DYE  PRODUCTION  BY  pH- INDICATORS 

An  electrolytic  reaction  may  change  the  pH-value  of  the  fluid  in  the  vicinity  of  the 
electrodes.  Due  to  the  migration  of  positive  hydrogen  ions  to  the  cathode  and  the  forma¬ 
tion  of  gas,  the  liquid  becomes  basic  near  the  cathode  and  consequently  acidic  near 
the  anode  (positive  electrode) .  A  number  of  pH-indicators  are  known  that  can  be  solved  in 
water  and  change  color  upon  a  specific  change  of  the  pH-value,  see  [4.14]  and  Table  4.2. 
For  the  methods  discussed  here,  indicators  are  of  interest  that  produce  a  visible  dye  for 
the  range  of  pH-values  reached  near  the  cathode.  With  a  thin  cathode  wire,  visible  time 
lines  (see  subsection  4.1)  can  be  produced  by  pulsing  the  voltage  between  the  electrodes. 
In  order  to  avoid  the  simultaneous  production  of  hydrogen  bubbles  at  the  cathode,  the 
applied  voltage  must  be  sufficently  low. 

Table  4.2 


pH-indicators  and  range  of  color  change  (data  taken  from  [4.14]  and  [4.17] 


indicator 


pH ( 1 ) /color 


pH ( 2) /color 


bromo  cresol  green 
bromo  phenol  red 
bromo  thymol  blue 
meta  cresol  purple 
thymol  blue 
phenolphtalein 


4 .0/yellow 
5 . 2/yellow 
6 . 0/yellow 
7 .6/yellow 
8 .0/yellow 
8 . 3/clear 


5 . 4 /blue 
7.0 /red 
7. 6 /blue 
9. 2 /purple 
9 . 6 /blue 
1o  /red 


The  indicator  most  commonly  used  in  this  application  is  thymol  blue.  Its  use  for 
flow  visualization  and  the  measurement  of  velocity  profiles  have  been  described  by  Ba¬ 
ker  [4.15].  The  working  fluid  is  an  aqueous  solution  of  thymol  blue  which  is  orange- 
yellow  in  an  acidic  environment  (pH  <  8.0)  and  turns  its  color  to  blue  if  the  solution 
becomes  basic  (pH  >  9.6).  The  working  fluid  is  prepared  by  solving  0.01  to  0.04  %  by 
weight  of  thymol  blue  in  distilled  water.  A  few  drops  of  an  acid,  e.g.  HC1,  are  added 
so  that  the  solution  is  close  to  the  point  of  becoming  basic.  The  electrode  system  is 
similar  to  that  used  for  the  hydrogen  bubble  method.  Applied  voltages  are  of  the  order 
of  1 0  V  at  currents  of  about  1o  mA  or  less.  The  optical  contrast  between  the  orange- 
yellow  liquid  and  the  blue  dye  produced  at  the  cathode  can  be  enhanced  by  employing  yel¬ 
low-light  illumination  from  a  sodium  lamp. 

The  aqueous  solution  of  thymol  blue  is  a  true  solution,  and  the  marked  fluid  ele¬ 
ments  are  neutrally  buoyant.  The  technique  is  therefore  used  with  preference  for  flows 
that  are  dominated  or  driven  by  mass  forces,  e.g.  thermal  convection,  stratified  flows, 
and  flows  in  rotating  systems.  Since  it  is  costly  to  fill  a  water  tunnel  or  towing  tank 
with  distilled  water,  and  because  it  is  difficult  to  control  the  pH-value  in  such  large 
reservoirs,  the  applications  are  restricted  to  relatively  small  volumes  of  the  working 
fluid.  Another  limitation  is  the  maximum  velocity  or  Reynolds  number,  because  the  dif¬ 
fusion  of  the  dye  into  the  surrounding  fluid  increases  with  increasing  Reynolds  number. 
The  reported  maximum  velocity  for  the  production  and  identification  of  time  lines  is  a- 
bout  4  cm/ s . 
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Variations  of  the  pH-value  take  place  in  the  mixing  region  of  two  fluid  flows  that 
were  given  different  pH-values  before  they  mix.  Color  changes  of  a  pH-indicator  solved 
in  the  fluid  may  occur  in  the  mixing  region.  This  technique  enabled  Breidenthal  [4.16] 
to  visualize  the  structure  of  the  mixing  layer  between  two  parallel  fluid  flows  of  dif¬ 
ferent  value  of  the  initial  velocity  (Fig.  4.3).  Extinction  measurements  allow  for  de¬ 
termining  the  fluid  concentration  in  the  mixing  layer.  Gardner  [4.17]  showed  that  the 
use  of  color  is  another  means  for  such  concentration  measurements  in  mixing  problems.  He 
uses  two  pH-indicators  (thymol  blue,  bromocresol  green)  simultaneously,  so  that  several 
color  changes  can  be  observed  over  a  wide  range  of  the  pH-value. 


Development  of 
instabilities  in 
the  mixing  of  two 
fluid  flows  of 
different  mean  ve¬ 
locities,  visuali¬ 
zed  by  a  pH-indi- 
cator.  Bottom:  si¬ 
de  view;  top:  plan 
view  (R.  Breiden¬ 
thal,  University 
of  Washington,  Se¬ 
attle,  WA,  USA). 


4 . 3  PHOTOCHROMIC  DYE  PRODUCTION 


Popovich  and  Hummel  [4.18]  reported  on  the  use  of  photochromic  materials  in  liquid 
solution  for  the  production  of  dye  by  irradiation  of  appropriate  light.  The  dye  can  be 
produced  along  a  straigt  line,  e.g.  by  using  a  respective  laser  source,  and  the  dye  pro¬ 
ducing  light  beam  can  be  controlled  externally  in  time  and  position,  so  that  no  distur¬ 
bing  mechanical  element  has  to  be  introduced  in  the  flow.  Photochromism  is  defined  as  the 
reversible  transition  of  a  chemical  substance  between  two  states,  A  and  B,  exhibiting  no¬ 
ticeably  different  absorption  spectra.  The  usual  situation  is  that  state  A  has  high  ab¬ 
sorption  rates  in  the  ultraviolet  (UV) ,  while  B  is  highly  absorbing  in  the  visible  range 
of  wavelengths.  Therefore,  the  substance  is  invisible  (transparent)  in  A  and  visible 
(opaque)  in  B.  The  transition  A  B  is  stimulated  by  UV  radiation,  and  it  is  desired  that 
the  reverse  transition  B  A  is  also  possible;  for  the  substances  which  are  of  interest 
here,  the  reverse  transition  is  spontaneous  and  can  be  enhanced  by  thermal  energy. 


Nb 

Na 


Fig.  4.4:  Principle  of  the 
transition  of  the 
photochromic  mate¬ 
rial  between  two 
states  of  different 
absorption,  A  and  B 
Na,  Nb  =  Number  of 
molecules  in  states 
A  and  B,  respective 
ly;  t  =  time 


The  principle  of  the  photochromic  transition  process  is  shown  in  Fig.  4.4,  where  the 
ratio  of  the  numbers  of  tb  nolecules  in  the  two  states  A  and  B,  NB/NA, is  plotted  as  a 
function  of  time.  The  value  of  this  ratio  increases  rapidly  during  the  short  period  of 
incidence  of  the  stimulating  UV  radiation,  and  it  decreases  slowly  during  the  reverse 
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transition  when  the  fluid  regains  its  initial  transparency.  Substances  exhibiting  the 
described  behaviour  are  components  known  as  "pyridines" /  e.g.  [4.18  -  4.20],  and  Mspy- 
ranes" ,  e.g.  [4.21  -  4.24].  These  substances  are  unsoluble  in  water.  Solutions  can  be 
prepared  with  a  number  of  organic  liquids,  e.g.  ethanol,  toluol,  kerosene.  Concentration 
of  the  photochromic  substance  should  be  between  0.1  and  0.2  %  by  weight  and  must  be  cho¬ 
sen  accordingly  to  the  special  experimental  conditions:  Too  low  a  concentration  results 
in  a  weak  dye  production?  too  high  a  concentration  can  cause  an  unwanted  temperature  in¬ 
crease  in  the  fluid,  because  only  a  fraction  of  the  incident  energy  is  used  for  the  sti¬ 
mulation  of  state  B,  while  the  major  portion  of  this  energy  is  converted  into  heat. 

Most  suitable  as  a  source  for  the  stimulating  UV  radiation  are  a  ruby  laser  equiped 
with  a  frequency  doubler,  or  a  N2  laser  which  can  be  built  at  relatively  low  cost  [4.25] . 
The  photochromic  dye,  i.e.,  the  opaque  state  B,  is  produced  along  the  laser  beam.  When 
the  laser  is  pulsed  at  equally  separated  time  intervals,  a  system  of  time  lines  can  be 
made  visible  (Fig.  4.5).  The  flow  channel  must  be  equipped  with  quartz  windows,  because 
normal  glass  is  not  transparent  for  the  UV  radiation. 


Fig.  4.5:  Photochromic  visualization  of  ve¬ 
locity  profiles  in  the  develop¬ 
ment  of  a  drop  (E.  Marschall, 
University  of  California,  Santa 
Barbara,  CA,  USA) . 


Experiments  performed  with  photochromic  dye  are  restricted  by  the  particular  proper¬ 
ties  of  the  solving  liquid.  One  consequence  is  the  limitation  of  the  volume  of  the  flow 
facility.  Reported  applications  concern  the  study  of  the  near  wall  region  in  boundary  la¬ 
yers  [4.18,  4.22],  or  the  measurement  of  velocity  in  the  interior  of  growing  drops  [4.24]. 
D'Arco  et  al.  [4.23]  produce  in  the  interior  of  the  fluid  a  rectangular  grid  by  means  of 
an  appropriate  illumination  pattern,  and  they  observe  the  deformation  of  this  grid  due  to 
the  flow;  this  allows  for  determining  shear  stresses  in  the  observed  plane. 

4.4.  *  SMOKE  WIRE 


The  hydrogen  bubble  method  has  its  equivalent  counterpart  for  air  flows  in  form  of 
the  smoke  wire,  which  has  been  proposed  first  by  Yamada  [4.26],  and  then  refined  and  de¬ 
veloped  further  by  several  investigators  [4.27  -  4.30].  Not  really  smoke  but  oil  fog  is 
evaporated  from  an  electrically  heated  wire  that  serves  as  the  line  source  of  the  tracer 
material,  like  the  electrode  wires  as  described  in  the  subsections  4.1  and  4.2.  A  stain¬ 
less  steel  wire,  approximately  0.1  mm  in  diameter,  is  oriented  vertically  or  horizontal¬ 
ly  and  normal  to  the  main  flow.  It  can  be  mounted  on  a  probe  or  span  the  whole  cross- 
section  of  a  wind  tunnel.  The  wire  is  coated  with  a  very  thin  layer  of  kerosene.  The  coa¬ 
ting  is  done  manually  by  wiping  the  wire  with  a  fine,  small  brush  or  a  cotton-tipped  app- 


licator.  Corke  et  at.  [4.27]  report  on  a  device  with  which  a  drop  of  oil  is  forced  down 
the  wire  by  applying  pressure  to  a  small  oil  reservoir. 
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black  wall  strobe  tight 


flow 


Fig.  4.6:  Schematic  diagram  of  a  smoke  wire  system  with  synchronized  strobe  light 

illumination  and  camera.  The  system  is  operated  by  pushing  the  camera  button 
(after  [4.29] ) . 

With  a  uniform,  fine  coating  and  with  the  application  of  an  equally  spaced,  pulsed 
electric  heating  it  is  possible  to  generate  time  lines  equivalent  to  those  obtained  with 
the  hydrogen  bubble  or  thymol  blue  method.  This  is  indicated  in  Fig.  4.6,  where  a  schema¬ 
tic  block  diagram  of  such  a  system  with  synchronized  strobe  lieght  illumination  and  came¬ 
ra  is  shown.  The  applications  of  the  smoke  wire  for  generating  time  lines,  however,  are 
rare  due  to  the  rapid  diffusion  of  the  smoke  into  an  air  stream.  More  frequent  is  the  use 
of  the  wire  for  producing  very  fine  smoke  lines  with  exact  spatial  and  temporal  control. 
Then,  the  smoke  line  is  either  observed  from  a  direction  parallel  to  the  wire;  or  the 
coating  is  done  a  little  thicker,  so  that  small  oil  beads  or  minute  droplets  form  along 
the  wire,  and  at  each  of  the  beads  a  smoke  filament  originates  when  the  wire  is  heated. 
The  spacing  between  such  lines  is  influenced  by  the  wire  size,  the  amount  of  oil  in  the 
coating,  and  the  oil  surface  tension. 

The  smoke  wire  is  a  device  for  introducing  smoke  into  an  air  stream  with  much  grea¬ 
ter  precision  than  with  an  injection  tube,  and  it  has  already  found  a  wide  application; 
see,  e.g.  [4.31],  The  smoke  wire  also  allowed  for  the  (so  far  only  reported)  generation 
of  colored  smoke  [4.32], 

4.5.  SPARK  TRACER  TECHNIQUE 


The  application  of  the  smoke  wire  is  restricted  to  low  air  velocities.  At  higher 
velocities,  a  mapping  of  the  velocity  field  is  possible  by  making  use  of  an  electric  dis¬ 
charge  in  the  air  flow.  The  principle  of  the  method  is  explained  with  the  aid  of  Fig.  4.7. 
An  electric  spark  discharge  is  generated  between  the  two  electrodes  upon  the  application 
of  a  high  voltage  (several  kV) .  The  spark  creates  in  the  gas  an  ionized,  luminous  air  co¬ 
lumn  (plasma) , whose  luminosity  persists  for  a  period  of  time  of  up  to  loo  p,sec ,  depending 
on  the  discharge  conditions.  The  plasma  column  can  be  observed  while  it  is  swept  down¬ 
stream  by  the  main  flow.  A  series  of  time  lines  can  be  produced  by  repeating  the  spark 
discharge:  Since  the  luminous  column  has  a  smaller  electric  resistance  than  the  surroun¬ 
ding  neutral  gas,  a  second  spark  produced  during  the  lifetime  of  the  plasma  column  will 
rather  follow  this  preionized  path  than  take  the  shortest  (straight)  way  between  the  elec¬ 
trodes.  In  a  series  of  sparks,  each  spark  will  follow  and  thereby  re-illuminate  the  column, 
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Fig.  4.7:  Principle  of  spark-tracer  technique  with  spiked  electrodes,  applied  to  three 

different  flows:  uniform,  parallel  flow  (a) ;  wake  flow  behind  sphere  (b) ; 
flow  over  an  airfoil  profile  (c) . 

which  was  generated  by  the  first  spark  and  swept  downstream  by  the  air  flow.  Photogra¬ 
phing  the  series  of  discharge  with  an  open-shutter  camera  results  in  the  visualization 
of  a  system  of  time  lines.  From  the  known  frequency  of  the  spark  series  and  the  measured 
displacement  one  may  derive  the  local  velocity  field. 

The  form  of  the  electrodes  must  be  such  that  it  fulfills  two  requirements:  It  should 
cause  a  minimum  of  interference  with  the  flow,  and  it  should  provide  the  first  spark  a 
straight  path  between  the  electrodes,  normal  to  the  main  flow.  The  result  is  a  compromise 
in  form  of  electrode  rods,  and  the  discharge  of  each  spark  starts  from  a  different  point 
along  the  rod,  according  to  the  instantaneous  position  of  the  preionized  gas  column.  The 
possible  distance  between  the  rods  is  several  cm,  with  the  maximum  distance  depending  on 
the  air  pressure  and  the  value  of  the  applied  voltage.  The  appropriate  spark  frequency  de¬ 
pends  on  the  mean  flow  velocity.  It  can  be  up  to  loo  kHz.  The  problem  then  is  to  recharge 
the  high-voltage  capacitors  within  the  intervals  of  a  few  micro-seconds  each.  The  value 
of  the  lower  limit  of  the  spark  frequency  is  determined  by  the  diffusion  of  the  ionized 
gas  into  the  ambient  gas. 


Fig.  4.8:  Spark-tracer  technique  with  rod  electrodes,  applied  to  an  air  flow  with 
strong  disturbances  ( Impulsphysik  GmbH,  Hamburg,  F.R.  Germany). 
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The  spark  tracer  technique  has  been  developed  and  adapted  to  the  study  of  boundary 
layer  flow  by  Bomelburg  et  al.  [4.33].  Improvements  of  the  system  are  due  mainly  to  the 
work  of  Friingel  and  his  co-workers;  see,  e.g.  [4.34,  4.35].  It  has  been  used  for  measu¬ 
ring  gas  velocities  in  high-Mach  number  flow  facilities,*  see,  e.g.  [4.38,  4.39],  in  mo¬ 
del  studies  of  turbo-machines  [4.40,  4.41]  and  other  flow  machinery  with  swirling  flow 
[4.42];  see  also  Fig.  4.8. 


The  application  of  the  spark  tracer  method  may  suffer  from  a  number  of  error  sources, 
particularly  in  accelerated  flow  fields.  These  errors  have  been  discussed  extensively  by 
Rudinger  and  Somers  [4.43]  and  by  Matsuo  et  al.  [4.44].  The  plasma  column  as  the  tracer 
is  lighter  than  the  working  fluid  (air)  and  therefore  has  a  different  inertia.  A  rela¬ 
tive  motion  of  the  ionized  column  with  respect  to  the  mean  flow  transforms  the  tracer  in¬ 
to  a  vortex.  This  vortex  consumes  energy  from  the  main  flow  and  it  may  induce  an  additio¬ 
nal  velocity  component  to  the  tracer.  The  velocity  lag  between  tracer  and  flow  increases 
with  the  magnitude  of  the  acceleration  and  with  the  density  difference  between  tracer  and 
fluid.  In  both  references  it  is  found  that  the  relative  errors,  in  extreme  situations, 
can  exceed  5o  %. 
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5.  LIGHT  SHEET  METHODS 


In  the  previous  chapters,  no  particular  requirements  were  made  with  respect  to  the 
illumination  system.  In  this  chapter,  methods  are  discussed  for  which  the  illumination  is 
provided  in  form  of  a  plane,  thin  sheet,  so  that  only  the  tracer  material  within  thin 
sheet  or  slice  is  illuminated.  Illumination  in  form  of  such  a  sheet  can  be  applied  to  the 
methods  already  discussed,  but  it  also  allowed  for  the  development  of  new  techniques,  e.g. 
speckle  velocimetry  and  laser-induced  fluorescence,  both  of  which  are  included  in  this 
chapter. 

5.1.  GENERATION  OF  A  LIGHT  SHEET 


A  light  sheet  can  be  produced  with  a  conventional  lamp  (mercury  lamp,  halogen  lamp, 
spot  light) .  A  much  better  way  is  to  expand  the  beam  from  a  powerful  laser  by  means  of  a 
cylindric  lens  in  one  plane.  The  thickness  of  such  a  sheet  can  be  made  much  smaller  than 
1  mm.  Observation  of  the  flow  pattern  in  the  illuminated  plane  is  normal  to  this  sheet  or 
under  a  certain  angle  to  it  (Fig.  5.1) .  Illumination  by  a  light  sheet  can  be  applied  to 
both  dye  and  smoke  visualization,  and  it  is  most  appropriate  for  recognizing  vortical 
structures  in  these  flows  (Fig.  5.2).  The  light  sheet  technique  and  its  use  in  convential 


Fig.  5.1:  Generation  of  a  light  sheet  in  a  flow  by  expanding  a  laser  beam  in  one  plane 
by  means  of  cylindrical  lenses. 

visualization  systems  have  been  reviewed,  e.g.,  by  Veret  [5.1],  and  the  reported  applica¬ 
tions  are  numerous.  Koga  et  al.  [5.2]  describe  how  the  laser  light  for  the  sheet  can  be 
directed  to  the  appropriate  position  by  a  fiber  optics. 

Flow  facilities  have  been  designed  particularly  for  the  application  of  a  light  sheet 
normal  to  the  main  flow  direction.  These  special  wind  tunnels  are,  to  a  certain  degree, 
filled  with  either  smoke  or  vapor,  and  the  technique  is  called  a  smoke  screen  or  vapor 
screen,  respectively.  They  allow,  e.g.,  for  the  visualization  of  trailing  vortices,  which 
separate  from  test  models.  A  schematic  view  of  the  test  section,  light  source,  vapor  (or 
smoke)  screen,  and  recording  camera  is  shown  in  Fig.  5.3. 

While  the  application  of  smoke  tunnels  is  restricted  to  subsonic  flow  [3.4],  vapor 
screens  are  very  often  produced  to  observe  supersonic  flows.  The  supersonic  tunnel  must 
be  run  with  moist  air; the  air  cools  as  it  expands  through  the  nozzle,  and  the  moisture 
condenses  to  form  a  fog  in  the  test  section.  The  uniform  distribution  of  fog  in  a  cross 
section  normal  to  the  tunnel  axis  is  disturbed  due  to  the  wake  from  a  test  model.  McGre- 


32 


gor  [5.5]  has  investigated  the  rate  of  humidity  required  to  produce  a  satisfactory  screen 
at  Mach  numbers  up  to  M  =  2,  and  as  a  function  of  tunnel  pressure  and  temperature.  Snow 
and  Morris  [5.6]  systematically  investigated  the  operational  conditions  of  the  vapor 
screen.  They  also  describe  a  periscope  system  for  observation  so  that  it  is  not  necessa¬ 
ry  to  have  the  camera  inside  the  tunnel. 


Fig.  5.2:  Smoke  visualization  of  leading-edge  vortices  on  a  swept  delta  wing;  illumina¬ 
tion  by  light  sheet  (R.C.  Nelson,  University  of  Notre  Dame,  IN,  USA  [5.4]) . 


Fig.  5.3:  Illumination  by  light  sheet  in  a  vapor  tunnel  or  smoke  tunnel  (adopted  from 
Snow  and  Morris  [5.6]). 

The  light  sheet  plays  an  important  role  in  measuring  the  velocity  of  tracer  partic¬ 
les  in  a  flow  from  a  photographic  record.  In  this  case  it  is  mandatory  to  use  a  pulsed 
light  source  (flash  lamp,  electric  spark,  ruby  laser)  with  a  definite  pulse  length,  and, 
if  needed,  with  the  possibility  of  generating  a  series  of  such  pulses  (stroboscope) .  With 
the  light  sheet,  one  visualizes  the  two  velocity  components  in  the  plane  of  the  sheet. 
Emrich  [5.7]  calls  the  record  of  such  a  plane  a  "time  window",  within  which  the  velocity 
can  be  determined  simultaneously  for  many  points  (this  is  also  called  a  whole-field  re¬ 
cord)  .  Shifting  the  sheet  in  a  direction  normal  to  its  plane  enables  one  to  scan  a  three- 
dimensional  flow,  provided  that  the  flow  is  stationary  over  the  test  period.  The  third 
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velocity  component,  normal  to  that  plane,  is  not  recovered  in  this  way.  This  would  need  to 
illuminate  light  sheets  perpendicular  to  the  first  set  of  sheets. 

The  tracer  particles  in  the  light  sheet  are  observed  or  photographed  in  a  direction 
approximately  normal  to  the  illuminated  plane.  If  one  takes  a  time  exposure,  each  partic¬ 
le  appears  in  the  image  plane  as  a  streak  whose  length  is  proportional  to  the  average  lo¬ 
cal  velocity  (Fig.  5.4).  This  simple  way  of  measuring  a  velocity  vector  is  associated  to 
a  number  of  technical  shortcomings  and  possible  errors,  e.g.,  precision  and  control  of 
the  exposure  time;  optical  parallax  and  distortion  in  viewing  the  illuminated  plane;  in¬ 
ability  to  determine  the  flow  direction.  Stroboscopic  illumination  results  in  a  series 
of  particle  images,  usually  three  or  four  [5.7,  5.8],  indicating  the  particle  velocity  by 
magnitude  and  direction. 

A  new  means  for  improved  data  evaluation  is  image  processing  of  the  particle-streak 
record.  Dimotakis  et  al .  [5.1o]  have  scanned  and  digitized  the  particle-streak  pattern 

from  the  photographic  negative  so  that  this  pattern  can  be  reproduced  as  a  digital  image. 
It  is  then  possible  to  apply  special  computer  programs  for  displaying  the  velocity  distri¬ 
bution,  or  for  manipulating  the  data,  e.g.  superimposing  a  constant  velocity,  deriving 
stream  function,  vorticity,  or  other  flow  parameters,  see,  e.g.  [5.11  -  5.13]. 


Fig.  5.4:  Time  exposure  of  the  vortex  formation  in  the  water  flow  around  a  rotating 
cylinder.  The  arrow  indicates  the  direction  of  rotation.  Tracer  particles: 
lycopodium?  Re  =  1o3  (M.  Coutanceau,  University  of  Poitiers,  France) . 


5 . 2  SPECKLE  VELOCIMETRY 


The  techniques  of  speckle  photography,  which  are  used  in  solid  state  mechanics  for 
measuring  surface  deformations,  have  been  adopted  to  the  measurement  of  flow  velocity  in 
fluids  seeded  with  tracer  particles.  This  particular  application  of  speckle  photography 
requires,  in  most  cases,  the  use  of  a  light  sheet.  The  technique  is  a  special  recording 
method,  rather  than  a  direct  visualization  method.  It  allows  for  both  quantitative  mea¬ 
surement  of  the  two  velocity  components  in  the  illuminated  plane,  and  the  generation  of 
a  visual  pattern  in  form  of  iso-velocity  contours.  This  method  has  been  developed  in  a 
number  of  laboratories,  and  it  is  interesting  to  note  that  the  first  results  from  the 
different  groups  have  been  reported  apparently  all  in  the  same  year  [5.14  -  5.17]. 
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Fig.  5.5:  Displacement  of  an  object  point  by  amount  ds  and  recording  by  a  photographic 
double  exposure. 

The  principle  of  the  method  is  explained  by  considering  an  object  point  which  is  im¬ 
aged  onto  a  recording  plane  (Fig.  5.5).  The  object  point  is  displaced  by  a  small  distance 
ds  and,  under  the  assumption  of  an  imaging  ratio  of  1:1,  the  image  point  is  displaced  in 
the  same  direction  and  by  the  same  amount  ds.  In  a  photographic  double  exposure,  one  ex¬ 
posure  taken  before  and  the  second  after  the  displacement,  the  image  point  appears  twice, 
separated  by  ds  as  shown  in  Fig.  5.5.  In  the  application  to  a  fluid  flow,  the  object 
point  represents  a  single  tracer  particle  moving  along  the  flow  with  the  velocity  w.  The 
displacement  of  the  particle  between  the  two  photographic  exposures  is  ds  =  w  •  dt,  where 
dt  is  the  time  interval  separating the  two  exposures.  With  known  time  interval  dt,  the  de¬ 
termination  of  the  velocity  w  requires  the  measurement  of  the  distance  ds.  Speckle  photo¬ 
graphy  is  a  means  for  measuring  both  direction  and  magnitude  of  small  in-plane  displace¬ 
ments  for  a  whole  field  of  view.  The  principle  of  such  measurement  will  be  explained  next- 


Fig.  5.6:  Illumination  of  the  photographic  double  exposure  of  Fig.  5.5  ( " specklegram" ) 
by  a  thin  laser  beam  and  formation  of  Young's  fringes  in  a  plane  of  observa¬ 
tion  . 


After  development  of  the  double  exposure  the  two  image  points  appear  on  the  plate  as 
two  bright  spots  on  a  dark  background.  These  two  points  are  illuminated  from  behind  with 
a  thin  laser  beam  (Fig.  5.6) .  For  the  propagation  of  the  light  in  the  space  to  the  right 
of  the  plate,  the  two  illuminated  spots  can  be  considered  as  two  point  sources;  they  emit 
light  into  two  cones  ("halo") .  In  the  overlapping  regime  of  the  two  cones  the  light  from 
the  two  sources  can  interfere.  Parallel  and  equally  spaced  interference  fringes  ("Young's 
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fringes")  can  be  observed  in  a  plane  normal  to  the  optical  axis  and  at  a  distance  Z  from 
the  two  sources.  The  fringes  are  normal  to  the  direction  of  ds  and  separated  by  the  spa¬ 
cing  A: 


A  is  the  laser  wave  length.  The  measurement  of  the  vector  ds  is  thereby  reduced  to  the 
measurement  of  the  fringe  direction  and  fringe  spacing  in  a  system  of  Young's  interfe¬ 
rence  fringes. 


The  laser  beam  in  the  set-up  of  Fig.  5.6  illuminates  on  the  developed  plate  a  small 
circular  area  of  diameter  j<  1  mm.  If,  as  it  has  been  assumed  here,  only  one  pair  of  par¬ 
ticle  images  is  present  in  that  circular  area,  the  flow  has  been  seeded  obviously  at 
a  low  seeding  rate  so  that  it  is  possible  to  identify  or  to  image  individual  tracer  par¬ 
ticles.  This  is  not  a  common  situation.  At  higher  seeding  rates  the  circular  area  may  in¬ 
clude  a  number  of  pairs  of  particle  images.  The  recording  by  the  double  exposure  might 
not  even  be  an  imaging  process,  but  a  process  of  multiple-interference  caused  by  the  scat¬ 
tering  of  light  from  the  tracer  particles.  Then,  the  result  of  the  recording  is  a  real 
speckle  pattern,  a  granular  structure  consisting  of  a  distribution  of  "speckles".  In  the 
following,  the  term  "speckle"  is  used,  and  it  is  understood  that  the  speckle  is  a  parti¬ 
cle  image  when  the  seeding  rate  was  low. 


Hence,  the  small  circular  area  illuminated  by  the  laser  beam  includes  a  number  of 
pairs  of  speckles  which  are  distributed  randomly  in  that  area  (Fig.  5.7) .  If  all  these 
speckles  were  displaced  by  the  same  amount  and  in  the  same  direction,  the  resulting  sys- 


specklegram 


Fig.  5.7:  Random  distribution  of  speckle 
pairs  in  the  specklegram.  The 
circular  area  of  diameter  <  1  mm 
is  illuminated  by  the  scanning 
laser  beam. 


Fig.  5.8:  Young's  fringe  pattern  with  high  (a)  and  low  (b)  degree  of  correlation 
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tem  of  Young's  fringes  will  be  the  same  as  if  only  one  pair  had  been  illuminated.  But  if 
the  illuminated  speckles  have  not  all  experienced  the  same  displacement,  different  sys¬ 
tems  of  Young's  frinqes  will  be  generated , and  they  overlap  each  other.  The  pattern  is  to 
a  certain  degree  decorrelated, the  fringe  system  is  of  reduced  contrast  or  visibility 
(Fig.  5.8) .  If  it  is  still  possible  to  identify  a  dominating  frequency  in  the  fringe  pat¬ 
tern,  i.e.,  to  measure  a  mean  value  of  the  fringe  spacing  A,  then  one  attributes  a  mean 
value  of  the  displacement  ds  to  all  speckles  appearing  in  the  circular  area.  A  reduction 
in  size  of  this  area,  e.g.,  by  focusing  the  laser  beam  with  a  lens,  is  desirable  for  im¬ 
proving  the  local  resolution. 

Single  speckles  illuminated  near  the  edge  of  the  circular  area  (Fig.  5.7)  contribute 
to  the  noise  level  of  the  signal,  i.e.,  they  reduce  the  fringe  visibility. 

The  discussion  of  a  single  pair  of  speckles  in  the  field  can  serve  to  explain  some 
limitations  of  the  method.  Each  speckle  has  a  certain  diameter  d^  whose  value  can  be  de¬ 
rived  from  diffraction  theory  (see  e.g.  Frangon  [5.18]).  d^  is  equivalent  to  the  particle 
image  diameter  when  the  method  is  explained  in  terms  of  an  imaging  process.  The  lower  li¬ 
mit  for  the  displacement  ds  that  can  be  measured  in  this  way  is  given  by  the  requirement 
ds  >  dp.  There  exists  also  an  upper  limit  when  the  fringe  spacing  A  becomes  too  small,  i. 
e.,  when  the  fringe  system  becomes  too  narrow.  A  rule  of  thumb  is  that  this  upper  limit 
is  reached  for  ds  ~  2o  •  d^.  From  these  two  values  follows  the  range  of  tracer  velocities 
that  can  be  measured  in  an  experiment  with  given  time  interval  of  the  double-exposure. 

Most  of  the  flow  experiments  are  performed  with  a  light  sheet  as  shown  in  Fig.  5.1. 
Very  suitable  for  taking  the  double  exposure  is  the  use  of  a  double-pulsed  ruby  laser  beam 
expanded  in  one  plane  by  a  cylindrical  lens.  The  width  of  the  single  pulses  as  well  as 
the  time  interval  between  pulses  can  be  controlled  short  enough  to  allow  for  a  tolerable 
value  of  the  displacement  to  be  measured.  The  concentration  of  tracer  particles  In  the 
flow  determines  whether  the  method  is  working  in  the  particle-image-mode  (low  concentra¬ 
tion)  or  as  real  speckle  photography  (high  concentration);  see,  e.g.  [5.19].  Particle 
concentration  and  particle  size  have  an  essential  influence  on  the  amount  of  light  avail¬ 
able  for  taking  the  double  exposure,  and  this  problem  is  probably  the  most  servere  re¬ 
striction  for  the  application  of  speckle  velocimetry.  The  light  intensity  distribution  is 
quite  different  from  laser  Doppler  velocimetry,  where  all  the  light  is  concentrated  in 
the  focal  point,  i.e.  the  measuring  volume.  Here,  however,  the  laser  light  is  expanded  in 
a  sheet,  and  the  incident  intensity  which  is  then  scattered  from  an  individual  particle 
is  by  orders  of  magnitude  lower  than  in  laser  Doppler  systems.  The  low  sensitivity  of 
photographic  material  to  radiation  in  the  red  (particularly  the  ruby-laser  wavelength!) 
contributes  to  this  intensity  problem.  Consequently,  the  reported  applications  of  speckle 
velocimetry  were  either  cases  with  relatively  large  tracer  particles,  or  flows  of  low  ve¬ 
locity,  so  that  the  duration  of  each  light  pulse  was  relatively  long. 

Besides  by  using  powerful  laser  sources,  the  intensity  problem  can  be  improved  by 
taking  multiple-exposures  instead  of  a  double-exposure  [5.2o].This  is  also  of  advantage 
for  the  clarity  of  the  Young's  fringes,  but  the  method  is  then  restricted  to  flow  fields 
with  smooth  and  moderate  velocity  changes.  Another  possibility  is  to  operate  the  method 
in  forward  scattering  [5.21],  but  this  is  a  possible  solution  for  plane  flow  fields  only. 

Illumination  for  the  recording  can  be  provided  also  with  a  white  source  [5.22,5.23]? 
since  the  use  of  such  a  source  merely  excludes  the  possibility  of  interference,  the  re¬ 
cording  is  then  strictly  a  particle-image  process. 

The  analysis  of  the  developed  double-exposure  (which  is  called  a  " specklegram")  with 
the  method  described  in  Fig.  5.6  is  referred  to  as  the  "point-by-point "  reconstruction  of 
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Fig.  5.9:  System  of  fully  automated  evaluation  of  a  specklegram  (Erbeck  [5.25]). 

the  plane  velocity  field.  The  specklegram  is  scanned  point-by-point  with  the  laser.  For 
each  illuminated  point  one  obtains  a  pattern  of  Young's  fringes.  From  the  measured  fringe 
spacing  and  fringe  direction  one  determines  the  two  velocity  components  in  the  respective 
point  of  the  illuminated  plane.  The  information  stored  in  the  specklegram  is  a  very  dense 
distribution  of  data  on  displacement  vectors  or  velocity.  The  rational  analysis  of  such  an 
extended  data  field  requires  the  use  of  an  automated  processing  system.  Meynart  [5.24] 
describes  an  interactive  system  that  needs  the  help  of  an  operator.  A  fully  automated  sys¬ 
tem  for  the  point-by-point  scanning  has  been  developed  by  Erbeck  [5.25] :  Each  system  of 
Young's  fringes  is  digitized;  fringe  spacing  and  direction  are  determined  by  a  computer 
which  also  controls  the  scanning  by  the  laser  beam  (Fig.  5.9).  In  this  system,  optical 
elements  are  incorporated  that  compensate  for  the  radial  decrease  of  mean  light  intensi¬ 
ty  in  the  fringe  system  (halo)  and  for  the  formation  of  secondary  speckles. 

From  the  fringe  pattern  it  is  not  possible  to  determine  the  particle  direction  (po¬ 
sitive  or  negative) .  This  difficulty  is  overcome  by  means  of  an  additional  displacement 
of  the  recording  plate  between  the  two  exposures;  see,  e.g.  [5.26].  This  technique  of 
image  shifting  is  equivalent  to  the  frequency  shift  in  laser  Doppler  anemometry. 

As  for  any  light  sheet  method,  the  velocity  component  normal  to  the  illuminated 
plane  is  not  recovered.  If  a  particle  moving  in  this  third  direction  is  in  the  plane  of 
the  sheet  at  the  instant  of  one  of  the  two  exposures,  it  will  appear  as  a  single  image 
(or  speckles  in  the  double-exposure.  The  light  scattered  from  such  particles  contributes 
to  the  ground  noise  and  reduces  the  visibility  of  the  fringes,  like  a  speckle  at  the  edge 
of  the  circular  area  in  Fig.  5.7.  Signal  noise  or  decorrelation  is  also  generated  by  tur¬ 
bulent  fluctuations  in  the  velocity  vector,  as  well  as  by  a  nonuniform  distribution  in 
particle  size  and  concentration.  At  the  same  time,  this  noise  carries  a  certain  informa¬ 
tion,  and  attempts  have  been  made  to  decode  this  information  for  obtaining,  additionally 
to  the  measurement  of  mean  velocity,  data  on  turbulence  characteristics  [5.27]  or  particle 
size  distribution  [5.28]. 

In  addition  to  the  measurement  of  the  velocity  vector,  the  speckle  technique  allows 
for  a  direct  visualization  of  the  velocity  field.  This  analysis  is  performed  by  spatial 
filtering  of  the  developed  double-exposure  ("specklegram") .  A  schematic  representation 
of  the  set-up  for  this  reconstruction  process  is  shown  in  Fig.  5.10.  The  specklegram  is 
illuminated  with  an  expanded  laser  beam,  and  it  is  imaged  by  means  of  a  lens  in  a  plane 
("Fourier  plane")  where,  in  the  aforementioned  point-by-point  analysis,  the  pattern  of 
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Fig.  5.1o:  Arrangement  for  the  spatial  filtering  of  a  specklegram. 


Young's  fringes  had  been  observed.  An  opaque  screen  is  placed  in  this  plane  that  has  a 
transparent  hole  at  an  off-axis  location.  A  small  amount  of  the  light  can  pass  through 
the  hole,  and  a  second  lens  forms  an  image  of  the  specklegram  in  the  final  image  plane. 
Since  every  point  of  the  specklegram  is  illuminated  by  the  expanded  laser  light,  the  Fou¬ 
rier  plane  is  covered  by  an  infinite  number  of  overlapping  fringe  systems.  Light  can  pass 
through  the  hole  if  the  hole,  in  the  Fourier  plane,  is  located  on  a  bright  interference 
fringe.  It  turns  out  that  this  requirement  is  met  by  such  systems  which  are  normal  to  the 
radius  combining  the  hole  and  the  axis,  and  for  which  the  ratio  of  the  distance  from  the 
hole  to  the  optical  axis  and  the  fringe  spacing  is  an  integer  number.  As  a  consequence, 
only  such  points  appear  illuminated  in  the  (final)  image  plane  for  which  the  flow  veloci¬ 
ty  in  the  respective  object  points  of  the  flow  plane  is  of  a  certain  value  and  direction. 
The  points  illuminated  in  the  image  plane  lay  on  fringes  or  contours  of  equal  velocity 
component.  The  velocity  difference  from  fringe  to  fringe  (not  to  confuse  with  the  Young's 
fringes!)  is  constant.  The  amount  and  the  direction  of  the  velocity  component,  visualized 
by  means  of  these  "equi-velocity  fringes"  depends  on  the  location  of  the  hole  in  the  Fou¬ 
rier  plane,  see,e.g.  [5.21,  5.23,  5.29]. 


Fig.  5.11:  Equi-velocity  fringes  in  the  field  of  an  internal  gravity  wave  behind  a  cy¬ 
linder  moving  vertically  through  stratified  saltwater.  Visualization  by  spa¬ 
tial  filtering  of  the  specklegram  (Gartner  et  al.  [5.21]). 

Fig.  5.11  shows  a  system  of  equi-velocity  fringes  in  the  field  of  a  plane  internal 
gravity  wave  which  has  been  produced  by  a  cylinder  moving  upwards  through  density-strati¬ 
fied  salt-water  [5.21].  The  observable  velocity-fringes  are  at  the  same  time  curves  of 
equal  phase  in  the  wave  system.  Locating  the  hole  at  a  different  position  in  the  Fourier 
plane  results  in  visualizing  a  different  component.  This  reconstruction  by  spatial  fil¬ 
tering  suffers  from  the  low  light  intensities  that  have  to  be  managed. 


Application  of  speckle  photography  for  velocity  measurements  is  not  yet  widespread. 
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An  apparent  limitation  is  the  availability  of  an  ultra-short  and  intensive  light  source 
for  generating  the  double-exposure.  This  problem  is  the  smaller  the  lower  the  flow  velo¬ 
city  to  be  measured; this  fact  explains  successful  applications  of  the  method  to  Benard 
convection  [5.29,  5.30]  and  the  low-speed  flow  in  a  towing  tank  [5.31,  5 . 32]  . However ,  the 
general  suitability  of  speckle  velocimetry  for  gas  flows  with  velocities  being  not  very 
small  has  also  been  verified  [5.33  -  5.35]. 

5.3  LASER  INDUCED  FLUORESCENCE 


Fluorescent  tracer  materials  have  been  used  since  quite  a  number  of  years  for  enhan¬ 
ced  flow  visualization.  These  materials  emit  a  characteristic,  fluorescent  radiation  upon 
excitation  by  light  of  an  appropriate  wavelength.  The  wavelength  of  fluorescence  is  dif¬ 
ferent  from  the  wavelength  necessary  for  excitation,  which  is  typical  for  inelastic  scat¬ 
tering.  The  fluorescent  radiation,  when  suitably  analyzed,  may  contain  information  on  the 
local  density,  species  concentration,  pressure,  temperature,  and  even  velocity  (Doppler 
effect) .  When  the  exciting  laser  light  is  provided  in  form  of  a  plane  sheet,  a  selected 
plane  in  the  flow  can  be  visualized,  and  data  on  the  two-dimensional  (plane)  distribution 
of  one  or  several  of  the  mentioned  quantities  can  be  generated.  Taking  records  or  photo¬ 
graphs  through  a  narrow  optical  filter  that  passes  only  the  wavelengths  of  fluorescence 
allows  for  blocking  off  all  the  noise  generated  by  the  incident  (exciting)  light. 

Laser-induced  fluorescence  measurements  with  a  light  sheet  from  an  argon  ion  laser 
have  been  performed  by  Koochesf ahani  and  Dimotakis  [5.36]  for  determining  local  concen¬ 
tration  values  in  the  mixing  of  two  liquid  streams.  However,  the  majority  of  applications 
of  planar  laser-induced  fluorescence  is  in  the  field  of  compressible,  usually  supersonic 
gas  flow  at  relatively  low  density  levels.  The  fluorescing  trace  material,  with  which  the 
gas  is  seeded,  is  either  molecular  iodine,  [5.37  -  5.44],  or  atomic  sodium,  Na  [5.45  - 
5.47]  .  Under  the  assumption,  that  the  tracer  molecules  are  uniformly  distributed  in  the 
flowing  gas,  one  may  conclude  from  the  recorded  light  intensity  onto  the  local  value  of 
the  gas  density  in  the  thin  volume  of  the  light  sheet.  McDaniel  et  al.  [5.40]  measured 
density  fluctuations  in  a  supersonic  nozzle  flow;  Ackermann  et  al.  [5.43]  studied  the 
possibility  of  measuring  pressure  fluctuations;  the  Doppler  shift  of  the  fluorescent  sig¬ 
nal  has  been  used  for  velocity  measurements,  e.g.,  by  Zimmerman  and  Miles  [5.46],  Hiller 
and  Hanson  [5.42],  and  Hassa  et  al.  [5.44].  Temperatures  have  been  measured  in  flames 
(see  below) . 

Iodine  is  an  aggressively  reacting  tracer  material,  incompatible  with  many  metals, 
and  it  has  been  used  as  a  tracer  in  supersonic  flows  of  helium  and  nitrogen  with  which  it 
can  be  mixed  easily  at  room  temperature.  The  emitted  fluorescence  is  yellow  and  can  be 
excited  with  the  green  (514.5  nm)  line  of  an  argon  ion  laser.  The  signal  intensity  is  not 
uniquely  dependent  on  the  gas  density  (or  number  of  I 2  molecules  in  the  control  volume) , 
but  it  is  also  influenced  by  pressure  and  temperature.  McDaniel  [5.41]  found  that  the  lat¬ 
ter  dependence  can  be  suppressed  by  detuning  the  laser  from  the  absorption  frequency, 
however  with  the  disadvantage  of  decreasing  the  mean  signal  amplitude.  Compared  with  io¬ 
dine,  sodium  has  a  much  higher  signal  intensity  and  the  additional  advantage,  that  the 
molecular  mass  of  sodium  vapor  is  approximately  the  same  as  that  of  the  gases  to  which  it 
is  applied  (nitrogen  or  helium) .  But  the  seeding  of  these  working  gases  with  sodium  is  a 
much  more  complicated  procedure  than  the  seeding  with  iodine  [5.45].  In  both  cases,  see¬ 
ding  rates  in  terms  of  mass  concentration  are  in  the  order  of  1o-3  to  1o”4. 

The  radical  OH  produced  in  many  flames  is  a  reactant  which  itself  is  fluorescent.  A 
number  of  experiments  have  been  reported  in  which  OH  has  been  excited  with  radiation  wa¬ 
velengths  in  the  visible  and  near  ultraviolet  part  of  the  spectrum  [5.48  -  5.51].  Since 
the  greatest  interest  is  in  turbulent  flames,  the  excitation  should  be  performed  with  a 
pulsed  dye  laser  in  order  to  freeze  the  turbulent  fluctuations.  The  signal  intensity  of 
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the  OH  fluorescence  is  much  lower  than  the  intensity  in  the  iodine  or  sodium  experiments, 
and  it  is  therefore  necessary  to  apply  an  image  intensifying  system  whose  output  can  be 
digitized  and  displayed  on  a  monitor  (Fig.  5.12).  Values  of  the  relative  OH-concentration 
can  be  extracted  from  such  a  display.  Determination  of  the  absolute  concentration  requi¬ 
res  at  least  one  single-point  calibration  with  another  method. 
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Fig.  5.12:  Set-up  for  measuring  a  spray  flame  with  laser-induced  fluorescence  (left). 

The  beam  from  a  pulsed  Nd:YAG-laser  is  expanded  in  a  plane.  A  photodiode 
array  camera  records  the  fluorescent  scattered  signal  which  is  displayed  on 
a  monitor  (right).  (R.K.  Hanson,  Stanford  University,  Ca,  USA) 

Besides  OH,  fluorescence  from  nitric  oxide  [5.52]  and  the  radical  CH  [5.53]  have 
been  used  for  concentration  measurements. 
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6.  METHODS  BASED  ON  REFRACTIVE  INDEX  MEASUREMENTS 


The  density  of  a  fluid  is  a  function  of  the  fluid's  refractive  index.  A  number  of 
optical  methods  are  known  by  which  variations  of  the  refractive  index  in  a  transparent  me¬ 
dium  can  be  visualized  and/or  measured.  These  methods  may  be  applied  to  the  study  of  such 
flows  that  are  associated  with  density  changes  caused  by  either  compressibility,  thermal 
effects,  or  mixing  of  fluids  of  different  densities.  Very  common  is  the  use  of  these  me¬ 
thods  in  high-speed  wind  tunnels,  ballistic  ranges,  and  shock  tubes.  As  a  matter  of  fact, 
this  is  historically  the  first  application  of  such  optical  methods  for  flow  visualization. 
Independently  of  this  field  of  application  developed  the  use  of  the  methods  for  study¬ 
ing  convective  heat  transfer.  Combustion  and  plasma  flows  are  further  fields  of  applica¬ 
tion  where  the  three  causes  of  density  changes,  compressibility,  mixing,  and  high  tempe¬ 
ratures,  may  all  be  present  simultaneously. 

Refractive  index  fields  are  visualized  or  measured  by  line-  of-  sight  methods,  and 
the  information  is  integrated  along  the  path  of  the  light  through  the  transparent  medium 
(fluid) .  When  passing  through  the  field,  the  light  is  affected  with  respect  to  its  direc¬ 
tion  of  propagation  (refraction)  and  its  optical  phase.  Either  one  of  these  effects  is 
used  in  the  methods  to  be  discussed  here.  An  extensive  literature  reviewing  this  field 
of  experimentation  is  available;  see,  e.g.  [6. 1-6. 6]. 

6.1  RELATIONSHIP  BETWEEN  FLUID  DENSITY  AND  REFRACTIVE  INDEX 

An  analysis  of  interaction  of  an  electromagnetic  wave  and  the  molecules  of  a  (trans¬ 
parent)  medium  delivers  the  relationship  between  the  density  of  a  fluid,  and  its  refrac¬ 
tive  index;  see,  e.g.  [1.1].  The  general  relationship  which  applies  to  both  gases  and  li¬ 
quids,  known  as  the  Clausius-Mosotti  equation  or  the  Lorenz-Lorentz  equation,  involves  a 
number  of  molecular  constants  of  the  fluid  and  is  not  handy  to  work  with.  A  simplified 
form,  the  Gladstone-Dale  equation,  applies  to  gases: 

n  -  1  =  K  •  p  (6.1) 

Thus,  the  refractive  index  of  a  gas,n,  is  proportional  to  the  density,  p.  The  Gladstone- 
Dale  constant  K,  having  the  dimension  [p  ^  ] ,  is  a  characteristic  quantity  for  each  gas, 
and  it  changes  slightly  with  the  light  wave  length  X,  since  the  refractive  index,  too, 
is  known  to  depend  on  wavelength.  Values  of  K  for  a  number  of  gases,  and  values  of  K  for 
air  at  various  wavelenths  are  given  in  Table  6.1  and  6.2,  respectively. 
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Table  6.1 


Gladstone-Dale  constants  for  different  gases  [1.1] 


gas 

K [cm3 / g] 

wavelength  [jim] 

temperature 

He 

0.196 

0.633 

295 

Ne 

0.075 

0.633 

295 

Ar 

0.157 

0.633 

295 

Kr 

0.115 

0.633 

295 

Xe 

0.119 

0.633 

295 

H2 

1.55 

0.633 

273 

°2 

0. 19o 

0.589 

273 

N2 

0.238 

0.589 

273 

co2 

0.229 

0.589 

273 

NO 

0.221 

0.633 

295 

H20  (vapor) 

0 . 31  o 

0.633 

273 

CF4 

0.122 

0.633 

3o2 

CH 

0.617 

0.633 

295 

«6 

0.113 

0.633 

295 

Table  6 . 2 

Gladstone-Dale  constant  for  air  (at  T  =  288  K)  [1.1] 


K  [cm3 / g] 

wavelenth 

0.2239 

0.9125 

0.225O 

0 . 7o34 

0.2259 

0.6o74 

0.2274 

0 . 5o97 

0 . 23o4 

0 . 4o79 

0.2330 

0.3562 

The  quantity  (n-1)  is  called  the  refractivity  of  a  fluid.  It  is  by  orders  of  magni¬ 
tude  larger  for  liquids  than  for  gases.  Therefore,  extremely  small  density  differences 
in  a  liquid  can  produce  considerably  large  signals  in  the  optical  methods  to  be  discussed. 
The  relationship  between  the  refractive  index  and  the  density  of  water,  as  the  most  in¬ 
teresting  working  liquid,  should  be  derived  from  the  full  Clausius-Mosotti  equation.  But 
since  the  necessary  molecular  data  is  not  completely  available,  the  relationship  is  usu¬ 
ally  provided  by  a  calibration.  The  refractive  index  then  is  determined  with  an  interfero¬ 
meter  or  refractometer.  Most  common  are  the  two  following  situations:  a)  n  is  calibrated 
against  the  temperature  of  water;  a  formula  describing  accurately  this  relationship  in  the 
temperature  range  between  2o  and  34°C  has  been  given  by  Dobbins  and  Peck  [6.7] ;b)  n  is  mea¬ 
sured  for  H^O-NaCl  solutions  against  the  salinity  at  constant  temperature.  These  solutions 
are  of  interest  for  performing  stratified  flow  experiments.  The  relationship  is  found  to 
be  of  the  form 

“saltwater  =  no  +  const  ’  p'  (6.2) 

p  being  the  density  of  the  solution  and  nQ  a  reference  value  of  the  refractive  index  of 
saltwater;  [6.8,  6.9] 


The  K  values  in  Table  6.1  are  given  for  pure  gases.  If  the  working  gas  is  a  mixture 
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of  several  components  with  constant  and  known  concentration  (e.g.  air) ,  the  Gladstone- 
Dale  constant  of  the  mixture  can  be  calculated  from  the  values  of  the  pure  components. 

Let  n,  K,  p  be  the  refractive  index,  Gladstone-Dale  constant,  and  density  of  the  mixture, 
respectively;  then,  eq.  (6.1)  applies  to  the  mixture,  and  K  can  be  expressed  by  the values 
(of  the  i-th  component)  through 


IN  p  ,  JN 

=  l  k  .  —  =  y 

.  t  +  i  p  .  L 


i=1 


i=1 


K.  c. 

l  l 


(6.3) 


where  the  mixture  is  assumed  to  consist  of  N  components  of  partial  density  p^  or  relative 
mass  concentration  c^?  i.e.  =  p  and  £c^  =  1.  Within  sufficient  accuracy,  the  Glad¬ 

stone-Dale  constant  of  air  (Table  6.2)  can  be  determined  from  the  respective  K-values 
for  oxygen  and  nitrogen  (Table  6.1). 


An  example  of  a  gas  mixture  with  unknown  and  local  and/or  temporal  variation  of  the 
concentration  values  of  the  components  is  a  flame.  This  means,  that  the  density  of  the 
mixture,  p,  cannot  be  determined  from  a  single  optical  measurement  of  the  refractive  in¬ 
dex  n,  because  the  Gladstone-Dale  constant  of  the  mixture,  K,  is  also  variable.  Or  in 
other  words:  the  number  of  unknowns  is  more  that  one  (it  is  exactly  N,  namely  p  and  (N-1) 
values  of  c^,  see  eq.  (6.3)) ,  and  the  information  from  one  measurement  of  n  is  not  suffi¬ 
cient  for  determining  the  N  unknowns.  A  possible  approach  to  this  problem  is  to  make  use 
of  the  dispersive  nature  of  n  and  K,  i.e.,  optical  measurements  would  have  to  be  perfor¬ 
med  with  several  wavelenths  X  ^ .  The  number  of  wavelengths  used  should  then  be  equal  to 
N,  and  this  would  yield  a  set  of  N  equations  of  the  form 


n  (X  j ) 


i=1 


W  ci  p 


(6.4) 


Unfortunately,  the  K^-values  for  most  gaseous  components  are  so  weakly  dispersive,  that 
the  optical  signals  for  measuring  n,  taken  at  different  wavelengths  in  the  visible  range, 
differ  not  enough  for  providing  independent  information. 


The  principle  of  taking  multiple-wavelengths  optical  data  has  been  successfully 
applied  to  plasma  flows,  i.e.,  to  the  flow  of  an  ionized  gas.  If  one  considers  a  mona¬ 
tomic  gas,  the  ionized  gas  is  a  mixture  consisting  of  neutral  atoms,  ions,  and  free 
electrons.  The  degree  of  ionization,  a^,  is  the  (relative)  mass  concentration  of  ions  in 
the  mixture.  Due  to  the  little  mass  of  an  electron,  the  mass  concentration  of  the  elec¬ 
tron  gas  component  can  be  neglected.  The  Gladstone-Dale  equation  for  this  mixture  can  be 
written  in  the  form 


n-1  =  p[(1-ax)KA  +  o^Kj]  +  NeK^  ,  (6.5) 

where  and  are  the  Gladstone-Dale  constants  of  the  atom  gas  and  the  ion  gas,  respec¬ 
tively;  is  the  electron  number  density,  i.e.,  the  number  of  electrons  per  unit  volume, 
and  the  Gladstone-Dale  constant  of  the  electron  gas,  ,  has  the  dimension  cm3,  different 

from  the  dimension  of  K  and  K  .  For  argon,  a  gas  frequently  used  in  plasma  experiments, 

A  •**  o 

Alpher  and  White  [6.1o]  have  found  =  -j  K^.  The  value  of  the  refractive  index  of  the 
electron  gas,  ng,  is  determined  in  plasma  theory,  and  it  is  found  to  be  proportional  to 
the  square  of  the  light  wavelength.  If  one  measure  X  in  cm,  the  Gladstone-Dale  equation 
for  the  electron  gas  can  be  written  as 

n-1  =  -  4.46-1o"14-  X2 -N 
e  e 


(6.6) 
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Thus,  K'  is  highly  dispersive.  The  contribution  of  at  high  enough  ionization  degrees, 
causes  the  refractive  index  n  of  the  ionized  gas  to  depend  strongly  on  the  wavelength. 
Measuring  the  optical  signal  at  three  different  wavelengths  X^  (see  eq.  (6.4))  therefore 
enables  one  to  determine  the  three  unknowns  in  eq.  (6.5),  p,  ,  and  Ne;  see,  e.g., 

[6.11]  . 


At  a  given  wavelength,  the  refractivity  (n-1)  per  electron  is  by  an  order  of  magni¬ 
tude  greater  than  that  per  atom.  The  optical  behavior  of  an  ionized  gas  is  therefore  do¬ 
minated  by  the  presence  of  free  electrons,  even  at  a  relatively  low  ionization  level. 
Since  the  Gladstone-Dale  constant  of  the  electron  gas  is  negative,  the  optical  response 
to  density  changes  is  opposite  to  the  respective  optical  response  caused  by  the  neutral 
or  the  atom  gas.  Ionization  is  often  produced  through  a  strong  shock  wave.  The  production 
of  free  electrons  and  the  simultaneous  increase  of  the  neutral  gas  density  cause  optical 
signals  of  opposite  sign,  in  the  worst  case  just  cancelling  one  another. 

A  dissociated  diatomic  gas  is  a  mixture  of  (neutral)  molecules  and  free  atoms. With 
being  the  degree  of  dissociation,  the  Gladstone-Dale  equation  for  a  dissociated  gas  is 

n-1  =  p[(1-aD)KM  +  <*dKa3  /  (6.7) 

where  and  K.  are  the  Gladstone-Dale  constants  of  the  molecule  and  the  atom  gas,  re- 
M  A 

spectively.  Values  of  for  oxygen  and  nitrogen  have  been  measured  in  shock  tube  experi¬ 
ments  [6.12,  6.13]. 

The  derivation  of  the  Gladstone-Dale  equation  (6.1)  involves  a  number  of  simplify¬ 
ing  assumption.  One  of  then  is  that  the  wavelength  of  the  light  in  the  measuring  system 
is  far  from  an  absorption  wavelength  of  the  test  gas.  This  is  usually  the  case  for  "nor¬ 
mal"  gases  and  at  "normal"  thermodynamic  conditions.  (For  this  reason,  the  gas  is  invi¬ 
sible.)  However,  the  assumption  might  be  violated  for  some  "unusual"  gases  or  vapors  and 
at  high  temperatures  when  the  gas  is  in  an  excited  state.  Then,  absorption  wavelengths 
may  exist  in  the  visible  range  of  the  spectrum.  Close  to  such  an  absorption  wavelength, 
the  gas  exhibits  anomalous  dispersive  behavior  with  the  refractivity  (n-1)  being  many  or¬ 
ders  of  magnitude  greater  than  in  the  nonresonant  regime.  The  result  is  that  the  same  ab¬ 
solute  density  change  Ap  will  produce  much  stronger  alterations  in  the  refractive  index 
of  the  gas,  An,  or  cause  much  higher  signal  amplitudes.  Enhancement  of  flow  visualization 
near  resonant  or  absorption  wavelengths  has  been  verified  by  seeding  a  gas  with  sodium 
vapor  and  operating  the  light  source  in  the  optical  system  with  the  sodium  D^  wavelength 
[6.14,  6.15] . 

6.2  DEPENDENCE  OF  THE  VISIBLE  SIGNAL  ON  THE  REFRACTIVE  INDEX  DISTRIBUTION 

A  light  ray  tensmitted  through  a  fluid  flow  of  variable  refractive  index  is  affec¬ 
ted  in  two  different  ways:  The  ray  is  deflected  from  its  original  direction  (this  is  usu¬ 
ally  understood  as  "refraction") ,  and  the  wave  describing  this  light  ray  or  beam  is  alte¬ 
red  with  respect  to  its  optical  phase.  Both  quantities,  light  deflection  and  optical  pha¬ 
se  change,  are  integrated  quantities,  i.e.,  they  can  be  measured  after  the  light  is  trans 
mitted  through  the  flow,  and  nothing  follows  from  these  measurable  quantities  on  how  the 
light  wave  was  affected  at  a  particular  position  inside  the  flow.  The  optical  methods  to 
be  discussed  depend  on  one  of  the  two  measurable  effects.  All  interferometric  methods  vi¬ 
sualize  differences  in  the  optical  phases  of  two  (or  more)  light  waves.  A  number  of  me¬ 
thods  are  available  whose  signal  is  a  measure  of  the  integrated  value  of  the  light  deflec 
tion  in  the  flow  field.  Most  of  the  latter  methods  provide  qualitative  results,  i.e.  flow 
picture,  while  the  interferometric  methods  also  allow  for  a  quantitative  determination  of 
the  refractive  index  or  density  distribution. 
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The  dependence  of  the  signal  of  each  method  on  the  refractive  index  distribution 
must  be  found  from  a  detailed  analysis  of  the  optical  path  in  the  flow  {see,  e.g.,  [1.1], 

also  for  further  references).  Tabel  6.3  lists  the  classes  of  optical  methods,  the  types 
of  signals  available  from  a  record  (photograph) ,  and  the  dependence  of  this  signal  on  the 
refractive  index  distribution  n(x,y,z) .  It  is  assumed  in  the  table  that  only  the  light  de¬ 
flection  or  the  change  of  n  in  the  y-direction  contribute  to  the  signal.  The  incident 
light  is  assumed  to  be  parallel  to  the  z-direction  (Fig.  6.1).  It  is  seen  that  the  methods 
are  grouped  in  the  table  in  a  systematic  way:  The  shadowgraph  depends  on  the  second  deri¬ 
vative  of  n;  schlieren,  moire,  speckle  and  schlieren  interference  technique  on  the  first 
derivative;  and  reference  interferometry  on  the  absolute  value  (zeroth  derivative)  of  n. 


Table  6 . 3 


Optical  methods  for  flow  visualization,  signal  pattern,  and  dependence  on  the  distribu¬ 
tion  of  the  refractive  index  n 


method  measured  signal  dependence  on  n 


shadowgraph 

intensity  constrast 

a  •  / 

zi 

Z  2 

32n  , 

w dz 

schlieren 

moire 

speckle  photography 

schlieren  interferometer 

intensity  constrast 

moire  fringe  shift 

displacement  of  speckles 

fringe  order  or 

fringe  shift 

z2 

; 

zi 

3n  j 

3y  dz 

reference  beam 

fringe  order  or 

interferometer 

fringe  shift 

z2 

; 

n  dz 

phase  contrast 

fringe  order 

zi 

incident  cross  section  of  plane  of 

light  test  flow  observation 


Fig.  6.1:  Deflection  of  a  light  ray  in  an  inhomogeneous  density  field,  z.,  z2  are  en¬ 
trance  and  exit  coordinates  of  the  ray  in  the  flow  field.  The  ray  is  de¬ 
flected  by  an  angle  e. 

In  the  general  case,  the  refractive  index  is  a  function  of  all  three  space  coordi¬ 
nates:  n  =  n(x,y,z).  As  shown  in  the  3rd  column  of  Table  6.3,  the  dependence  of  n  on  z 
is  not  recovered  because  of  the  integration  from  z^  to  z2  in  the  z-direction.  Here,  z1 
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is  the  coordinate  where  a  particular  ray  enters  the  flow  field,  and  z 2  is  its  exit  coor¬ 
dinate.  The  refractive  index  in  the  medium  outside  of  the  test  field  is  assumed  to  be 
homogeneous.  Determining  the  three-dimensional  distribution  n  (x,y,z)  requires  to  take 
more  than  one  record  of  the  visible  signal,  each  under  a  different  viewing  angle.  Such 
a  procedure  and  the  involved  computer-evaluation  of  the  data  is  known  as  tomography. 

The  results  in  the  third  column  of  Table  6.3  are  derived  in  the  respective  theory 
under  the  simplifying  assumption  that  the  light  propagates  along  a  straight  path,  parall¬ 
el  to  the  z-direction  through  the  flow.  This  is  somewhat  contradictory  to  the  fundamental 
principle  of  those  methods  which  are  based  on  the  deflection  of  light  in  the  test  field. 
The  theory  takes  into  account  the  finite  curvature  of  the  light  rays  when  they  leave  the 
field,  but  it  neglects  the  differences  in  the  heights  y  of  the  entrance  and  exit  position 
(Fig.  6.1).  The  assumption  is  very  reasonable  for  the  majority  of  experiments  performed 
in  gas  flows.  Violations  of  this  assumption,  referred  to  as  “strong  refraction  effects", 
occur  in  the  presence  of  large  gradients  of  the  refractive  index  or  density,  and  in  li¬ 
quids  where  the  refractivity  ( n— 1 )  is  not  very  small  when  compared  with  one. 

In  the  case  of  strong  refraction  one  is  confronted  with  the  problem  that  not  only 
the  refractive  index  along  the  light  path  is  unknown,  but  also  the  (geometrical)  light 
path  itself.  The  problem  has  been  investigated  for  a  number  of  special  cases  and  with  a 
two-dimensional  refractive  index  field,  n  =  n(x,y) .  Most  of  the  approaches  are  numerical 
studies  that  demonstrate,  by  means  of  optical  ray  tracing,  the  difference  between  strong 
refraction  and  the  refraction-free  situation.  The  thermal  boundary  layer  close  to  a  hea¬ 
ted,  vertical  flat  plate  is  a  typical  example  of  such  a  plane  flow  with  strong  refraction 
in  one  direction  due  to  a  large  temperature  gradient  normal  to  the  plate.  The  refractive 
bending  of  the  light  may  be  studied  in  a  plane  x  =  const  (the  plane  of  Fig.  6.1)  .  With  an 
exact  knowledge  of  the  distribution  of  the  refractive  index  in  this  plane,  n  (y) ,  one 
could  calculate  the  trajectory  of  a  light  ray  through  the  flow.  If  one  can  assume,  that  n 
is  distributed  linearly  in  this  plane,  it  follows  that  the  light  trajectory  is  a  parabola; 
see,  e.g.  [6.16  -  6.18].  The  refractive  index  field  can  then  be  reconstruced  from  the  in- 
terferogram  by  making  use  of  this  "correction  parabola" . 

Due  to  the  strong  refraction  effects,  the  image  of  the  test  field  in  the  recording 
plane  is  distorted.  This  distortion  can  be  minimized  by  focusing  onto  the  recording  plane 
a  plane  z  =  const  in  the  test  field  situated  between  the  center  and  a  position  at  one- 
third  the  field  width  away  from  the  exit  plane  into  the  test  section  [6.19]. 

For  the  purpose  of  a  comparison  of  the  optical  signal,  which  can  be  obtained  with 
the  various  methods,  Figs.  6.2  -  6.6  present  visualizations  of  the  same  test  object  by 
■different  methods.  The  flow  under  study  is  the  Rayleigh-Benard  convection  in  a  rectangu¬ 
lar  enclosure  with  a  slight  temperature  difference  (AT  ~  2K)  between  top  and  bottom. plate . 
The  test  fluid  is  silicon  oil.  Due  to  the  high  refractivity  of  the  liquid,  the  shadowgraph 
exhibits  a  considerable  optical  distortion.  The  flow  cells  are  clearly  distinguishable. 
Since  the  test  field  is  two-dimensional,  the  fringes  in  the  Mach-Zehnder  interf erogram 
are  curves  of  equal  temperature.  The  fringe  contours  in  the  schlieren  interf erogram  of 
Fig.  6.5  coincide  with  the  streamlines.  This  applies  only  to  the  special  adjustment  of 
the  interferometer  in  this  figure,  as  it  has  been  verified  by  Buhler  [6.20].  A  different 
adjustment  of  this  instrument  (Fig.  6.6)  results  in  a  completely  different  pattern.  The 
generation  of  the  patterns  visible  in  these  figures  will  be  explained  in  the  following 
chapters . 
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Fig.  6.2:  Shadowgraph  of  the  Ray- 
leigh-Benard  flow  cells 

in  a  rectangular  enclo¬ 
sure  filled  with  silicon 
oil.  Temperature  diffe¬ 
rence  between  top  and 
bottom  plate  is  about  2K 
The  photographic  imaging 
is  distorted  due  to  the 
strong  refraction  in  the 
test  field. 


Fig.  6.3:  Schlieren  record  of 

the  flow  shown  in  Fig. 
6.2.  Original  is  taken 
in  a  color  schlieren 
system. 


Fig.  6.4:  Mach-Zehnder  interfero- 
gram  of  the  flow  shown 
in  Fig.  6.2.  The  inter¬ 
ferometer  was  aligned 
to  infinite  fringe 
width.  The  fringes  are 
contours  of  constant 
temperature . 
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Fig.  6.5:  Schlieren  interfero- 

gram  taken  with  a  Woll¬ 
aston  prism  interfero¬ 
meter  of  the  flow  shown 
in  Fig.  6.2.  The  inter¬ 
ferometer  was  aligned 
to  finite  fringe  width 
with  the  undisturbed 
fringes  being  vertical. 


Fig.  6.6:  Schlieren  interfero- 
gram  taken  with  a 
Wollaston  prism  in¬ 
terferometer  of  the 
flow  shown  in  Fig. 

6.2.  The  interfero¬ 
meter  was  aligned  to 
infinite  fringe  width 
at  horizontal  undis¬ 
turbed  fringes. 
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7 .  SHADOWGRAPH 


The  shadowgraph,  scientificly  explored  first  by  Dvorak  [7.1],  is  the  simplest  of 
the  optical  visualization  methods.  It  does  not  require  any  optical  element  excdpt  a  light 
source,  and  the  shadow  effect  produced  by  refractive  index  fields  can  be  observed,  there¬ 
fore,  outside  a  laboratory  in  the  open  air  where  the  sun  serves  as  the  light  source.  Ex¬ 
amples  are  the  shadow  patterns  of  raising  warm  air  or  benzene  vapor,  projected  by  the  sun 
light  onto  a  rigid  surface.  Even  projections  of  shock  waves  have  been  observed  in  the  open 
air,  as  will  be  shown  below. 

In  a  laboratory  experiment,  either  a  divergent  or  a  parallel  beam  of  light  is  trans¬ 
mitted  through  the  test  section  of  the  flow  facility  (Fig.  7.1) .  The  shadow  pattern  is  ob¬ 
served  in  a  plane,  normal  to  the  optical  axis,  at  a  distance  l  behind  the  flow  field.  Most 
essential,  as  for  any  optical  visualization  method,  is  the  use  of  a  point  source.  In  order 
to  a.yoid  noise  by  diffraction,  a  thermal  light  source  should  be  used  and  not  a  laser  source 


3 


5 
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Fig.  7.1:  Shadowgraph  system  without  (a)  and  with  (b)  camera;  1:  point  light  source; 

2:  optical  axis;  3:  collimating  lenses;  4:  cross  section  of  flow  facility 
test  section;  5:  recording  plane  (photographic  film);  6:  camera  objective. 

It  is  of  interest  to  visualize  large  fields  of  view.  Because  it  is  then  impractical  to 
record  the  shadow  on  a  respectively  large  photographic  plate  (Fig.  7.1a) ,  a  camera  is 
focused  onto  the  position  (plane)  at  distance  £  from  the  test  field  (Fig.  7.1b) .  The 
choice  of  the  value  of  i  depends  on  the  flow  under  study.  The  sensitivity  of  the  method 
with  respect  to  resolving  small  density  changes  increases  with  £;  however,  edges  of  ri¬ 
gid  bodies  in  the  test  field,  e.g.  test  models,  will  be  the  more  out  of  focus,  the  lar¬ 
ger  the  value  of  £,  so  that  it  becomes  often  necessary  to  find  a  compromise  between  opti¬ 
cal  sensitivity  and  local  resolution. 


Visualization  of  a  flow,  which  is  invisible  when  viewed  without  any  optical  aid, 
means  that  the  plane  of  observation  (recording  plane)  contains  an  information  which  is  per¬ 
ceivable  for  the  eye.  Since  the  eye  can  perceive  information  only  in  the  form  of  changes 
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of  light  amplitude  or  intensity,  the  optical  analysis  of  the  shadowgraph  determines  the 
relative  changes  of  light  intensity,  AI/I,  which  exist  in  the  recording  plane  due  to  the 
shadow  effect.  I  is  the  uniform  light  intensity  for  a  homogeneous  refractive  index  field, 
and  AI  are  the  intensity  changes  caused  by  variation  of  the  refractive  index  n  in  the 
flow.  The  result  for  AI/I  (also  called  the  contrast)  is  (see,  e.g.  [1.1]): 


AI (x , y) 
I 


i 


9 

^  3 2  3 2 

•  f  +  -mr)  In  n(x,y,z)  dz . 

Z1 


(7.1) 


The  information,  here:  the  contrast  AI/I,  is  recorded  in  a  x-y  plane.'  The  dependence  of  n 
on  z  is  not  recovered  due  to  the  integration  in  z-direction. 

The  most  dramatic  change  in  the  second  derivative  of  n  or  the  density  p  occurs  in 
a  shock  wave.  The  shadowgraph  is  most  appropriate  for  visualizing  the  geometry  of  shock 
patterns,  either  a  steady  configuration  in  supersonic  flow  (Fig.  7.2) ,  or  unsteady  shock 
configurations  produced  in  a  shock  tube  (Fig.  7.3).  The  shock  appears  as  a  dark  contour, 


Fig.  7.2:  Shadowgraph  of  a  projectile  flying  at  supersonic  speed.  Clearly  seen  are  the 

formation  of  shock  waves  and  the  turbulent  wake  behind  the  projectile  (German- 
French  Research  Institute  ISL,  St.  Louis,  France) . 

the  "shadow"  of  the  shock.  No  light  is  arriving  at  this  contour  line  because  it  has  been 
deflected  due  to  the  strong  density  change.  In  many  cases  the  deflected  light  is  regained 
in  concentrated  form  as  a  bright  band  on  the  downstream  side  of  the  dark  shock  contour. 
Shadowgraphs  are  used  for  measuring  the  stand-off  distance  of  the  bow  shock  at  a  blunt  bo¬ 
dy  in  supersonic  flow.  For  such  a  measurement  it  is  important  to  know  that  the  exact  lo¬ 
cation  of  the  infinitesimally  thin  shock  surface  is  the  upstream  edge  of  the  visible  dark 
band.  While  in  conventional  systems  the  recording  of  these  data  is  performed  by  photogra¬ 
phic  means  [7.2],  modern  systems  use  CCD  cameras  or  linear  arrays  of  photodiodes  for  this 
purpose  [7.3,  7.4], 

Projections  by  the  sunlight  of  a  shock  wave  onto  a  surface  can  be  observed  in  the 
open  air.  The  shadow  of  the  shock  waves  giving  rise  to  the  buffeting  of  an  airplane  at 
high  subsonic  Mach  numbers  can  be  seen  from  the  inside  of  the  plane  on  the  wing  surface 
if  the  sun  is  at  an  appropriate  position  [7.5].  The  projection  of  the  spherical  blast 
shock  generated  by  a  large  scale  explosion  is  shown  in  Fig.  7.4.  This  photograph  has  been 
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Fig.  7.3:  Series  of  shadowgraphs,  taken  at 

different  instants  of  time,  of  the 
propagation  of  an  unsteady  shock 
wave  in  a  bending  channel.  The  di¬ 
rection  of  the  incident  shock  wave 
is  from  left  to  right.  This  shock 
wave  is  reflected  from  the  oblique 
channel  walls  and  diffracted  around 
the  convex  corner,  where  the  forma¬ 
tion  of  a  vortical  flow  structure  is 
observed  (W.  Heilig,  Ernst-Mach-In- 
stitut,  Freiburg,  F.R.  Germany) . 
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Fig.  7.4.:  Projection  by  sunlight  of  a  hemispherical  shock  wave  onto  the  ground  surface, 
as  seen  from  an  airplane.  The  shock  originated  from  a  5oo  tons  TNT  hemispheri¬ 
cal  explosion  (J.M.  Dewey,  University  of  Victoria,  B.C.,  Canada) 


taken  from  an  airplane.  The  center  of  the  explosion  is  in  the  center  of  the  circular  area 
of  which  only  a  part  is  included  in  the  photograph.  Since  the  projection  by  the  sunlight 
onto  the  surface  is  oblique,  the  shadow  of  the  shock  is  an  ellipse,  and  it  is  not  the 
true,  circular  position  of  the  shock  on  the  surface.  The  true  position  is  also  visible 
due  to  a  compression  and  respective  change  of  the  surface  (soil)  structure. 

The  light  rays  bended  and  deflected  in  the  shadowgraph  system  may  overlap  or  inter¬ 
sect  before  reaching  the  recording  plane,  sometimes  forming  a  caustic.  The  density  profile 
in  a  laminar  compressible  boundary  layer  can  cause  such  a  concentration  of  light,  forming 
a  bright  band  parallel  to  the  wall  (Fig.  7.5).  Shadowgraphs  taken  of  turbulent  flow  fields 


Fig.  7.5: 
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Bending  of  light  rays  as  caused  by  the  density  profile  in  a  compressible 
boundary  layer.  The  resulting  light  intensity  in  the  recording  plane  has  a 
maximum;  correspondingly,  the  shadowgraph  will  exhibit  a  bright  band  parallel 
to  the  wall. 


with  short  enough  an  exposure  time  show  a  granular  pattern,  which  carries  information  on 
the  structure  of  turbulence  in  the  fluctuating  density  field.  This  pattern  which,  again, 
is  the  result  of  an  integration  along  the  light  path,  can  be  analyzed  with  statistical 
methods.  Uberoi  and  Kovsznay  [7.6]  have  devised  a  way  which  allows,  in  principle,  for  de- 
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riving  the  statistics  of  the  three-dimensional  density  field  from  the  statistical  analy¬ 
sis  of  the  two-dimensional  shadow  pattern,  provided  that  the  turbulence  can  be  assumed  to 
be  isotropically  distributed.  Although  attempts  have  been  made  to  quantify  this  process 
by  digital  image  processing  [7.7],  the  procedure  is  rather  impractical,  and  the  analysis 
of  this  problem  by  speckle  photography  (see  chapter  9)  has  been  shown  to  deliver  more  re¬ 
liable  results. 

In  summary,  the  shadowgraph  is  not  a  method  suitable  for  quantitative  measurements 
of  the  fluid  density.  However,  it  is  a  convenient  method  for  obtaining  a  quick  survey  of 
flow  fields  with  varying  density,  particularly  shock  waves. 
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8.  SCHLIEREN  SYSTEMS 


With  only  a  slight  modification  of  its  optical  arrangement,  the  shadowgraph  can  be 
converted  into  a  system  which  is  much  more  sensitve  with  respect  to  density  changes:  the 
schlieren  system.  It  became  widely  known  through  the  detailed  discussion  in  Toepler 1 s 
book  of  1864  [8.1],  and  it  has  since  then  be  named  the  Toepler  method.  A  great  variety  of 
modifications  have  been  described,  the  newest  developments  have  continuously  been  surve¬ 
yed  [8.2  -  8.8].  Among  these  reviews,  the  survey  by  Schardin  [8.3]  has  the  reputation  to 
be  the  classical  reference  work,  particularly  in  the  fields  of  experimental  aerodynamics, 
gasdynamics,  and  ballistics. 

A  schlieren  system  with  parallel  light  through  the  test  field  is  considered  (Fig. 
8.1) .  The  parallel  light  beam  is  made  convergent  by  a  lens  or  a  spherical  (parabolic) 
mirror  called  the  "schlieren  head".  An  image  of  the  light  source  is  formed  in  the  focal 
point  of  the  schlieren  head.  At  this  position,  a  knife  edge  (perpendicular  to  the  plane 


schlieren 


Fig.  8.1:  Toepler  schlieren  system  with  parallel  light  passing  through  the  test  section 
of  the  flow  facility. 

of  Fig.  8.1)  is  cutting  off  a  certain  portion  of  the  light  source  image,  thereby  reducing 
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the  intensity  with  which  the  recording  plane  is  illuminated.  The  camera  lens  focuses  a 
plane  in  the  flow  field  onto  the  recording  plane,  so  that  shadow  effects  are  eleminated. 

An  essential  requirement  for  the  schlieren  system  is  the  use  of  a  point  light  sour¬ 
ce  or  a  slit  source,  the  slit  being  parallel  to  the  knife  edge.  The  image  of  the  point 
source  is  a  small  circle,  of  which  only  a  part  of  height  a  is  allowed  to  pass  over  the 
knife  edge  (Fig.  8.2).  With  a  homogeneous  test  field,  the  recording  plane  is  evenly  illu- 


disturbed  image 
of  source  slit 


Fig.  8.2:  Image  of  the  light  source  in  the  plane  of  the  knife  edge. 

minated  with  an  intensity  I(x,y)  =  const,  which  is  the  lower,  the  smaller  the  value  of  a. 
Light  rays  deflected  by  an  angle  e  due  to  a  disturbance  in  the  test  field  cause  a  (ver¬ 
tical)  shift  of  the  light  source  image  by  an  amount  Aa: 

Aa  =  f  •  tan  ,  (8.1) 

with  f  being  the  focal  length  of  the  schlieren  head.  The  corresponding  image  points  in  the 
recording  plane  receive  an  intensity  changed  by  AI;  the  relative  intensity  change  (con¬ 
trast)  is 


AI  Aa  f 

=  —  =  —  •  tan  e 

I  a  a  y 


(8.2) 


and  from  the  respective  optical  analysis,  which  relates  the  deflection  angle  to  the  re¬ 
fractive  index  variation  in  the  test  field,  it  follows  that 


AI (x,y)  _  f  ,  1  9n  „ 

I  “  a  ‘  z{  n  Sy  dz  ' 

which  reduces  for  a  gas  flow  to 


(8.3) 


AI(x,y)  _  K.f  ,  3d  . 

—i  — r  *  a? dz 


(8.4) 


If  the  knife  edge  is  turned  by  9o°,  the  deflection  angle  ex  will  be  measured,  and  the  re¬ 
corded  relative  intensity  change  then  is 


AI (x ,y)  _  f  ,  2  1  3n  „ 

I  "  i  J  n  TTx  dz 

~  Z1 


(8.5) 


This  brief  analysis  of  the  Toepler  system  shows  that  a  schlieren  photograph  exhibits 
changes  of  the  relative  illumination  intensity  which  are  a  measure  of  the  refractive  in¬ 
dex  or  density  gradient  in  the  direction  normal  to  the  knife  edge.  Determining  the  two 
,  do 

components  -g-j  and  -g^,  requires  to  take  two  records  with  the  knife  edge  turned  by  9o°  be¬ 
tween  the  two  exposures. 
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For  a  given  density  field,  the  relative  change  in  light  intensity  or  the  contrast 
in  the  photographic  record  is  the  larger,  the  smaller  the  ratio  (a/f) .  If  one  assumes, 
that  a  relative  intensity  change  of  Al/I  =  0.1  is  still  detectable,  the  smallest  deflec¬ 
tion  angle  that  can  be  measured  is  e  .  =0.1  (a/f).  While  the  focal  length  of  the  schlie- 

ren  head,  f,  usually  is  a  given  quantity,  one  could  conclude  that  the  sensitivity  of  the 
system  could  be  increased  without  limitation  by  decreasing  the  aperture  a.  There  are,  how¬ 
ever,  mainly  three  reasons  which  limit  the  possibility  of  reducing  the  value  of  as  1.  the 
overall  illumination  intensity  I(x,y)  is  reduced  with  a,  and  the  mean  signal  amplitude 
might  fall  below  a  minimum  determined,  e.g.,  by  the  film  speed;  2.  it  is  desirable  to  mea¬ 
sure  positive  and  negative  values  of  e  with  the  same  resolution;  the  greatest  possible 
intensity  change  for  negative  e  (Fig.  8.2)  is  absolute  extinction,  i.e.  AI  =  -I  for  e  = 

-e  ,  and  (a/f)  .  <  e  ;  3.  diffraction  effects  become  effective  at  small  values  of  a, 

thereby  blurring  the  image  of  the  test  field.  Since  diffraction  noise  increases  with  the 
coherence  of  the  light  source,  it  follows  that  a  laser  is  not  an  appropriate  source  for  a 
schlieren  system.  The  sensitivity  of  the  Toepler  system,  that  obviously  depends  on  a  num¬ 
ber  of  different  physical  phenomena,  has  extensively  been  studied  in  the  past  [8.9-8.14]. 


Fig.  8.3:  Toepler  schlieren  photographs  of  the  unstable  jet  of  a  foreign  gas  (helium) 
exhausting  in  the  ambient  air. 


The  black-and-white  schlieren  photograph  (Fig.  8.3)  can,  in  principle,  be  evaluated 
for  determining  the  distribution  of  the  deflection  angle,  e(x,y)  or,  after  a  respective 
integration,  the  distribution  of  the  fluid  density  p(x,y).  This  evaluation  process  requi¬ 
res  measuring  the  shades  of  grey  in  quantitative  form. 'A  serious  source  of  errors  is  in¬ 
troduced  by  the  photographic  recording  and  development.  The  use  of  a  television  camera 
for  the  recording  and  the  subsequent  digitization  of  the  image  can  make  this  evaluation 
independent  from  subjective  measurements;  see,  e.g.,  [8.15].  However,  the  determination 

of  gradual  intensity  changes  or  contrast  remains  an  unreliable  method,  and  schlieren  sys¬ 
tems  are  primarily  used  for  qualitative  investigations  of  density  fields,  but  with  a 
much  higher  degree  of  resolution  than  the  shadowgraph.  As  will  be  shown  in  the  following 
sections,  a  schlieren  setup  can  easily  be  converted  into  systems  being  more  appropriate 
for  quantitative  measurements  of  9n/9y. 


A  number  of  optical  errors  are  associated  with  the  use  of  lenses.  These  errors  af¬ 
fect  the  quality  of  the  light  source  image  in  the  focal  point  of  the  schlieren  head,  and 
therefore  the  precision  of  the  schlieren  effect.  Astigmatism  can  be  corrected  by  a  system 
of  piano-cylindrical  lenses  [8.16,  8.17].  Coma  is  avoided  if,  in  an  off^-axis  configuration 
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with  spherical  or  parabolic  mirrors,  the  optical  axis  forms  a  Z  (Fig.  8.4) 

light  source 


Fig.  8.4:  Z-conf iguration  of  a  schlieren  system  with  spherical  (or  parabolic)  mirrors 

In  a  "double-pass  schlieren  system" ,  which  serves  for  increasing  the  sensitivity 
towards  small  density  changes,  a  conical  light  beam  passes  twice  through  the  test  field. 
Such  a  system  is  often  applied  in  ballistic  ranges  to  study  the  flaw  around  hyperveloci¬ 
ty  projectiles  under  extreme  low  density  conditions  [8.18,  8.19].  With  the  use  of  a  beam 
splitter,  as  shown  in  Fig.  8.5,  the  system  is  "on-axis",  and  one  avoids  almost  totally 
coma  and  astigmatism. 

recording  plane 


Fig.  8.5:  Double-pass  schlieren  system  using  a  beam  splitter 

The  manipulation  of  the  light  in  the  focal  point  of  the  schlieren  head  could  be  con¬ 
sidered  as  the  essential  11  schlieren  effect".  Many  alterations  of  the  Toepler  method,  which 
uses  the  knife  edge  for  such  a  manipulation,  have  been  proposed  and  realized.  In  the  pro¬ 
posed  systems  the  light  wave  is  manipulated  with  respect  to  its  intensity  (like  in  the 
Toepler  system),  or  its  phase  ("dephasing  schlieren  systems"),  or  its  color  ("color  schlie¬ 
ren  systems") .  The  aim  of  these  methods  can  be:  increase  of  the  optical  sensitivity  in  re¬ 
solving  small  density  changes;  visualization  of  the  deflection  in  different  directions; 
providing  a  signal  which  is  more  suitable  for  a  quantitative  measurement. 


In  the  "dark  field"  system,  the  conventional  knife  edge  is  replaced  by  an  opaque 
diaphragm  or  stop  cutting  off  symmetrically  all  but  a  very  small  fraction  of  the  undis¬ 
turbed  light  source  image.  The  undisturbed  field  appears  with  a  very  reduced  intensity 
(dark  field) .  When  a  point  source  is  used,  the  schlieren  stop  has  to  be  a  small  circular 
diaphragm  (Fig.  8.6).  All  deflection  angles  of  the  same  absolute  value  produce  the  same 
contrast,  regardless  of  their  orientation  in  the  x-y  plane.  Dark  field  systems  have  been 
studied  extensively  by  Wolter  [8.5] .  The  opaque  schlieren  stop  can  be  made  by  taking  pho¬ 
tographs  of  the  light  source  image  at  different  exposure  times  and  then  using  the  appro¬ 
priate  negative  [8.15].  Since  the  sharp  edge  of  the  stop  enhances  the  diffractive  noise, 
it  has  been  proposed  to  employ  a  stop  with  a  gradual  variation  of  its  transparency  (from 
1  to  0) ;  such  gradual  filters  can  be  fabricated  either  by  coating  a  piece  of  glass  or  by 
photographic  means  [8.20]. 
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Fig.  8.6:  Dark-field  schlieren  stop 

Probably  the  most  spectacular  flow  photographs  can  be  obtained  with  color  schlieren 
systems.  These  systems  require  the  use  of  a  white  light  source.  Light  deflections  in  the 
test  field  are  converted  into  colors  by  respective  elements  in  the  focal  plane  of  the 
schlieren  head.  A  few  of  such  elements  are  listed  in  Table  8.1,  together  with  represen¬ 
tative  references.  Additional  possibilities  for  introducing  color  in  a  schlieren  system 
may  be  found  in  the  reviews  of  Schardin  [8.3]  and  Settles  [8.8].  Colors  allow  for  ea¬ 
sily  discriminating  between  positive  and  negative  deflection  angles  (e.g.  with  the  three- 
color-filter)  ,  or  for  recognizing  the  direction  of  light  deflection  in  the  x-y  plane 
(with  the  multi-color-filter) .  With  a  dispersion  prism  or  a  diffraction  grating  one  is 
restricted  to  the  natural  sequence  of  colors  in  the  visible  spectrum. 

Table  8.1 


Color  schlieren  system 


optical  element 

providing  color 

visualization  of 

reference 

tri-color  filter 

Kessler  and  Hill  [8.21] 

dispersion  prism 

3p 

3y 

Holder  and  North  [8.22] 

diffraction  grating 

Maddox  and  Binder  [8.23] 

multi-directional 

color-filter 

4#  and  |£ 

Settles  [8.24] 

rainbow  filter  Howes  [8.25] 

("bull's  eye  filter") 


A  number  of  modifications  of  the  schlieren  system  are  known,  in  which  a  certain 
portion  of  the  light  source  image  in  the  focal  plane  of  the  schlieren  head  is  shifted  in 
phase.  A  shift  of  one-half  of  this  light  by  an  angle  of  18o°,  realized  by  a  respective 
phase  plate,  results  in  a  dark  field  of  view  due  to  the  interference  of  the  light  having 
passed  through  the  phase  plate  with  the  unshifted  light  [8.26  -  8.28].  The  purpose  of  ma¬ 
nipulating  the  phase  in  the  focal  point  of  the  schlieren  head  is,  in  most  cases,  an  in¬ 
crease  of  the  optical  sensitivity  at  low  density  flow  conditions.  It  is  not  always  clear 
whether  one  should  regard  such  systems  as  a  schlieren  system  or  an  interferometer;  see, 
e.g.,  [8.29,  8.30]. 

The  latter  is  particularly  true  for  the  phase  contrast  method.  In  the  focal  point 
of  the  schlieren  head,  a  minute  phase  plate  that  coincides  with  the  zeroth-order  diffrac¬ 
tive  image  of  the  source  changes  the  phase  of  the  zeroth-order  light  by  an  angle  of  9o°. 
This  technique,  originally  developed  for  being  used  in  microscopy  of  phase  objects,  re- 
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solves  small  changes  of  the  optical  phase  with  much  higher  contrast  than  the  Toepler 
schlieren  method  or  any  interferometer  [8.10,  8.13,  8.31-8.33].  With  small  optical  phase 
change  or  at  a  low  density  level,  a  record  taken  with  the  phase  contrast  method  has  the 
appearance  of  a  schlieren  record  (Fig.  8.7).  At  an  elevated  density  level,  fringes  appear 
in  the  field  of  view,  and  the  record  is  equivalent  to  a  reference  beam  interf erogram,  the 
fringes  being  curves  of  constant  density  when  the  flow  is  two-dimenstional . 


Fig.  8.7:  Phase-contrast  recording  of  the  supersonic  flow  (M  =  3.74)  over  a  two-dimen¬ 
sional  profile  (M.  Philbert,  ONERA,  France  [ 8 . 1 o] ) . 


Fig.  8.8:  Schlieren  system  in  combination  with  holographic  recording,  (a):  Recording 
of  the  hologram  ;  (b)  :  reconstruction  of  the  hologram  and  visualization  in  a 
schlieren  system,  (see  also  [8.34]). 
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The  schlieren  method,  as  any  other  optical  visualization  system,  can  be  combined 
with  a  holographic  set-up.  If  one  takes  a  hologram  through  the  test  field,  one  may  freeze 
the  instantaneous  information  of  the  refractive  index  field.  When  the  test  field  is  re¬ 
constructed  from  the  hologram,  it  may  be  observed  with  the  suitable  optical  method.  Fig. 
8.8  shows  the  principal  arrangement  for  taking  such  a  hologram,  and  for  observing  the  re¬ 
constructed  test  field  with  a  schlieren  system  [8.34],  Such  a  postponed  visualization  is 
of  great  interest  when  the  test  flow  is  unsteady  and  does  not  allow  for  adjusting  the  op¬ 
tical  system  during  the  exeriment.  The  reconstruction  provides  practically  unlimited  time 
for  choosing  the  appropriate  method  and  the  most  suitable  adjustment  [8.35] . 
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9.  M0IR£  METHODS  AND  SPECKLE  PHOTOGRAPHY 


Many  attempts  have  been  made  in  modifying  a  schlieren  system  such  that  the  form  of 
the  signal  is  more  suitable  for  a  quantitative  evaluation.  A  common  feature  of  all  these 
attempts  is  that  one  generates,  in  the  undisturbed  field  of  view,  a  more  or  less  regular 
reference  pattern,  which  is  then  disturbed  in  the  presence  of  the  refract ively  disturbing 
flow.  A  comparison  of  the  disturbed  pattern  with  the  undisturbed  reference  pattern  allows 
for  measuring  and  mapping  the  deflection  angles  in  the  whole  field  of  view  ("moire  pat¬ 
tern")  .  The  pattern  most  frequently  used  is  a  system  of  parallel,  equidistant  fringes  pro¬ 
duced,  e.g.,  by  inserting  a  grid  or  a  system  of  equally  spaced  wires  in  the  path  of  the 
light  beam. 

Holder  andNorth  [8.7]  have  investigated  the  case  when  a  grid  is  placed  between  the 
schlieren  head  and  its  focal  point  (Fig.  9.1) .  An  imaging  lens  in  or  near  the  focal  point 


object  schlieren 

plane  head 


image 

plane 

(film) 


Fig.  9.1.:  Deflection  mapping  system  with  a  grid. 


focuses  a  plane  in  the  flow  onto  the  recording  plane.  If  a  density  variation  is  present  in 
the  flow,  the  fringes  appearing  in  the  recording  plane  are  not  straight  but  distorted,  the 
distortion  (fringe  shift)  being  proportional  to  the  amount  of  light  deflection.  Such  a 
fringe  pattern  is  equivalent  to  the  pattern  observed  with  the  "finite  fringe  width  align¬ 
ment"  of  a  shearing  or  schlieren  interferometer  (see  chapter  1o) .  If  the  fringe  distortion 
is  measured  normal  to  the  direction  of  the  undisturbed  fringes,  in  Fig.  9.1  assumed  to  be 
parallel  to  the  x-axis,  the  absolute  fringe  shift  As^  is  related  to  the  component  e  of 
the  deflection  angle  by  (see  [8.7]) 
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As 


y 


n-fs)  f- } 

s 


(9.1) 


The  notation  of  this  equation  follows  from  Fig.  9.1.  The  fringe  shift  does  not  de¬ 
pend  on  the  absolute  spacing  of  the  grid.  For  g  +  0  the  system  approaches  the  Toepler 
schlieren  set  up. 


Systems  like  that  in  Fig.  9.1  have  been  applied  to  flow  situations  with  not  too 
small  density  variations,  e.g.  flames  (an  example  is  given  in  Fig.  9.2) ,  shock  waves 
[9.1],  or  stratified  flow  in  salt  water  [9.2].  Dewey  [9.3]  has  demonstrated  the  usefull- 
ness  of  this  method  for  observing  large-scale,  openair  shock  configurations  (Fig.  9.3). 


Fig.  9.3:  Shock  waves  produced  by  the  explosion  of  two  5oo  kg  charges.  The  shock  wave 

is  seen  against  a  striped  background  (J.M.  Dewey  and  D.J.  McMillin,  University 
of  Victoria,  B.C.,  Canada;  see  also  [9.3]. 

Moire  patterns  generated  by  the  superposition  of  two  fine  grids  of  equal  spatial  frequen¬ 
cy,  but  rotated  at  a  certain  angle,  have  been  described  and  extensively  discussed  by  Kafri 
and  his  co-workers;  see,  e.g.,  [9. 4-9. 7].  These  systems,  which  they  used  mainly  for  the 

study  of  flames,  present  additional  degrees  of  freedom  in  alignment  in  comparison  to  the 
set-up  with  only  one  simple  grid  (Fig.  9.1) . 
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The  moire  pattern  is  a  large-scale  pattern,  in  contrast  to  the  micro-sized  pattern 
produced  by  a  method  based  on  speckle  photography.  The  principles  of  speckle  photography 
have  been  outlined  already  in  chapter  5.  Kopf  [9.8]  and  Debrus  et  al.  [9.9]  describe  the 
adoption  of  this  principle  to  the  measurement  of  refractive  index  fields. The  idea  is  that 
a  microsized  reference  pattern  is  produced  by  illuminating  and  imaging  a  plate  of  ground 
glass  on  the  recording  plane,  and  that  this  pattern  is  distorted  due  to  the  presence  of 
the  refractively  disturbing  flow  field  and  can  later  be  analyzed  by  the  respective  tech¬ 
niques  used  in  speckle  photography.  The  optical  systems  used  in  [9.8]  and  [9.9]  suffer 
from  a  possible  decorrelation  of  the  two  patterns  when  the  deflection  angles  of  the  trans¬ 
mitted  light  exceed  a  certain  value. 

An  optical  set-up  which  overcomes  this  defect  is  shown  in  Fig.  9.4,  see  [9.10]. 

The  principle  of  the  method  will  be  explained  with  this  system.  The  light  from  a  laser 
source  is  expanded  and ' transmitted  in  form  of  a  parallel  beam  through  the  flow  field.  A 
lens  focuses  a  plane  of  the  test  field  onto  the  ground  glass  plate.  A  second  imaging  lens 
focuses  a  plane  at  distance  1  from  the  ground  glass  onto  the  recording  plane  (photogra¬ 
phic  plate) .  Two  exposures  are  taken  on  the  same  plate:  First,  a  reference  exposure  in 
the  absence  of  the  flow;  and  a  second  exposure  with  the  flow  to  be  tested  in  place,  as 
shown  in  Fig.  9.4.  The  first  exposure  generates  on  the  photographic  plate  a  speckle  pat¬ 
tern  which  is  determined  by  the  structure  and  the  scattering  characteristics  of  the  ground 
glass.  This  pattern  is  distorted  in  the  second  exposure  due  to  the  light  deflection  in  the 
test  flow.  The  displacement  of  each  individual  speckle  between  the  two  exposures  is  a  mea¬ 
sure  of  the  local  deflection  angle  e(x,y);  x,y  being  the  coordinates  in  the  recording  pla¬ 
ne  (field  of  view) .  Due  to  the  imaging  of  the  test  field  on  the  ground  glass,  a  light  ray, 
deflected  by  an  angle  s  in  the  test  flow,  will  arrive  at  the  ground  glass  at  the  same 
point  where  the  corresponding,  undisturbed  ray  of  the  reference  exposure  has  arrived  (see 
Fig.  9.4) .  However,  the  two  rays  form  an  angle  e  between  each  other,  if  a  1 : 1  imaging  ra¬ 
tio  is  chosen.  In  the  plane  at  distance  1  from  the  ground  glass,  the  two  rays  are  sepa¬ 
rated  by  ds  -  e  •  1.  At  the  same  time,  ds  is  the  displacement  of  the  recorded  speckles 
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laserlight 


flow  field 


lens 


ground  lens  recording 

glass  plane  (film) 


Fig.  9.4:  Deflection  mapping  by  speckle  photography.  Speckles  are  produced  by  the 

ground  glass.  The  plane  at  distance  1  from  the  ground  glass  is  imaged  onto 
the  recording  plane  (see  [9.10]). 


in  the  double  exposure,  and  the  distribution  ds(x,y)  can  be  determined  with  the  point- 
by-point  reconstruction  technique  via  an  evaluation  of  the  respective  systems  of  Young's 
fringes  (see  chapter  5.2  and  Fig.  5.6).  The  system  of  Fig.  9.4  is  equivalent  to  the 
speckle  photographic  measurement  of  inclination  angles  [9.11].  The  distance  1  has  to  be 
chosen  such  that  the  recorded  speckle  patternsare  always  correlated,  i.e.,  1  must  be  de¬ 
creased  if  the  deflection  angle  z  increases,  and  it  is  thus  possible  to  investigate  exten¬ 
ded  flow  fields  producing  laxge  values  of  £/  which  was  not  possible  with  the  original 
arrangements  used  by  Kopf  [9.8]  and  Debrus  et  al.  [9.9]. 
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In  contrast  to  the  schlieren  method,  this  speckle  photographic  system  is  a  means 
for  quantitatively  measuring  the  distribution  of  the  deflection  angle  e(x,y).  In  one  ex¬ 
periment,  i.e.  in  one  double  exposure,  the  displacement  ds  is  measured  as  a  vector,  and 
the  two  components  of  the  deflection  angle,  ex  and  e^,  are  therefore  measured  simultane¬ 
ously.  They  can  be  converted  into  the  respective  derivatives  of  the  refractive  index  or 
fluid  density,  see  Table  6.3.  The  quantitative  results  are  comparable  to  those  which  can 
be  obtained  from  a  shearing  or  schlieren  interf erogram  (see  chapter  1o) .  But  the  number 
of  data  points  obtainable  per  square  unit  of  the  speckle  record  is  much  higher  than  in 
the  case  of  an  interf erogram.  For,  in  the  latter  case  the  derivation  of  data  is  restric¬ 
ted  to  the  existence  of  an  interference  fringe,  whereas  the  micro-sized  speckle  pattern 
allows  for  a  quantitative  evaluation  in  practically  any  point  of  the  field  of  view.  The 
usefulness  of  the  speckle  method  has  been  proven  with  applications  to  flames  [9.12],  con¬ 
vective  heat  transfer  [9.13]  and  mass  transfer  [9.14]. 

While  the  latter  references  are  applications  to  laminar  flow  problems,  it  has  been 
demonstrated  by  Wernekinck  et  al .  [9.15],  that  the  method  even  can  be  used  for  resolving 

the  light  deflection  in  turbulent  flow  with  fluctuating  density.  The  system  must  then  be 
operated  with  a  pulsed  light  source,  e.g.  a  ruby  laser.  It  is  thus  possible  to  map  the 
randomly  distributed  deflection  angles  and  to  provide  quantitative  information  from  a 
turbulent  flow  field,  which,  in  most  cases,  is  not  possible  with  an  interferometer  due 
to  the  usual  blurring  of  the  interference  fringes.  In  addition  to  this  principal  availa¬ 
bility  of  quantitative  data,  a  very  dense  distribution  of  the  information,  i.e.  values 
of  the  deflection  angle  e(x,y) ,  can  be  obtained  from  the  micro-sized  speckle  pattern. 

An  automated  evaluation  system  as  mentioned  in  chapter  5.2  is  of  great  advantage.  It  is 
therefore  possible  to  analyze  the  plane  distribution  e(x,y)  by  means  of  an  appropriate 
statistics . 

Erbeck  [9.16]  (see  also  [9.17])  has  determined  the  correlation  function  of  the  de¬ 
flection  angles  e  measured  in  a  slightly  heated  wind  tunnel  flow  with  grid  turbulence. 

Such  a  flow  is  assumed  to  have  homogeneous,  isotropic  turbulence.  An  interesting  result 
is,  that  the  correlation  coefficient  R*  is  always  positive  when  the  components  of  e, 
which  are  correlated  with  eachother,  are  orthogonal  to  the  direction  of  correlation,  while 
the  correlation  coefficient  can  become  negative  when  two  components  being  parallel  to  that 
direction  are  correlated  (Fig.  9.5).  This  result  is  verified  with  a  visualization  of  this 


Spatial  correlation  coefficient  R*  of  the  deflection  angle  e  of  light  passing 
through  a  turbulent  flow  with  temperature  fluctuations.  The  correlation  compo¬ 
nents  of  e  are  parallel  (|  |)  or  normal  (J_)  to  the  direction  of  correlation. 

T  is  the  distance  of  correlation  (see  [9.16,  9.17]). 


Fig.  9.5: 
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flow  by  spatial  filtering  (see  below).  In  [9.16,  9.17]  it  is  shown,  that  the  earlier  idea 
of  Uberoi  and  Kovasznay  [7.6]  for  providing  statistical  information  on  the  three-dimen¬ 
sional  density  field  can  be  realized.  A  mathematical  algorithm  and  the  assumption  of  ho¬ 
mogeneous,  isotropic  turbulence  allow  for  determining  the  spatial  correlation  of  the 
temperature  (or  density) ,  R*,  from  R*  (Fig.  9.6) .  This  result  is  a  true  spatial  correla¬ 
tion,  in  contrast  to  measurements  taken  with  a  cold  wire  probe  in  the  same  flow,  which 
need  the  application  of  Taylor's  hypothesis  for  the  spatial  correlation.  At  the  same  time 
this  result  is  a  clear  demonstration  of  the  advantages  of  an  optical  whole  field  method 
in  application  to  a  turbulent  flow. 


Fig.  9.6:  Spatial  correlation  coefficient  of  the  temperature  as  calculated  from  the  re¬ 
sult  shown  in  Fig.  9.5.  M  =  4o  mm  is  the  mesh  size  of  the  grid  used  for  gene¬ 
rating  grid  turbulence.  The  indicated  data  points  refer  to  various  probe  mea¬ 
surements  for which  it  was  necessary  to  assume  a  constant  mean  velocity  U  (see 
[9.16,  9.17])  . 

The  specklegram  taken  from  a  refractive  index  field  can  also  be  analyzed  with  the 
spatial  filtering  technique  (see  chapter  5.2  and  Fig.  5.1o).  The  result  compares  to  a 
schlieren  pattern  (Fig.  9.7).  The  reconstruction  of  the  aforementioned  flow  of  homogene¬ 
ous,  isotropic  turbulence  by  means  of  spatial  filtering  in  y-direction  (Fig.  9.8)  veri¬ 
fies  the  different  patterns  of  the  two  correlation  functions  in  Fig.  9.5:  The  figure 
shows  more  frequent  changes  in  sign  of  intensity  (dark/bright)  when  it  is  traced  in  ver¬ 
tical  direction  (symbol  ||  in  Fig.  9.5),  i.e.  possible  negative  values  of  R* ,  while  there 
are  only  few  changes  in  sign  when  the  figure  is  traced  in  horizontal  direction  (symbol  j) 
again  a  proof  of  the  possibility  of  directly  visualizing  statistical  characteristics  of 
a  turbulent  flow. 


Spatial  filtering  of  a 
specklegram  taken  of  a 
turbulent  jet  of  helium 
exhausting  vertically  into 
the  ambient  air. 


Fig.  9.7: 
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Fig.  9.8: 


Spatial  filtering  of  a 
specklegram  taken  of  an 
air  flow  with  fluctua¬ 
ting  temperature  and 
homogeneous,  isotropic 
turbulence  (grid  turbu¬ 
lence)  (K .Oberste-Lehn , 
Universitat  Essen,  F.R. 
Germany)  . 
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10.  WOLLASTON  PRISM  SCHLIEREN  INTERFEROMETER 


Instruments  belonging  to  the  class  of  shearing  interferometers  provide  the  same  kind 
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of  information  on  the  density  field  as  the  methods  discussed  in  the  previous  sections, 
i.e.,  the  signal  obtained  from  a  shearing  interf erogram  is  a  function  of  the  refractive 
index  gradient  or  the  density  gradient  in  the  test  flow.  The  description  of  a  system  is 
based  on  the  assumption  that  two  slightly  sheared  light  beams  are  transmitted  through  the 
refractive  index  field  (Fig.  lo.l).  The  optical  set-up  of  many  shearing  interferometers 
is  similar  to  that  of  a  schlieren  system?  therefore,  the  name  "schlieren  interferometer" 
has  become  popular.  Among  a  great  number  of  interferometric  systems  (see,  e.g.,  [lo.l, 
1o.2] ,  an  arrangement  with  a  Wollaston  prism  as  the  shearing  element  is  most  frequently 
used  in  experimental  fluid  mechanics,  gasdynamics,  and  aerodynamics.  Only  this  system  will 
be  discussed  in  this  chapter. 


Fig.  lo.l:  Principle  of  a  shearing  interferometer 

In  contrast  to  the  methods  described  in  chapter  8  and  9,  the  signal  obtained  with 
a  shearing  interferometer  is  caused  by  the  changes  of  the  optical  phase  of  light  waves 
passing  through  the  flow,  and  not  by  light  deflection.  Two  parallel,  laterally  sheared 
light  rays  passing  through  the  test  flow  at  positions  (y+~)  and  (y~)  are  considered 
(Fig.  lo.l).  The  "shearing  element"  positioned  in  or  near  the  focal  point  of  the  schlie¬ 
ren  head  provides  that  the  two  rays  coincide  after  they  have  passed  through  this  element. 
Under  the  assumption  that  the  conditions  of  optical  coherence  are  fulfilled,  the  two 
coinciding  light  rays  can  interfere  with  eachother.  If  the  refractive  indices  n  (or  the 
fluid  densities  p)  at  the  coordinates  (y+^)  and  (y-^1)  in  the  flow  are  not  equal,  the  two 
light  rays  (or  waves)  exhibit  a  difference  in  their  optical  path,  Al,  given  by 


A  l(x,y)  =  /  n(x,y+"2,z)  dz  -  /  n(x,y-^,z)  dz  .  (lo.l) 

Z1  Z1 

This  equation  is  derived  under  two  assumptions:  First,  the  coordinates  where  the  rays  en¬ 
ter  and  leave  the  flow;  z^  and  z 2/  are  equal  for  the  two  rays?  i.e.,  the  flow  is  bounded 
by  plane  walls  (deviations  from  this  situation  are  easy  to  handle) .  Second  and  more  im¬ 
portant:  The  propagation  of  light  in  the  flow  is  straight  and  undeflected.  This  assump¬ 
tion  is  reasonable  for  most  gas  flows,  but  it  might  be  violated  in  liquid  flows  and  in 
gases  with  very  strong  density  gradients  (e.g.  in  a  plasma  flow) .  Exceptions  from  the 
validity  of  this  assumption,  known  as  "strong  refraction  effects",  will  not  be  discussed 
here.  The  situation  then  is  much  more  complicated,  because  not  only  the  refractive  index 
or  the  density,  but  also  the  path  of  the  light  in  the  flow  is  unknown. 

In  a  shearing  interferometer,  the  shear  or  the  separation  of  the  two  considered 
light  rays,  d,  is  assumed  to  be  very  small,  e.g.  in  comparison  to  the  diameter  of  the 
field  of  view.  The  right  hand  side  of  eq.  (lo.l)  can  be  developed  into  a  Taylor  series 
around  n(x,y,z) .  Taking  into  account  only  the  linear  term  (the  quadratic  term  vanishes) , 
because  d  is  small,  yields 
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A l(x,y)  =  d 


8y' 


:n  ( x ,  y ,  z )  dz 


(1o. 2) 


The  quantitiy  A]_(x,y)  is  measured  in  the  recording  plane  in  form  of  interference  fringes. 
Bright  fringes  appear  where 


Y  •  A  l(x,y)  =  0,  +  1,  t  2, 


(1o. 3) 


with  X  being  the  light  wave  length.  The  equation  of  the  dark  fringes  is 


X  •  A 1 (x ,y)  =  + 


(1o.4) 


It  is  seen  that  the  measured  signal,  the  left  hand  side  of  eqs.  (1o. 3)  and  (1o. 4) ,  is 
two-dimensional,  whereas  the  quantity  to  be  determined,  n(x,y,z) ,  is  three-dimensional  in 
the  general  case.  A  direct  measurement  of  n  is  possible  for  a  plane  refractive  index 
field,  n(x,y).  Methods  for  deriving  the  three-dimensional  distribution  of  n  (or  p)  from 
the  plane  information  ( interf erogram,  specklegram)  will  be  discussed  in  chapter  12. 


The  dependence  of  the  measured  signal  on  the  gradient  of  the  refractive  index,  like 
in  the  methods  which  rely  on  the  measurement  of  the  deflection  angle  e,  is  obvious  from 
eq.  (10.2).  The  measurement  of  3n/3x  is  possible  by  rotating  the  optical  system  by  9o° 
around  the  optical  axis  (z-axis) . 

From  eq.  do. 2)  it  follows  that  an  interf erogram  with  3n/3y  =  const  will  appear  uni 
formly  illuminated;  all  light  waves  arrive  in  the  recording  plane  with  the  same  optical 
phase.  This  case  is  called  the  "infinite  fringe  width"  alignment  of  the  interferometer. 

It  can  be  aligned  in  a  different  way,  so  that  a  system  of  equidistant,  parallel  interfe¬ 
rence  fringes  appears  for  the  object  with  3n/3y  =  const.  This  "finite  fringe  width"  align 
ment  is,  in  principle,  the  superposition  of  an  artificial  object  with  (3n/3y)  =  a  *  x  + 
b  *  y  (a,b  determine  the  direction  of  these  fringes  in  the  x-y  plane) .  A  density  varia¬ 
tion  in  the  flow  will  distort  this  regular  fringe  pattern,  and  the  deviation  or  shift  of 
a  fringe  from  its  undisturbed  position  is  a  measure  of  the  density  disturbance.  A  fringe 
shift  by  one  fringe  width  is  equivalent  to  a  difference  in  optical  path  length  by  one 
wavelength.  It  follows  that  a  fringe  shift  AS  in  a  point  (x,y) ,  measured  in  terms  of  the 
fringe  width  S  of  the  undisturbed  equidistant  pattern,  is 

z2 

1  •  AS  (x , y)  =  y  f  -|~n(x,y,z)  dz  *  (1o. 5) 

Z1 

The  fringe  shift  AS(x,y)  is  now  the  signal  to  be  measured  in  the  interf erogram.  The  cur¬ 
ves  connecting  the  points  in  which  the  disturbed  fringes  intersect  with  the  undisturbed 
fringes  of  the  finite  fringe  width  case  are  identical  with  interference  fringes  of  the 
infinite  fringe  width  alignment. 


Which  of  the  two  alignments  is  preferable  depends  on  the  particular  test  situation. 
For  a  comparison  of  the  two  cases.  Fig.  1o.2  to  1o.7  show  the  same  test  object,  the  lami¬ 
nar  plume  raising  from  a  candle  flame,  in  the  infinite  and  the  finite  fringe  width  mode. 
This  flow  has  a  much  stronger  density  (or  temperature)  gradient  in  horizontal  (radial) 
direction  than  in  vertical  (axial)  direction.  From  these  figures  it  becomes  obvious  how 
the  different  alignments  of  the  interferometer  in  terms  of  rotation  of  the  shearing  ele¬ 
ment  (Wollaston  prism)  react  to  these  gradients.  In  the  finite  fringe  width  alignment, 
the  density  gradient  is  always  measured  in  a  direction  perpendicular  to  the  undisturbed, 
parallel  fringes.  The  Wollaston  prism  interferometer  was  originally  developed  for  use  in 
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Fig.  1o.2:  Schlieren  interf erogram  of  a  Fig.  1o.3:  Same  as  Fig.  1o.2;  infinite 

laminar  candle  flame;  fringes  fringe  width, 

oblique  finite  fringe  width; 


Fig.  1o.4:  Laminar  candle  flame;  fringes 

horizontal,  finite  fringe  width. 


Fig.  1o . 5 :  Same  as  Fig.  1o.4;  infinite 
fringe  width. 
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Fig.  1o.6:  Laminar  candle  flame;  fringes  Fig.  1o.7.:  Same  as  Fig.  I0.6,  infinite 
vertical,  finite  fringe  width.  fringe  width. 

interference  microscopes  [ 1 o . 3 ] .  Various  authors  have  adopted  this  principle  for  the  vi¬ 
sualization  of  compressible  flows  and  have  given  extented  descriptions  of  the  system,  e. 
g.  [1o.4  -  lo.lo].  The  optical  set-up  is  very  similar  to  a  schlieren  system  (Fig.  I0.8). 


schlieren 

head 


cross  section 
of  f(ow  facility 


recording 
plane  (film) 


Fig.  I0.8:  Schlieren  interferometer  using  Wollaston  prism  as  the  shearing  element. 


The  knife  edge  known  from  the  Toepler  system  is  replaced  here  by  a  prism  unit  consisting 
of  the  Wollaston  prism  and  two  crossed  polarizers.  In  the  fundamental  or  infinite  fringe 
width  alignment,  the  center  of  the  Wollaston  prism  coincides  with  the  focal  point  of  the 
lens  or  spherical  mirror  that  was  previously  called  the  "schlieren  head".  The  role  of  the 
Wollaston  prism  (which  actually  consists  of  two  prism  halfs  glued  together)  as  the  shea¬ 
ring  element  is  explained  with  the  aid  of  Fig.  1o.9.  A  light  ray  entering  the  prism  is 
separated  into  two  components  forming  an  angle  8  with  eachother.  This  angle  is  deter¬ 
mined  by  the  characteristics  of  the  prism,  and,  within  the  accuracy  of  a  first-order  ana¬ 
lysis,  8  is  independent  of  the  direction  of  incidence.  The  two  components  leaving  the 
prism  are  linearly  polarized,  the  directions  of  polarization  being  perpendicular  to  each- 
other.  In  Fig.  1o.9,  the  incident  ray  1  is  separated  into  the  two  components  I1  and  1°, 
where  the  symbols  '  and  °  indicate  a  linear  polarization  in  the  plane  of  the  figure  and 
normal  to  this  plane,  respectively.  In  the  incident  converging  light  beam,  there  is  a 
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Fig.  1 o . 9 :  Reparation  of  incident  light  rays  in  the  Wollaston  prism.  The  symbols  1  and 
designate  linear  polarization  in  the  plane  and  normal  to  the  plane  of  the 
figure. 

ray  2  including  the  angle  3  with  1,  and  being  separated  into  2'  and  2°.  The  separated  ray 
2'  coincides  with  1°,  and  the  two  rays  could  interfere,  if  they  had  equal  polarization. 
This  requirement  is  fulfilled  by  means  of  the  second  polarizer  P2  of  the  prism  unit  (Fig. 
1o.8) ,  which  provides  a  polarization  under  45°  with  respect  to  the  directions  1  and  °. 

The  optical  axis  of  the  polarizer  is  either  parallel  or  rotated  at  9o°  to  the  axis  of 
P^/  and  its  purpose  is  to  provide  equal  intensity  to  the  two  interfering  components.  The 
same  procedure  applies  to  rays  1  and  3;  i.e.,  V  and  3°  can  interfere,  etc. 

The  distance  d  by  which  conjugate  rays  are  sheared  or  separated  in  the  test  section 
is  d  =  3  •  fg,  with  f  being  the  focal  length  of  the  "schlieren  head".  The  choice  of  the 
Wollaston  prism,  which  is  commercially  available,  usually  depends  on  the  desired  size  of 
the  separation  d.  In  many  practical  applications,  a  reasonable  size  range  of  shear  is 
1  mm  <  d  <  3  mm,  and  with  focal  lengths  of  the  schlieren  head  ranging  between  1  and  3  m, 
the  prism  should  provide  separation  angles  of  a  few  minutes  only.  A  number  of  different 
criteria  determine  the  choice  of  d.  This  quantity  should  be  small  enough  so  that  equa¬ 
tion  ( 1 o . 2)  applies.  On  the  other  hand,  if  the  density  gradient  to  be  resolved  decreases, 
the  shear  d  should  increase.  The  edges  of  rigid  bodies,  which  are  normal  to  the  direc¬ 
tion  of  shear  d,  appear  as  a  double  image  or  gray  zone  of  width  d  (or  d  •  cosa,  where  a 
is  the  angle  between  the  direction  of  shear  and  the  normal  to  the  wall) .  The  formation  of 
this  double  image  is  due  to  the  blocking  of  one  of  the  two  conjugate  rays  by  the  wall. 

The  exact  position  of  the  wall  edge  is  in  the  middle  of  the  gray  zone. 

If  the  prism  unit  is  displaced,  along  the  optical  axis,  by  a  distance  w  from  the  fo¬ 
cal  point  of  the  schlieren  head  (Fig.  1o.9),  parallel,  equidistant  fringes  will  appear  in 
the  field  of  view.  The  fringe  width  S  of  this  "finite  fringe  width"  mode  is  inversely  pro¬ 
portional  to  w: 

\.fs  __  X-f* 

S  3*w  d-w  (1o. 6) 

Rotation  of  the  prism  unit  around  the  optical  axis  rotates  the  parallel  fringe  system. 

As  it  was  mentioned  above,  the  prism  unit  can  be  implanted  into  an  existing  schlie¬ 
ren  system  in  place  of  the  knife  edge. This  allows  for  producing  interferograms  in  a  rela¬ 
tively  large  field  of  view  (e.g.,  0,5  m  in  diameter  as  in  [1o.7]) .  If  necessary,  use  can 
be  made  of  the  advantages  of  a  double-pass  system  (Fig.  8.5).  The  optical  set-up  shown 
in  Fig.  1o.8  requires  a  coherent  light  source,  i.e.,  a  laser  whose  beam  is  expanded  and 
spatially  filtered.  Since,  in  most  cases,  the  light  from  a  laser  has  linear  polarization, 
is  it  necessary  either  to  adjust  the  direction  of  this  polarization  with  the  axes  of  the 
polarizers  of  the  prism  unit,  or  to  pass  the  laser  beam  through  a  A/4-phase  plate  to  pro- 
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vide  circular  polarization  to  the  beam.  Many  schlieren  systems  are  operated  with  a  white 
source.  In  converting  such  a  system,  a  second  Wollaston  prism  (with  the  same  separation 
angle  3)  must  be  used  on  the  side  of  the  light  source  for  improving  the  spatial  coherence 
of  the  light  (Fig.  lo.lo).  The  fringesof  higher  order  now  appear  colored,  and  only  the 
fringe  of  zeroth  order  is  either  white,  if  the  axes  of  the  two  polarizers  are  parallel, 
or  black,  if  the  two  polarizers  are  rotated  at  9o°.  In  earlier  descriptions  of  the  inter¬ 
ferometer,  the  role  of  this  second  Wollaston  prism  was  interpreted  as  to  give  to  the  light 
beam  the  necessary  shear. 


Wollaston 

prfsnn 


Fig.  lo.lo:  Operation  of  the  schlieren  interferometer  with  two  Wollaston  prisms;  here: 

left  (light-source)  side  of  the  system.  The  right  side  is  as  in  Fig.  1o.8, 
but  without  the  polarizer  . 

A  great  number  of  compressible  flow  situations  have  been  made  visible  with  the  Wol¬ 
laston  prism  interferometer;  see,  e.g. ,  [1o.11]*  Of  particular  interest  is  the  visuali¬ 
zation  or  measurement  of  a  shock  wave.  Since  the  density  gradient  through  the  shock  is 
infinite,  the  linearized  equation  (1o.2)  cannot  describe  the  fringe  pattern  in  this  case. 
Instead,  the  complete  equation  ( 1 o . 1 )  must  be  considered.  Only  those  pairs  of  conjugate 
rays,  which  have  one  ray  passing  behind,  and  the  secondray  passing  in  front  of  the  shock 
wave  can  contribute  to  the  formation  of  the  respective  pattern  (Fig.  1o.11) .  The  infini¬ 
tesimally  thin  shock  surface  appears  with  finite  thickness  d  or  d  •  cosa,  in  analogy  to 
the  formation  of  a  double  image  of  the  edges  of  rigid  bodies  (see  above) .  In  the  finite 
fringe  width  mode,  the  shock  cannot  be  visualized  when  the  fringes  are  normal  to  the 
shock  (one  has  then  a  =  9o°l).  Also  it  is  not  appropriate  to  have  the  fringes  parallel 
to  the  shock  surface  (Fig.  lo.l) .  If  the  fringes  are  oriented  oblique  to  the  shock,  one 
measures  a  relative  fringe  shift,  AS/S,  which  is  given  by 


[as]  _  K  .z2 

'■S  shock  ^  z, 


Ap  dz 


(1o. 7) 


where  K  is  the  Gladstone-Dale  constant,  and  Ap  the  density  jump  across  the  shock.  This 
relationship  is  identical  with  that  derived  for  a  reference  beam  interferometer  (see 
chapter  11),  however,  the  respective  patterns  are  quite  different,  and  the  sensitivity 
of  the  two  methods  in  resolving  the  density  change  through  a  shock  wave,  therefore,  is 
not  the  same . 


eaocfc  position 
""  oi  shock  wav? 


Fig.  1o.11:  Visualization  of  a  (plane)  shock 
front  with  the  schlieren  inter¬ 
ferometer  aligned  at  finite  frin¬ 
ge  width?  above:  fringes  parallel 
to  shock  front;  below:  fringes 
oblique  to  shock  front. 
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The  Wollaston  prism  interferometer  is  particularly  attractive  for  heat  transfer 
studies,  because  the  heat  transfer  coefficient,  e.g.  in  natural  convective  flow,  is  pro¬ 
portional  to  the  temperature  gradient,  and  this  quantity  can  be  measured  directly  from 
the  fringe  shift,  provided  that  the  pressure  can  be  taken  as  constant.  The  principles 
of  this  interferometer  have  been  revised  several  times  under  the  particular  aspects  of 
determining  heat  transfer  rates,  [1o.12  -  1o.15].  The  numerous  reported  applications  re¬ 
fer,  e.g.,  to  the  measurement  of  convective  heat  transfer  [I0.I6  -1o.193,  mass  transfer 
in  free  jet  flows  [1o.2o,  1o.21],  plasma  flows  [1o.22] ,  or  stratified  flow  in  liquids 
[6.9,  1o.23] . 
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1  1  .  REFERENCE-BEAM- INTERFEROMETER 

a  reference  beam  interferometer/  a  light  wave  passing  through  the  test  flow  in— 
terferes  with  a  second,  undisturbed  wave,  the  reference  wave.  In  analyzing  the  principle 
of  this  interferometer/  onecan  make  use  of  eq.  (lo.l)  if  one  assumes  that  the  separation 
between  two  conjugate  rays,  d,  is  so  large,  that  one  of  the  two  rays  passes  outside  of  the 
test  section  in  a  uniform  environment  with  constant  refractive  index,  n^.  All  light  rays 
in  this  beam  {the  reference  beam) are  propagating  along  optical  paths  of  equal  lengths,  so 
that  eq.  (lo.l)  now  becomes 
z2 

Al(x,y)  =  f  [n (x,y , z)  -  n^]  dz  .  (11.1) 

Z1 

The  equations  for  the  bright  and  dark  interference  fringes,  eqs.  (1o.3),  (1o.4),  also  hold 
in  this  case,  and  it  becomes  evident,  that  the  measured  signal,  either  in  form  of  the 
fringe  order  or  a  fringe  shift,  is  directly  related  to  the  absolute  value  of  the  refrac¬ 
tive  index  n  (or  density  p) ;  see  Table  6.3. 

The  two  modes,  "infinite  fringe  width"  and  "finite  fringe  width"  are  known  also  for 
the  reference  beam  interferometer.  The  latter  situation  is  realized  by  the  superposition 
of  an  artificial  density  field  described  by  n  (x,y)  =  ax  +  by  (a,b,  =  const.). 

Among  the  many  possible  arrangements  of  a  reference  beam  interferometer,  only  two 
will  be  discussed  in  this  context.  An  important  criterion  for  the  application  of  such  an 
interferometer  to  aerodynamic  testing  is,  that  the  two  beams,  test  beam  and  reference 
beam,  be  widely  separated,  so  that  the  required  form  of  the  instrument  is  consistent 
with  the  geometry  of  the  test  facility. 

11.1  MACH- ZEHNDER- INTERFEROMETER 

This  instrument,  originally  designed  by  Ernst  Mach's  son  Ludwig  and  L.  Zehnder, 
combines  a  wide  separation  of  the  test  beam  and  the  reference  beam  with  a  relatively  lar¬ 
ge  field  of  view  (up  to  25  cm  in  diameter) .  Theory  and  practice  of  the  Mach-Zehnder-inter- 
ferometer  (MZI)  have  been  reviewed  several  times;  see  [11.1  -  11.4].  In  its  basic  arrange¬ 
ment  (Fig.  11.1)  ,  light  from  a  point  source  is  made  parallel  with  a  lens  or  spherical  (pa¬ 
rabolic)  mirror  L1 .  The  essential  components  of  the  interferometer  are  the  plane,  fully 
reflecting  mirrors  and  ,  and  the  plane,  semi-reflecting  mirrors  (beam  splitters)  M,j 
and  .  These  four  mirrors  are  arranged  so  that  their  centers  form  the  four  corners  of 
either  a  rectangle  (like  in  Fig.  11.1)  or  a  parallelogram.  In  the  latter  case,  a  larger 

test  section 


glass  plates 

Fig.  11.1:  Basic  arrangement  of  a  Mach-Zehnder  interferometer.  M  ,  M0  are  fully  reflec¬ 
ting  mirrors,  M^j  ,  beam  splitters.  1  2 

percentage  of  the  mirror  surface  is  effectively  used  in  the  optical  set-up;  or:  the  paral¬ 
lelogram  requires  smaller  mirrors  for  the  same  diameter  of  the  field  of  view.  The  test  sec- 
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tion  of  a  flow  facility  is  brought  in  the  path  of  the  test  beam.  If  this  test  section  is 
bounded  by  viewing  windows,  it  is  advisable  to  insert  two  identical  glass  plates  into  the 
path  of  the  reference  beam,  in  order  to  compensate  for  the  large  optical  path  length  tra¬ 
velled  by  the  test  beam  in  the  glass.  The  beam  splitter  M^  plays  the  role  for  the  "shearing 
element"  in  Fig.  lo.l,  providing,  however,  a  large  amount  of  shear  or  separation.  A  plane 
in  the  test  section  is  focused  onto  the  recording  plane  (film)  by  means  of  the  camera  ob¬ 
jective. 

Each  of  the  four  mirrors  is  required  to  permit  rotation  around  a  horizontal  and  a 
vertical  axis  for  the  basic  adjustment  of  the  interferometer,  in  which  all  mirrors  are 
parallel.  This  is  the  "infinite  fringe  width"  alignment,  i.e. ,  in  the  absence  of  a  test 
flow  and  with  all  optical  components  being  of  perfect  quality,  the  light  rays  arriving  in 
the  recording  plane  are  in  phase.  The  "finite  fringe  width"  mode  can  be  realized  by  til¬ 
ting  one  plate  or  mirror  by  a  small  angle  y.  The  resulting  width  of  the  parallel,  undis¬ 
turbed  fringe  system  is 


Certain  problems  arise  here  in  having  both  the  fringe  system  and  the  test  object  in  fo¬ 
cus.  Kinder  [11.2]  has  developed  an  arrangement  that  fulfills  this  requirement  and  can  be 
adjusted  by  controlling  only  one  mirror. 

ded  for  constructing  a  MZI.  Tolerances  of  surface  flatness,  surface  parallelism,  trans¬ 
lational  and  rotational  displacement,  control  of  tilt  must  be  in  the  order  of  one  tenth 
of  a  wavelength.  The  same  applies  to  the  optical  homogeneity  of  the  beam  splitters  and 
tie  windows of  the  test  section.  The  instrument  is,  therefore,  expensive,  the  costs  growing 
rapidly  with  the  desired  diameter  of  the  field  of  view.  The  mirrors  and  plates  are  moun¬ 
ted  on  a  U-shaped  steel  frame.  The  fourth  side  of  the  frame  is  open  for  inserting  the 
test  section  of  a  flow  facility.  The  required  mechanical  stability  of  the  frame  limits  the 
size  of  the  U  and  thereby  the  possible  separation  of  test  and  reference  beam.  It  will  be 
shown  below  that  many  of  these  requirements  of  mechanical  precision  do  not  exist  for  a 
holographic  interferometer,  which  has  replaced  the  MZI  in  many  fields  of  application. 

Modifications  of  the  MZI  have  been  described,  which  try  to  combine  a  large  diame¬ 
ter  of  the  test  beam,  i.e.  a  large  field  of  view,  with  a  reduction  of  the  reference 
beam's  dimensions,  in  order  to  minimize  the  mechanical  complexity  of  the  instrument.  The¬ 
se  modifications  require  the  use  of  a  coherent  laser  source.  In  the  arrangement  proposed 
by  Grigull  and  Rottenkolber  [11.5],  the  reference  beam  is  simply  the  thin,  collimated  la¬ 
ser  beam,  which  is  only  expanded  in  the  immediate  vicinity  of  the  second  beam  splitter 
^2  ^^-9*  11*2)  •  Howes  [11.6]  describes  a  set-up  which  is  basically  a  schlieren  system  of 


cross  section 
of  flow  facility 


Fig.  11.2:  Mach-Zehnder  interferometer  with  laser  light  source.  The  major  parts  of  the 
reference  beam  and  the  test  beam  are  fed  as  thin  laser  beams  (from  [11.5]) . 
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large  aperture  followed  by  a  relatively  small  MZI-like  arrangement.  The  reference  beam 
is  fed  through  a  spatial  filter,  where  all  the  phase  information  from  the  schlieren  field 
is  removed. 

It  is  difficult  to  oversee  the  large  number  of  reported  applications  of  Mach-Zehn- 
der  interferometry  to  the  study  of  compressible  flows,  and  a  list  of  references  would 
in  any  case  remain  incomplete.  Two  interferograms  are  included  here  as  typical  examples 
for  the  visualization  of  the  flow  pattern  in  a  supersonic  cascade  wind  tunnel  {Fig.  11.3) , 
and  the  diffraction  and  reflection  of  an  air  shock  wave  in  a  shock  tube  (Fig.  11.4)  .  In 


Fig.  11.3:  Mach-Zehnder  interferogram  of  the  supersonic  flow  in  a  cascade  wind  tunnel. 

Mach  number  of  incident  flow  is  Mra  =  1.24  (DFVLR-Institut  for  Antriebstechnik , 

Koln,  F.R.  Germany) . 

these  applications  there  is  a  particular  interest  in  determining  the  density  jump  across 
a  shock  front. These  quantity  is  directly  available  from  the  fringe  shift  measured  in  a 
finite  fringe  width  interferogram  and  the  application  of  eq.  (1o.7) ,  which  also  holds  for 
reference  beam  interferometry.  The  two  sample  figures  make  evident,  that  the  fringe  pat¬ 
tern  across  a  shock  front  is  different  here  from  the  pattern  obtained  with  a  schlieren 
interferometer  (see  Fig.  1o.12). 

Analyzing  a  turbulent  flow  with  density  fluctuations  by  means  of  interferometry  is 
difficult,  if  not  impossible,  because  the  interference  fringes  are  blurred  and  lose  their 
identity.  Flow  instabilities  can  be  observed  and  evaluated  up  to  a  certain  degree,  as  ve¬ 
rified  by  Anderson  et  al.  [11.7].  An  extended  literature  is  available  on  the  use  of  the 
MZI  for  convective  heat  transfer  studies  and  the  investigation  of  flames;  for  an  overview 
see  [6.1]  and  [6.2].  A  few  representative  applications  are  included  in  the  list  of  refe¬ 
rences  [11.8  -  11.12].  Error  sources  arising  from  the  special  situation  of  heat  transfer 
measurements  have  been  particularly  investigated  [11.3,  11.14].  A  modification  of  the  MZI 
has  been  used  for  measurements  in  a  spacelab  environment  [11.15]. 
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Fig.  11.4:  Mach-Zehnder  interferogram  of  the  diffraction  of  a  shock  wave.  Shock  wave 
propagation  is  from  left  to  right.  Compressible  vortices  separate  from  the 
edges  of  the  vertical  walls  (H.  Reichenbach,  Ernst-Mach-Institut ,  Freiburg, 
F.R.  Germany) 


11.2  HOLOGRAPHIC  INTERFEROMETER 


Holography  is  a  means  for  freezing  all  the  information  contained  in  a  light  wave, 
and  for  reconstructing  this  wave  at  any  desired,  later  instant  of  time.  In  contrast  to 
normal  photography,  which  gives  information  only  on  the  amplitude  pattern  of  a  wave,  ho¬ 
lography  also  provides  the  complete  phase  information  of  such  a  light  wave.  It  was,  the¬ 
refore,  a  consequent  step  to  make  use  of  this  principle  in  the  field  of  interferometry 
where  two  (or  more)  light  waves  are  made  to  interfere  with  oneanother.  The  result  is 
known  as  holographic  interferometry  [11.16],  and,  as  noted  already  in  the  previous  para¬ 
graph,  this  relatively  new  technique  is  taking  over  the  role  that  the  "classical"  Mach- 
Zehnder  interferometry  has  played  over  several  decades  of  testing  flows  with  varying  fluid 
density . 

A  holographic  interferometer  is  operated  with  a  laser  light  source  whose  beam  is 
separated  by  a  beam  splitter  into  two  beams.  The  principal  set-up  (Fig.  11.5)  shows  that 


cross  section  of 
flow  facility 


Fig.  11.5:  Principal  arrangement  for  taking  a  holographic  interferogram  with  a  double 
exposure. 
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one  of  the  two  beams,  after  beeing  expanded  and  then  collimated,  passes  through  the  test 
flow  and  arrives  at  a  holographic  plate  in  the  recording  plane.  This  beam  will  be  called 
the  signal  beam  or  holographic  test  beam  (not  to  be  confused  with  the  test  beam  of  a 
MZIi) .  The  second  beam  is  fed  outside  of  the  test  object  where  it  remains  undisturbed. 

This  beam  is  called  the  holographic  reference  beam.  It  is  also  expanded  and  arrives  at 
the  recording  plane  where  it  overlaps  with  the  signal  beam.' The  result  is  a  fine-scale 
interference  pattern,  the  "hologram",  which  allows  for  reconstructing  the  light  wave  that 
has  passed  through  the  test  flow,  i.e.,  the  signal  beam.  The  reconstruction  consists  of 
illuminating  the  developed  hologram  with  only  the  holographic  reference  beam,  i.e.,  the 
hologram,  after  being  developed,  is  returned  in  its  original  position,  and  the  signal 
beam  is  blocked  off. 

So  far,  the  system  does  not  yet  allow  for  seeing  a  flow,  because  the  reconstructed 
wave  as  well  as  the  original  signal  wave  are  not  visible  to  the  naked  eye.  Interferometry 
is  achieved  with  the  following  procedure:  A  first  hologram  is  produced  as  described  above 
with  no  flow  in  the  test  section.  The  developed  hologram  is  replaced  in  its  original  po¬ 
sition,  the  compressible  flow  is  turned  on,  and  the  hologram  illuminated  with  both  the 
signal  beam  and  the  holographic  reference  beam.  Due  to  the  existence  of  the  latter  beam, 
the  signal  wave  or  beam  of  the  first  exposure  is  reconstructed.  This  wave  can  interfere 
with  the  now  existing  signal  wave  that  carries  information  on  the  existing  refractive  in¬ 
dex  field.  It  is  obvious  that  the  first  exposure  is  equivalent  to  the  reference  beam  of 
a  MZI ,  while  the  signal  beam  with  the  flow  being  turned  on  is  equivalent  to  the  test  beam 
of  a  MZI.  The  result  is  a  reference  beam  pattern  in  the  infinite  fringe  mode,  which  be¬ 
comes  visible,  if  one  views  from  behind  through  the  holographic  plate. 

When  such  a  system  is  operated  as  it  has  been  described,  i.e.  with  the  first  holo¬ 
gram  developed  and  replaced  in  its  original  position,  it  is  possible  to  observe  the  in¬ 
terference  pattern  during  the  experiment.  It  is  said  that  the  interferometer  is  operated 
in  "real  time".  A  time-dependent  flow  could  then  be  recorded  with  a  movie  camera,  as  it 
can  be  done  with  a  MZI.  The  technical  problem  here  is  that  the  developed  hologram  has  to 
be  returned  exactly  in  its  original  position,  i.e.  with  a  precision  of  a  small  fraction 
of  a  wave  length.  Since  this  is  difficult  to  realize,  holographic  plate  holders  have  been 
developed  in  which  the  plate  can  be  developed  in  place,  i.e.  without  removing  it  from  the 
holder.  An  alternative  is  to  record  two  exposures  in  the  system,  one  without  flow  and  the 
second  with  flow,  on  the  same  holographic  plate  ("holographic  double  exposure") .  When  the 
plate  is  developed,  the  interference  pattern  becomes  visible.  However,  this  operation  does 
not  allow  for  a  direct  or  "real  time"  observation  of  the  flow. 

From  this  description  of  the  method,  the  dif ferencesbetween  holographic  and  Mach- 
Zehnder  interferometry  become  obvious:  In  the  MZI,  the  test  beam  and  the  reference  beam 
exist  simultaneously,  but  they  are  separated  in  space.  In  a  holographic  interferometer, 
the  test  beam  (or  wave)  and  the  reference  beam  (wave)  exist  at  different  instants  of  ti¬ 
me,  but  they  coincide  spatially.  This  spatial  coincidence  of  the  signal  beams  in  the  two 
holographic  exposures  (without  and  with  flow)  results  in  an  important  technical  advantage, 
which  makes  the  holographic  interferometer  superior  to  the  MZI:  The  optical  paths  of  the 
two  signal  beams  differ  only  due  to  the  difference  in  fluid  density  between  the  flow  and 
the  no-flow  situation.  Optical  disturbances  or  impurities  in  the  glass  of  test  section 
windows  or  in  any  other  optical  component  of  the  system  are  cancelled  out,  and  a  hologra¬ 
phic  interf erogram  can  be  taken  through  glass  windows  of  relatively  low  quality. 

Holographic  interferometry  was  first  applied  for  flow  visualization  and  the  mea¬ 
surement  of  density  fields  by  Hef linger  et  al.  [11.17]  and  Tanner  [11.18].  When  the  many 
advantages  over  Mach-Zehnder  interferometry  became  evident,  the  new  method  was  soon  picked 
up  by  a  number  of  laboratories  [11.19  -  11.21].  The  further  development  is  well  reflected 
in  the  book  of  Vest  [11.16].  Today,  holographic  interferometers  are  commercially  available, 
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but  an  engineer  who  is  not  totally  unfamiliar  with  the  practice  of  optics  can  well  build 
his  own  system.  A  practical  realization  of  such  a  system  is  shown  in  Fig.  11.6.  While  ma¬ 
ny  systems  are  built  flat  on  a  socalled  holographic  table  (a  heavy,  plane  granite  plate 
protecting  the  system  against  mechanical  vibrations) ,  this  system  is  mounted  vertical  on 
an  optical  bench.  The  beam  expander  included  in  the  set-up  of  Fig.  11.6  provides  an  ex¬ 
pansion  of  the  thin  laser  beam.  Such  a  unit  consist  of  an  objective  of  short  focal  length 
(microscope  objective)  and  a  small  pinhole  placed  in  the  focal  point  of  the  objective 
(Fig.  11.7).  The  diameter  of  this  pinhole  is  of  the  size  of  the  zeroth  diffraction  order 
of  the  light  in  the  focal  point,  and  it  serves  for  cleaning  the  expanding  beam  from  dif¬ 
fractive  noise. 


Fig.  11.6:  Set-up  of  a  holographic  interferometer;  M. :  mirrors?  S:  beam  splitter;  E: 

beam  expander;  L,:  lenses;  G:  glass  windows  of  test  section;  D:  diffuser 
glass;  H:  holographic  plate  (after  [11.22]). 


microscope 

objective 


Fig.  11.7:  Optical  element  for  expanding  and  spatial-filtering  of  a  laser  beam. 

The  signal  beam  in  the  arrangement  of  Fig.  11.6  is  passed  through  two  test  section 
windows.  The  mentioned  compensation  of  window  imperfections  in  the  holographic  double  ex¬ 
posure  is  achieved  only  when  the  state  of  the  glass  does  not  change  between  the  two  expo¬ 
sures.  A  deviation  from  this  situation  has  been  described  by  Koster  [11.14].  He  shows 
that,  in  performing  heat  transfer  measurements,  a  temperature  field  may  build  up  in  the 
test  section  windows,  thereby  causing  a  time  change  in  the  refractive  index  distribution 
in  the  windows.  A  similar  problem  may  arise  when  an  optical  method  is  applied  to  a  cryo¬ 
genic  wind  tunnel.  Snow  et  al.  [11.23]  report  on  the  observation  of  a  tremendous  loss  of 
signal  quality  when  optical  diagnostic  techniques  are  used  with  such  a  facility.  The  ex¬ 
planation  are  refractive  index  variations  induced  by  strong  temperature  gradients  in  the 
environment  of  the  (cool)  test  section  windows.  This  refractive  index  disturbance  outsi¬ 
de  of  the  test  section  would,  of  course,  not  be  compensated  for  in  a  holographic  double 
exposure. 

The  arrangement  shown  in  Fig.  11.6  includes  a  diffuser  plate  or  ground  glass  in 
the  path  of  the  signal  beam.  With  this  element  one  introduces  a  certain-  kind  of  spatial 
(three-dimensional)  information  in  the  interferogram.  Since  the  diffuser  scatters  the 
light  in  various  directions,  light  beams  are  traversing  the  test  field  in  a  variety  of 
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directions,  and  the  reconstructed  scene  may  be  examined  under  different  viewing  angles. 

In  order  to  observe  interference  fringes  it  is  necessary  that  the  optical  path  lengths  of 
all  rays  passing  with  different  direction  through  the  same  image  point  are  independent  of 
the  viewing  angle.  The  respective  formulation  of  this  condition  (see,  e.g.,  [  1 1 . 1 6]  , [  1  . 1  ]  ) 

defines  a  surface  in  the  test  field  onto  which  the  recording  plane  (camera)  has  to  be  fo¬ 
cused  for  obtaining  fringes.  If  these  requirements  are  fulfilled,  interf erograms  can  be 
taken  in  various  directions  from  the  same  holographic  double  exposure,  the  possible  vie¬ 
wing  directions  covering  an  angle  in  the  order  of  -2o°  <  a  <  +  2o° .  Problems  of  fringe 
clarity  arising  from  this  arrangement  have  been  discussed  by  Tanner  [11.24]. 

Several  technical  solutions  exist  for  producing  a  holographic  interferogram  in  the 
finite  fringe  width  mode.  When  the  system  is  operated  in  real  time,  fringes  of  finite 
fringe  width  will  appear  with  an  undisturbed  test  field,  if  the  developed  first  hologra¬ 
phic  exposure  is  not  replaced  precisely  in  its  original  position.  A  more  accurate  control 
of  fringe  spacing  and  direction  of  the  undisturbed  fringe  pattern  is  possible  with  the 
dual-plate  holder:  Reference  (first)  and  test  (second)  hologram  are  taken  on  two  separate 
plates,  which,  after  being  developed,  can  be  rotated  and  sheared  with  respect  to  each 
other  in  the  holder  during  the  simultaneous  reconstruction  [11.25,  11.26].  The  same  result 
is  obtained  by  a  small  change  in  the  adjustment  of  the  holographic  reference  beam  between 
the  two  exposures  [11.27]. 

Several  authors  describe  the  realization  of  the  finite  fringe  width  mode  by  utili¬ 
zing  in  the  system  two  separate  holographic  reference  beams  [11.28  -  11.31],  see  Fig. 

11.8.  One  reference  beam  serves  for  taking  the  reference  (first)  exposure,  the  second 
beam  is  used  for  the  test  exposure.  In  the  reconstruction  process,  the  hologram  is  illu¬ 
minated  by  the  two  holographic  reference  beams  simultaneously,  and  the  interferogram  ap¬ 
pears  in  the  desired  fringe  mode.  The  appropriate  fringe  pattern  is  selected  by  control¬ 
ling  the  angle  between  the  two  beams. 


flow 


mirrors 


laser 


Fig.  11.8:  Holographic  interferometer  with  two  reference  beams  (after  [  1 1  .  3  o] ) . 

Since  about  197o,  the  number  of  applications  of  holographic  interferometry  has  ex¬ 
ploded.  Many  laboratories  have  taken  over  the  new  technique  for  experiments  with  compres¬ 
sible  flows  in  wind  tunnels,  ballistic  ranges,  shock  tubes  (Fig.  11.9)  or  for  heat  trans¬ 
fer  studies,  and  laboratory  reports  describing  details  of  the  new  technique  have  been  pro¬ 
duced  ;  see,  e.g.,  [11.32  -  11.34].  Only  a  few  representative  references,  which  stand  for 

many  others,  can  be  given  [11.35  -  11.41].  The  resolution  of  three-dimensional  density 
fields  requires  taking  more  than  only  one  holographic  interferogram,  because  the  recon¬ 
struction  of  such  a  field  needs  a  total  viewing  angle  greater  than  the  mentioned  range 
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of  +  2o°  (see  chapter  12). 


Fig.  11.9;  Holographic  interferogram  of  a  focusing  shock  wave  taken  in  a  shock  tube. 

The  focusing  shock  wave  is  seen  in  the  center  in  form  of  a  small  triangle. 
Initial  shock  Mach  number  M  =  1.26  in  nitrogen  at  a  pressure  of  1  bar  (H. 
Gronig,  RWTH  Aachen,  F.R.  Germany). 
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12.  AXISYMMETRIC  AND  3-D  FLOW  FIELDS 


Due  to  the  integration  along  the  light  path,  the  signal  obtained  by  the  optical  vi¬ 
sualization  methods  from  a  three-dimensional  flow  is  plane  or  two-dimensional.  This  is  ob¬ 
vious,  e.g.  from  the  equations  (1o.2)  or  (11.1)  where  the  signal  Al  is  given  as  a  function 
of  x,y,  whereas  the  refractive  index  n  depends  on  the  three  Carthesian  coordinates  x,y,z. 
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In  order  to  resolve  the  dependence  of  n  on  z ,  it  is  necessary  to  direct  the  optical  light 
beam  in  different  directions  through  the  flow  ("tomography").  As  it  will  be  discussed  be¬ 
low,  the  number  of  necessary  viewing  directions  is  reduced  if  the  flow  has  some  kind  of 
symmetry.  In  the  case  of  rotational  symmetry,  one  projection  or  viewing  direction  is  suf¬ 
ficient  for  determining  n(x,r) ,  where  it  is  assumed  that  the  x-direction  is  the  axis  of 
rotational  symmetry. 


12.1  AXISYMMETRIC  FLOW 


The  light  propagates  again  in  z-direction,  and  the  refractive  index  depends  on  the 

1/2 

two  coordinates  x  and  r=(y2+z2)  ,  see  Fig.  12.1.  The  equations  relating  the  optical  sig¬ 

nal  with  the  refractive  index  distribution  must  be  transformed  accordingly.  These  equa¬ 
tions  have  on  their  left  side  a  signal  or  data  function,  D  (x,y) ,  with  x,y  being  the  co¬ 
ordinates  in  the  recording  plane.  D(x,y)  stands  for  either  the  fringe  order  or  the  rela¬ 
tive  fringe  shift  in  interferograms ,  or  the  light  deflection  angle  s(x,y)  as  measured 
with  speckle  photography  or  moire.  The  right  side  is  a  refractive  index  function  R(x,y,z) , 
whose  form  follows  from  the  respective  equations  in  chapters  8,  1o  and  11.  After  transfor- 


Fig.  12.1:  Cross-section  of  axisymmetric  flow  field  with  outer  radius  r^.  The  light 
propagates  in  the  z-direction. 

mat ion ,  the  fundamental  equation  for  the  axisymmtric  case  is 
r 

D ( x ,y )  =  /  R(x,r)  (r2-y2)  2  *  d(r2).  (12.1) 

r=y 

The  axisymmetric  flow  field  is  bounded  by  the  outer  radius  r^.  The  forms  of  D  and  R  for 
the  three  quantitative  methods,  schlieren  interferometry,  reference  beam  interferometry, 
and  speckle  photography,  are  given  in  Table  12.1. 


Eq.  (12.1)  is  an  integral  equation  of  the  classical  Abel  type  and  it  can  be  direct¬ 
ly  inverted,  yielding 


R(x,r) 


1  a 

7T  d  (r2  ) 


rco 

/  D(x,y) 
r 


_1 

(y2 -r2 )  2  d(y2)  . 


(12.2) 


Usually,  D(x,y)  is  not  available  as  a  continuous  function  but  in  form  of  a  set  of  discrete 
data  points.  Solignac  [12.1]  and  South  [12.2]  have  described  schemes  for  approximating  the 
set  of  discrete  data  by  appropriate  analytic  expressions.  Computer  programs  and  algorithms 
for  solving  the  Avel  inversion  (12.2)  have  been  presented  in  various  laboratory  reports, 
e.g.  [12.3,  12.4]. 
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Table  12.1 


Data  function  D  and  refractive  index  function  R  in  the  evaluation  of  axissymmetric  re¬ 
fractive  index  fields. 


reference  beam 

shearing 

speckle 

interferometer 

interferometer 

photography 

infinite 

D  _  Al(x,y) 

X 

r>  _  1  .  Al(x,y) 

y*d  X 

D  -  £ 

X  X 

fringe  width 

alignment 

R  =  j  {n(x,r)-noo} 

_  1  3n(x,r) 

R  ‘  X  "HPT ' 

D  =  e 

y  y 

finite 

D  _  AS(x,y) 

D  =  1  .  As(*>y> 

y  •  d  S 

R  =  in 

Rx  3x 

fringe  width 

alignment 

R  =  j  {n(x,r)-noo} 

„  _  1  .  8n(x,r) 

R  I  3 (r2 ) 

R  -  3n 

Ry  "5y 

A  disadvantage  of  the  Abel  inversion  is  the  needed  dif f erentiation  ,d/d  (r2 )  ,  of  the 
experimental  data  D(x,y),  as  shown  in  eq.  (12.2).  Such  a dif ferentiation  is  associated 
with  the  generation  of  errors.  The  disadvantage  is  avoided  with  optical  methods  depending 
on  the  measurement  of  the  refractive  index  gradient,  like  schlieren  method,  deflection 
mapping  by  speckle  photography  or  moire,  and  schlieren  interferometry.  The  two  differenti¬ 
ations  d/d(r2)  on  both  sides  of  the  respective  equations  are  then  eliminated  by  a  simple 
integration,  and  the  result,  written  for  the  case  of  def lectometry ,  where  the  deflection 
angel  e  is  measured,  becomes 


n (x, r)  -  n(x,rj 


x 

7T 


ey (x,y) 

y 


d(y2) 


(12.3) 


This  result  is  surprising,  because  eq.  (12.3)  states  that  the  methods,  whose  signal  de- 
lends  on  the  refractive  index  gradient,  deliver  the  absolute  value  of  the  refractive  in¬ 
dex  or  density  when  they  are  applied  to  an  axisymmetric  field  [1o.2o,  12.3,  12.5  -  12.7]. 


Eq.  (12.1)  and  Fig.  12.1  show  that  all  refractive  index  values  with  r>y  contribute 
to  the  variation  in  optical  phase  of  a  ray  propagating  at  position  y.  This  fact  has  been 
used  by  Schardin  [8.3]  for  developing  a  scheme  for  the  approximation  and  evaluation  of 
the  integral  in  eq.  (12.1).  At  the  same  time,  this  scheme  clearly  visualizes  that  one  pro¬ 
jection  is  sufficient  for  measuring  an  axisymmetric  refractive  index  field.  The  principle 
of  the  scheme  is  explained  with  Fig.  12.2.  A  circular  cross-section  of  the  flow  is  subdi¬ 
vided  into  N  annular  zones  of  width  h.  The  refractive  index  function  R  is  assumed  to  be 
constant  within  each  zone.  The  integral  is  solved  in  a  step-wise  procedure  starting  with 
a  ray  that  passes  through  the  most  outer  zone  (in  Fig.  12.2  ray  no.  1).  The  optical  signal 
generated  by  this  ray  is  sufficient  for  determining  the  respective  value  of  RJJ,  N  being 
the  number  of  the  most  outer  zone.  With  R^  being  determined,  ray  no.  2  in  the  next  step 
also  passes  through  only  one  zone  of  unknown  value  of  R,  namely  .  An  equation  can  be 

written  relating  the  optical  signal  of  ray  no.  2  to  the  known  contribution  of  RN  and  the 
contribution  of  R^  ,  which  can  now  be  determined.  The  procedure  continues  until  one  ar- 
rivesat  the  ray  through  the  center  of  the  circular  field.  Computer  algorithms  for  this 
procedure  have  been  described  in  laboratory  reports  [12.4,  12.8]. 

Of  particular  interest  for  supersonic  aerodynamic  studies  is  the  treatment  of  a 
discontinuity  in  the  refractive  index  distribution,  as  caused,  e.g.,  by  the  shock  around 
a  cone  without  angle  of  attack.  The  data  function  is  then  split  into  two  portions,  one 
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taking  account  of  the  continuous  density  field  and  the  second  of  the  density  discontinui¬ 
ty.  The  problem  has  been  discussed,  e.g. ,  by  Bennett  et  al.  [12.9],  Solignac  [12. 1o],  and 
Maruyama  et  al.  [12.11]. 


y 


2 


Fig.  12.2:  Subdivision  of  the  axi symmetric  flow  field  in  annular  zones  of  width  h. 


12.2.  TOMOGRAPHY 


The  problem  of  reconstructing  a  three-dimensional  object  field  from  a  set  of  two- 
dimensional  (plane)  projections  arises  in  various  areas  of  experimentation,  besides  flow 
visualization  also  in  X-ray  diagnostics,  electron  microscopy  or  radio  astronomy.  A  mathe¬ 
matical  formulation  of  the  problem  was  given  already  by  Radon  in  1917.  A  great  amount  of 
special  literature  of  this  aspect  is  available,  as  well  as  reviews  and  discussions  of  the 
general  aspects  [12.12  -  12.14].  As  in  the  axisymmetric  case,  there  are  ,  in  principle, 
two  ways  of  reconstructing  the  three-dimensional  density  field  without  symmetry:  inver¬ 
sion  schemes  and  discretization  of  the  problem. 

Projections  are  taken  in  a  number  of  different  directions  where  forms  the 
angle  ou  with  the  z-axis  (Fig.  12.3).  The  projection  generates  the  plane  data  field 
D(x,y,a.).  The  equation  to  be  inverted  is 

^i2 

D  (x,y ,a . )  =  /  R(x,y,z)dC.  ,  (12.4) 

where  and  £ ^  are  the  entrance  and  exit  coordinates  of  the  light  rays  at  the  test 

field  in  the  i-th  projection.  The  refractive  index  function  R(x,y,z)  must  be  presented 
by  series  of  generating  functions  which  allow  for  inverting  this  integral  equation.  In¬ 
version  schemes  with  appropriate  functional  series  have  been  described  by  Matulka  and 
Collins  [12.15]  and  Sweeney  and  Vest  [12.16].  Applications  of  these  schemes  include  the 
supersonic  flow  over  circular  cones  at  angle  of  attack  [12.17],  transonic  corner  flow 
[12.18],  mixing  in  a  buoyant  plume  [12.19],  or  the  transonic  flow  near  the  tips  of  a  ro¬ 
tor  blade  [12. 2o,  12.21].  The  referred  applications  are  all  based  on  holographic  inter¬ 
ferometry.  It  has  been  mentioned  in  chapter  11  that  the  spatial  information  obtainable 
from  one  holographic  interf erogram  is  not  sufficient  for  the  reconstruction  of  an  arbi¬ 
trary  3-D  density  field.  The  analysis  shows  that  it  is  necessary  to  use  projections  co¬ 
vering  an  angular  range  of  -9o°  £  ou  £  +  9o° . 

Inversion  schemes  which  are  based  on  Fourier  transforms  have  been  described  by  Row- 
ley  [12.22]  and  by  Junginger  and  van  Haeringen  [12.23].  Zien  et  al.  [12.24]  applied  this 
method  to  the  supersonic  flow  over  a  cone  at  angle  of  attack.  They  also  discuss  the  pro- 
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blem  that  arises  when  part  of  the  test  beam  is  blocked  off  by  the  model  in  the  wind  tun¬ 
nel,  i.e.  when  the  projections  are  incomplete ; see  also  [12.25]. 


Fig.  12.3:  Projection  of  light  in  ^-direction  through  a  three-dimensional  flow  field. 

Only  a  few  applications  of  discretisized  methods  for  reconstructing  a  three-dimen¬ 
sional  flow  field  have  been  reported  [12.26  -  12.28].  The  situation  is  represented  in  Fig. 
12.4,  and  this  figure  may  also  serve  for  providing  an  easy  description  of  the  general  pro¬ 
blems  of  reconstructing,  a  3-D  field  from  plane  projections.  The  figure  shows  a  cross-sec- 


Fig.  12.4:  Subdivision  of  a  cross-section  of  the  three-dimensional  flow  field  into  ele¬ 
ments,  in  each  of  which  the  refractive  index  function  is  assigned  a  constant 
value  R^.  Different  projections  of  the  light  through  the  field. 

tion  of  the  flow  in  the  y-z  plane  (see  also  Fig.  12.3)  ,  which  is  subdivided  into  N  seg¬ 
ments,  in  each  of  which  the  refractive  index  function  R  assumes  a  constant  value  R ^ ,  with 
1  £  j  £  N.  The  problem  is  to  determine  N  unknown  quantities  in  this  particular  cross-sec¬ 
tion,  R^  •  .  . . R^,  and  this  requires  defining  and  solving  a  set  of  N  independent  equations. 

A  light  ray  in  the  i-th  direction  carries  information  on  the  R^-values  of  m  segments,  with 
m  <  N.  This  information  is  expressed  by  a  respective  value  D^  of  the  data  function  D(x,y,a^) . 
The  equation  relating  D^  to  the  m  values  R^ , 

D.  =  l  R.s . . 
l  d 

is  just  one  equation  of  the  mentioned  set.  is  the  length  covered  by  the  light  ray 

within  the  j-th  segment  in  the  i-th  direction,  and  the  equation  is  seen  to  be  linear.  The 
light  beam  in  the  i-th  direction  cannot  provide  enough  equations,  and  it  becomes  necessary 
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to  direct  the  light  in  different  directions  through  the  refractive  index  field,  until  a 
sufficient  number  (N)  of  independent,  linear  equations  has  been  generated. 

The  number  of  necessary  viewing  directions  or  projections  increases,  the  finer  the 
subdivision  of  the  respective  cross-section,  i.e.  the  higher  the  number  N.  The  question 
then  is,  whether  the  visual  pattern  of  each  projection  is  sufficiently  different  from 
its  neighbouring  projection,  or:  whether  the  signal  is  sufficiently  sensitive  to  chan- 

gesof  the  directional  angle  cu  (Fig.  12.3).  Fig.  12.4  also  gives  an  indication  of  how  the 
situation  is  simplified  by  a  certain  degree  of  symmetry  of  the  flow:  Such  a  symmetry  simp¬ 
ly  reduces  the  number  of  unknowns  to  be  determined.  This  discussion  of  the  problem  with 
considering  only  one  cross-section  requires  that  the  light  rays  are  not  bent  out  of  the 
y-z  plane,  an  assumption  which  has  been  made  already  in  chapter  6. 
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13.  VISUALIZATION  OF  LOW-DENSITY  FLOWS 


The  optical  visualization  methods  discussed  in  chapters  7  to  1 1  have  a  lower  sensi¬ 
tivity  limit  in  resolving  changes  of  the  fluid  density.  It  is  somewhat  difficult  to  give 
general  number  values  for  these  limitations  in  sensitivity,  because  such  numbers  are  de¬ 
pendent  on  the  particular  flow  conditions.  There  is,  however,  a  class  of  fLuid  flows  whose 
density  level  is  in  any  case  below  the  range  of  sensitivity  of  the  optical  (refractive 
index)  methods.  These  flows  are,  in  most  cases,  governed  by  supersonic  or  hypersonic 
speeds  at  low  average  density,  sometimes  referred  to  as  “low-density  gas  dynamics".  Ex¬ 
perimental  means  for  visualizing  these  flows  have  been  developed  already  in  the  early 
stages  of  this  field.  These  methods  are  based  on  the  collision  of  electrons  with  the  mo¬ 
lecules  in  these  flows.  The  colliding  electrons  are  provided  either  by  an  electron  beam 
or  by  an  electric  discharge.  The  collision  causes  a  luminous  reaction,  with  the  intensi¬ 
ty  of  the  luminous  signal  depending  on  the  local  gas  density. 

13.1  ELECTRON  BEAM  FLOW  VISUALIZATION 


A  narrow  beam  of  high-energy  electrons  (as  generated,  e.g.,  in  a  TV  tube)  is  direc¬ 
ted  through  the  gas  flow  under  study.  Due  to  inelastic  collisions  between  the  fast  elec¬ 
trons  and  the  gas  molecules,  a  certain  portion  of  the  gas  molecules  are  excited.  Subse¬ 
quently  they  return  to  the  ground  state  with  emission  of  a  characteristic  radiation.  If 
this  returning  to  the  ground  state  is  prompt  ("spontaneous") ,  the  radiation  is  emitted 
more  or  less  at  the  same  position  where  the  gas  is  excited,  i.e.  at  the  position  of  the 
electron  beam,  which  appears,  therefore,  as  a  column  of  bright  fluorescent  light,  some¬ 
times  called  a  "fluorescence  probe".  Under  certain  conditions,  the  intensity  of  this 
direct  radiation  is  proportional  to  the  local  gas  density.  By  moving  the  electron  beam 
with  constant  speed  in  a  particular  plane  through  the  flow,  and  by  taking  a  photographic 
time  exposure  (with  open  shutter)  while  the  beam  is  moving,  one  obtains  a  representation 
to  the  local  distribution  of  the  gas  density  in  this  plane.  The  visible  signal  (Fig. 13.1) 
is  similar  to  the  pattern  of  a  schlieren  record,  however,  its  physical  information  is 
quite  different,  because  the  information  on  the  density  distribution  is  local  (one  ob¬ 
serves  the  scattered  radiation!)  and  not  integrated.  There  is  a  closer  analogy  to  the  me¬ 
thods  of  laser  induced  fluorescence  (chapter  5) . 

The  emission  of  the  visible  radiation  may  also  occur  from  an  excited  metastable  sta¬ 
te  whose  lifetime  is  relatively  long.  The  light  emission  then  takes  place  after  the  mole¬ 
cule  is  swept  a  certain  distance  with  the  flow,  i.e.  radiation  is  emitted  at  a  position 
downstream  of  the  original  beam  position  ("afterglow  radiation").  For  utilizing  this  af¬ 
terglow  radiation,  which  is  also  appropriate  for  visualizing  density  differences  in  the 
flow,  it  is  not  necessary  to  move  the  electron  beam  as  for  visualization  by  the  direct 
radiation,  but  the  intensity  of  the  afterglow  is  much  lower  than  that  of  the  former. 
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Fig.  13.1s  Electron  beam  visualization  of  a  supersonic  nitrogen  jet  (K.A.  Biitefisch, 

DFVLR  Gottingen,  F.R.  Germany) . 

(  Experiments  on  flow  visualization  by  electron  beams  have  been  reported  first  by 

Schumacher  in  1953  [13.1,  13.2].  It  became  a  frequently  used  technique  in  the  196o's 
when  hypersonic  flow  research  had  its  first  boom.  Many  details  of  the  method  have  been 
described  by  Lewy  [13.3].  In  addition  to  qualitative  visualization  of  the  flow,  the  elec¬ 
tron  beam,  in  combination  with  a  spectroscopic  analysis,  may  serve  for  quantitative  mea¬ 
surements  of  flow  temperature  and  density?  see,  e.g.  [13.4].  This  application  is  beyond 
the  scope  of  the  present  discussion. 

An  electron  beam  of  approximately  1  mm  in  diameter  is  generated  by  an  electron  gun 
which  is  operated  at  voltages  of  about  2o  kV  and  1  mA  current.  For  making  use  of  the  di¬ 
rect  radiation,  the  beam  is  either  moved  mechanically  parallel  to  itself  [13.5],  or  de¬ 
flected  over  a  certain  angular  section  by  means  of  deflection  coils  [13.3,  13.6].  When 
applied  to  the  flow  of  air  or  nitrogen,  the  direct  radiation  is  due  to  a  transition  of  the 
molecule  resulting  in  the  emission  of  light  at  a  wavelength  of  391,4  nm.  The  propor¬ 
tionality  between  the  observed  radiation  intensity  and  the  local  gas  density  is  only  the 
result  of  a  simplified  analysis,  in  which  a  number  of  real  effects  have  been  neglected. 

The  sensitivity  of  the  method  with  respect  to  resolving  density  variations  at  a  reduced 
level  of  the  average  density  is  about  one  order  of  magnitude  better  than  a  schlieren  sys¬ 
tem. 


The  afterglow  radiation  is  less  intense  than  the  direct  radiation.  It  can  be  ob¬ 
served  in  the  flow  of  nitrogen  and  noble  gases,  or  mixtures  of  both,  but  it  disappears 
almost  completely  in  air  due  to  collisional  quenching  between  excited  N2  and  nonexited 
02  molecules  [13.7]. 
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13.2  GLOW-DISCHARGE  FLOW  VISUAL I ZATION 


The  electric  discharge  at  several  hundred  volts  applied  between  two  electrodes  in 
a  low-density  gas  produces  a  plasma  column  in  the  gas  volume  between  the  two  electrodes. 
The  free  electrons  and  ions  in  the  plasma  are  accellerated  by  the  external  electric  field, 
and,  upon  collisions  with  the  neutral  gas  molecules,  they  can  excite  these  molecules, 
viiich  subsequently  emit  radiation  due  to  a  spontaneous  transition  into  the  ground  state. 
This  regime  of  visible  radiation  is  called  the  positive  column.  Within  a  certain  range  of 
the  gas  density,  usually  at  values  of  about  1o  ^  to  1o  ^  of  the  density  at  normal  condi¬ 
tions,  the  intensity  of  the  positive  column  is  a  direct  function  of  the  gas  density. 

The  visible  radiation  emitted  from  the  positive  column  has  been  used  for  visuali¬ 
zing  the  flow  in  low-density  wind-tunnels,  particularly  at  hypersonic  speeds.  The  poten¬ 
tials  of  this  technique  have  been  reviewed  by  Fisher  and  Bharathan  [13.8].  The  flow  in 
the  wind-tunnel  test  section  is  passed  through  the  volume  between  two  electrodes  to  which 
voltages  of  3oo  to  looo  V  ac  are  applied.  By  a  suitable  choice  of  the  geometry  of  the  e- 
lectrodes,  the  field  of  the  positive  column  can  be  varied  so  as  to  cover  the  desired  por¬ 
tion  of  the  flow  field.  In  some  cases,  the  test  model  itself  can  be  made  one  of  the  elec¬ 
trodes  . 

Density  changes  in  the  flow,  e.g.  produced  by  a  shock  wave,  appear  as  a  change  in 
intensity  and  sometimes  in  color  of  the  emitted  radiation  (Fig.  13.2)  .  The  technique  has 
been  applied  to  the  visualization  of  flows  of  nitrogen  [13.9],  helium  [13. 1o],  and  air 
[13.11].  The  information  on  the  density  distribution  is  not  local,  as  from  the  electron 
beam  method,  but  it  is  to  a  certain  degree  integrated,  because  the  visible  signal  origi¬ 
nates  from  a  particular  volume  in  the  flow  (the  positive  column) . 


Fig.  13.2:  Glow  discharge  visualization  of  the  flow  interference  between  two  free  jets 
expanding  into  a  vacuum  chamber;  chamber  pressure  is  0.18  mbar  (G.  Koppen- 
wallner,  DFVLR  Gottingen,  F.R.  Germany) . 
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14.  IMAGE  PROCESSING  AND  THE  GENERATION  OF  DIGITAL  FLOW  PICTURES 


In  this  context,  image  processing  is  the  evaluation  of  a  flow  picture  by  means  of 
a  computer.  For  this  purpose,  the  visual  record  obtained  with  one  of  the  visualization 
methods  must  be  digitized.  Such  a  record  is  a  plane  distribution  of  certain  visible  struc 
tures,  which  can  be  expressed  in  form  of  a  data  function  D(x,y) ,  with  x,y  being  the  coor¬ 
dinates  in  the  image  plane.  The  aim  of  image  processing  is  to  identify  these  structures, 
to  store  the  values  of  D(x,y)  in  the  memory  of  the  computer,  and  then  to  determine  the 
relevant  flow  parameters  from  the  data;  e.g.  the  flow  velocity  field  from  the  visualized 
time  line  pattern,  or  the  fluid  density  distribution  from  the  pattern  of  interference 
fringes.  The  visible  data  or  signal  is  often  superimposed  by  high-amplitude  noise.  Com¬ 
puter  programs  are  available  which  serve  to  separate  the  signal  from  the  noise,  i.e.  to 
recognize  the  desired  pattern  in  a  noisy  environment.  It  follows  that  image  processing  is 
not  a  visualization  method  but  a  means  for  making  an  evaluation  procedure  independent  of 
the  subjective  decisions  of  an  operator. 

An  abundance  of  literature  exists  on  the  general  aspects  of  image  processing  and 
on  the  development  of  special  computer  programs.  The  main  purpose  of  this  chapter  is  to 
give  a  number  of  relevant  references  which  describe  the  adaption  of  image  processing  me¬ 
thods  to  flow  visualization.  The  essential  step  for  this  adaption  is  the  digitization  of 
the  visual  record.  A  visualization  experiment  may  result  directly  in  a  digital  image  when 
an  appropriate  recording  system  is  applied.  In  the  simplest  case,  this  is  a  TV  camera, 
whose  output  is  digitized,  and  which  allows  for  taking  a  movie  picture  of  a  time-depen¬ 
dent  flow.  Winter  et  al.  [14.1]  discuss  the  use  of  vidicon-based  imaging  systems  and  sys¬ 
tems  of  photodiode  arays  for  flow  visualization  and  the  recognition  of  large-scale  struc¬ 
tures  in  the  flow.  Even  in  their  finest  form,  i.e.  with  the  maximum  of  pixels  now  avai¬ 
lable  per  square  unit,  these  systems  cannot  yet  reach  the  resolution  of  photographic  film 
(several  hundred  lines  per  mm) .  This  limits  their  application  to  the  recording  of  fine- 
scale  patterns,  e.g.  in  speckle  photography.  However,  it  is  expected  that  the  rapid  tech¬ 
nical  development  will  change  the  situation  in  the  near  future. 

When  the  flow  picture  is  digitized,  it  can  be  subjected  to  a  computer  analysis  in 
which  the  specific  structures  of  the  pattern  are  identified,  e.g.  the  interference  frin¬ 
ges,  time  lines,  or  the  images  of  tracer  particles.  The  recognition  and  identification 
usually  is  based  on  a  spectral  analysis  of  the  digital  pattern.  The  coordinates  of  the 
structures  or  of  their  edges  are  determined,  and  the  structures  may  be  approximated  by 
geometrical  curves,  which  finally  describe  the  data  distribution  D(x,y) . 


The  automatic  evaluation  of  interf erograms  has  been  described  for  the  various  fields 
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of  application  of  interferometry,  not  only  flow  visualization;  see,  e.g.  [14.2  -  14.4]. 

Kim  [14.5]  reports  on  the  approximation  of  the  fringe  contours  by  polynomial  fits.  Once 
such  curves  have  been  determined,  it  is  possible  to  interpolate  the  data  in  between  these 
curves  in  order  to  obtain  a  higher  density  of  the  data  per  square  unit  and  improve  the 
resolution  of  the  data  distribution  D(x,y).  Since  optical  interferometry  is  applied  often 
to  gas  dynamic  problems,  particular  attention  has  to  be  paid  to  the  abrupt  changes  of  a 
fringe  contour  as  caused  by  a  shock  wave  [14.6].  Among  all  applications  to  the  flow  visu¬ 
alization  records,  image  processing  is  best  developed  for  interferometry,  because  inter¬ 
ference  fringes,  particularly  those  of  a  laminar  flow,  make  a  well  defined  pattern. 

Time  lines  which  can  be  produced  by  various  visualization  techniques  (see  chapter 
4)  are  contours  that  can  be  subjected  to  an  analysis  similar  to  that  used  for  interfe¬ 
rence  fringes.  Image  processing  of  hydrogen-bubble  time  lines  recorded  with  a  video  ca¬ 
mera  has  been  described  by  Lu  and  Smith  [14.7]  and  Ongoren  et  al .  [14.8],  in  both  cases 
with  the  particular  emphasis  of  identifying  particular  structures  in  a  turbulent  flow. 
Nagib  et  al.  [14.9]  report  on  the  analysis  of  time  lines  generated  with  the  smoke-wire 
technique . 

The  streaks  generated  by  tracer  particles  in  a  photographic  time  exposure  from  a 
pattern  that  is  more  difficult  to  be  analyzed  by  a  computer  than  the  aforementioned  time 
lines  and  interference  fringes.  Therefore,  some  systems  for  the  processing  of  these  pat¬ 
terns  are  interactive,  i.e.  identification  of  the  streaks  is  performed  by  an  operator 
[14.1o] .  A  number  of  fully  automated  systems  and  procedures  have  been  described,  e.g. 
[14.11  -  14.13],  A  major  problem  in  the  analysis  of  the  streak  pattern  is  the  identifi¬ 
cation  of  the  flow  direction.  The  final  result  of  the  processing  is  the  display  of  a  field 
of  velocity  vectors  or  isovelocity  contours.  When  these  processing  methods  are  applied  to 
turbulent  flows,  it  is  of  interest  to  identify  large  scale  structures,  sometimes  called 
"coherent  structures",  which  are  masked  by  small-scale  turbulent  noise.  These  structures, 
characterized  by  low  wave  numbers  in  a  respective  Fourier  transform,  can  be  recognized 
when  the  digitized  image  is  filtered  with  a  2-D  low-pass  filter  function  [14.14].  The 
pattern  recognition  thus  is  made  independent  of  a  subjective  interpretation. 

Digital  flow  pictures  can  be  generated  easily  as  the  result  of  a  computation  or  a 
survey  of  the  flow  field  with  a  measuring  probe  .These  pictures  can  be  compared  with  the 
patterns  obtained  from  a  flow  visualization  experiment,  and  conclusions  my  be  drawn,  e.g., 
on  the  quality  of  a  computation.  Flow  visualization  by  means  of  computer  graphics,  how¬ 
ever,  is  beyond  the  definition  of  the  term  "flow  visualization"  in  this  review. 
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